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The Initial Stages of Growth of Oriented Copper Nuclei 
on Single Crystal Surfaces of Silver} 


By E. Grunpaumt, R. C. Newmans and D. W. PASHLEY|| 
Physics Department, Imperial College, London 


[Received July 12, 1958] 


ABSTRACT 


Layers of metallic copper of less than 1 A in average thickness have been 
deposited in a controlled manner on to smooth (111) surfaces of silver inside an 
electron diffraction camera. It is found that an electron diffraction pattern 
from oriented copper nuclei appears suddenly at some critical thickness. 
This critical thickness decreases with increasing substrate temperature, 
being about 0-9 A at room temperature and 0-2 4 at 300°c. Possibilities for 
the initial structure of the deposit, and for the change which occurs at the 
critical thickness, are considered, and the effect of substrate temperature is 
discussed. 


$1. LyTRoDUCTION 


IN a previous paper (Newman and Pashley 1955) results were given for 
the growth of extremely thin evaporated layers of copper on (111) faces 
of silver. It was found that, for an average deposit thickness of less than 
about 0-8 A, no visible change occurred in electron diffraction reflection 
patterns from the substrate surface. As the layer was increased in 
thickness beyond this value, a diffraction pattern due to oriented copper 
crystallites appeared, and rapidly increased in intensity. It was supposed 
that this initial appearance of a diffraction pattern from copper was 
controlled by a limiting sensitivity of detection. 

This paper reports the results of a more careful study of the rate of 
increase in intensity of the spots due to diffraction by copper, carried 
out in 1955. This work has shown that the pattern appears extremely 
rapidly, and that the change from no spots due to copper, to quite a 
strong visible pattern occurs for a thickness increase of only about 
0-054. This did not seem to be consistent with a limiting sensitivity, 
and it was considered that a more likely explanation is that a transition 
in the mode of growth occurs at a certain stage. In order to investigate 
the effect further, the following experiments were carried out: (1) the 
diffraction patterns have been studied over a range of electron energies; 
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(2) the effect of varying the substrate temperature has been investigated. 
The results of both of these experiments, reported below, support this 
latter view. 
§ 2, EXPERIMENTAL TECHNIQUE 

The copper was deposited on (111) faces of silver crystals in position 
in the electron diffraction camera; radioactive copper was used in order 
to allow accurate thickness measurements to be made. In the previous 
work (Newman and Pashley 1955) a large number of specimens of different 
thicknesses was prepared and examined, in order to investigate thickness 
effects. As this method is rather laborious, a method was developed to 
allow a series of measurements to be made on one single specimen which 
remains in situ in the diffraction camera. The automatic recording 
technique described by Kehoe et al. (1956) was used, and the rate of 
deposition was determined in the following manner (see fig. 1). 


Fig. 1 
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Radioactive copper (C) was deposited simultaneously on to the 
specimen (S) under examination, and on to a second similar specimen (‘T) 
which was one of a number mounted on the disc D. During the experiment, 
D was rotated at intervals so that each specimen on the disc received 
a copper deposit for a certain period during the deposition. At the end 
of the experiment, the thickness of the copper deposit on each specimen 
such as T was determined from its radioactivity, and a growth curve 
was plotted. A check on the rate of growth on the specimen (S) was 
made by means of a determination of the total thickness on this specimen. 
The technique was, in effect, used to check that the rate of growth was. 
constant, so that accurate interpolations could be made to obtain the 
thickness of the layer at any given instant at which a diffraction pattern 
was recorded. In all results reported below, the rate of growth was 
accurately constant during the entire deposition. 


§ 3. ExXpERIMENtTs at Vartous ELEcTRON ENERGIES 


Several runs were carried out at a number of different accelerating 
voltages in the range 10-60 kv. In each case, the mean thickness (hereafter 
called the critical thickness) at which the copper deposit suddenly appeared 
on the diffraction pattern was determined. This critical thickness was 
found to be the same in all cases, within experimental accuracy, and to be 


ees in reasonable agreement with previous results (Newman and Pashley 
00). si 
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§ 4. Experiments on HpatTep SuBstTRATES 


The deposition of the copper was studied on substrates which were 
maintained at temperatures ranging from 20-300°C, by means of the 
furnace described previously (Kehoe et al. 1956). It was found that 
on all specimens, the pattern due to oriented copper appeared suddenly 
at some critical thickness, which depended upon the substrate temperature 
as shown in fig. 2. As the substrate temperature is increased, the critical 
thickness becomes less, reaching a value of about 0-2 4 at 300°c. 


Fig. 2 
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§ 5, Discussion 


It seems unlikely that the sensitivity of detection of the copper deposits 
should be independent of the electron energy over the range 10-60 kv, 
since the inelastic and elastic scattering cross sections both vary 
considerably and differently over this range. This variation causes 
marked changes in the depth of penetration of the electrons, which 
would be expected to have a considerable effect on the proportion of the 
electrons diffracted by any surface film. The first experiment therefore 
gives strong support to the hypothesis that the critical thickness arises 
from some other cause than that of a limit of detection. Da 

The second experiment confirms this view, since the critical sepaat 
is found to depend upon temperature. It would be expected that ae 
limit of detection (in terms of average thickness/unit area) would e 
independent of the substrate temperature, provided the limit of Bt cies 
is not influenced by changes in the form of the nuclei. If the nuclei are 
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sufficiently small to allow complete penetration, without any appreciable 
absorption, the sensitivity of detection will have a maximum value. For 
nuclei which are sufficiently large to give rise to appreciable absorption 
effects, the larger the nuclei, the lower the sensitivity of detection. However, 
in the second experiment, the nuclei formed at the elevated temperatures 
are the larger, as indicated by the sharpness of the diffraction pattern, 
and yet they are detected at lower average thicknesses. This proves 
that the initial appearance of the diffraction patterns due to oriented 
copper nuclei is not related to a limiting sensitivity of detection, but arises 
from some property of the growth of the deposit. It is therefore necessary 
to consider (a) what form is taken by the initial deposits, such that they 
give rise to no observable diffraction patterns, and (6) what change takes 
place at the critical thicknesses indicated in fig. 2. 
There are two possibilities which would explain (a): 


(1) The initial deposit is in the form of small crystalline islands of 
monatomic thickness which have the same lateral spacing as that of the 
substrate silver. This is the mechanism for epitaxial growth given by 
Frank and Van der Merwe (1949 a, b) except that these authors consider 
the possibility of stability of the nuclei up to much greater thicknesses. 


(2) The initial deposit occurs as a distribution of single copper atoms 
which are able to diffuse over the silver surface. At any instant the 
position of a particular atom is determined by the substrate only and is 
not influenced by adjacent copper atoms. The copper atoms may thus be 
regarded as a two-dimensional gas (see Frenkel 1923). 


In both cases, it is reasonable that re-crystallization to give copper 
nuclei with their normal spacing should occur at some critical average 
surface density denoted by n,. The expected dependence of m», on 
temperature will now be examined for each case. 

In case (1) the pseudomorphic islands will re-crystallize when dislocations 
are generated in the islands. Since an activation energy is required for 
this process, re-crystallization is expected to occur more readily at higher 
temperatures. This temperature dependence will be enhanced by the 
growth of larger islands at higher temperatures due to the higher surface 
mobility of the depositing copper atoms; these larger islands would be 
inherently less stable than small islands as the probability of a dislocation 
being generated in a particular island is dependent on its size. The 
experimental observations would therefore be explained qualitatively by 
this mechanism. 

In case (2), it is assumed that the silver substrate exerts a periodic 
binding force on each copper atom. Ifthe energy barrier between adjacent 
sites is #, the number of mobile atoms at any instant will be given by the 
total number present on the surface, multiplied by the appropriate 
Boltzmann factor. Crystallization will occur when the local surface 
density exceeds a critical value such that the nucleus is formed and does 
not immediately re-evaporate. The frequency with which this occurs 
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will obviously depend on the mobility of the copper atoms. Thus, the 
value of n, when crystallization occurs will be given by 


n,=N, exp {+(E/kT)} 


where JN, is related to the critical number of mobile atoms required to 
forma nucleus. This mechanism would explain the observations provided 
N, does not vary too rapidly with temperature in opposition to the 
exponential factor. 

Unfortunately, it is not possible to calculate the activation energy for 
the generation of dislocations in pseudomorphic metal layers, as the force 
laws involved are not known. Also, the energy # above cannot be computed 
with any accuracy, due to the presence of adsorbed gases on the substrate 
surface during the deposition. It would thus appear that the experimental 
evidence is insufficient to allow any distinction to be made between the 
two possible mechanisms which have been suggested. 

The present results must lead to a modification of a conclusion concerning 
sensitivity which was reported previously (Newman and Pashley 1955). 
It was stated that the higher sensitivity of electron diffraction for detecting 
silver bromide (a minimum average layer of 0-4 A was detected) compared 
with copper deposits was due to differences in the scattering factors of the 
elements present in the layers. This is now contradicted by the present 
results where a copper deposit of only 0-2 A in mean thickness has been 
detected. It therefore seems likely that the first appearance of the silver 
bromide is controlled by the mechanism of growth rather than by the 
limiting sensitivity of detection. 
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On the Thermal Instability of a Rotating Fluid Sphere; 


By F. E. BissHopp 
Department of Physics and Yerkes Observatory, The University of Chicago 


{Received July 14, 1958] 


SUMMARY 


In this paper the problem of the thermal instability of a self-gravitating 
rotating fluid sphere of constant density with internal heating is treated. A 
method of solving this problem is described which correctly takes into account 
the full set of boundary conditions. The cases when the bounding sphere is 
rigid and when it is free are both treated. The Rayleigh—Ritz method of 
approximating lowest characteristic values is used in a variation principle 
presented by Chandrasekhar to calculate the critical Rayleigh number as a 
function of the Taylor number. 


§ 1. IyTRODUCTION 


In a recent paper Chandrasekhar has studied the problem of the onset of 
thermal instability in a rotating fluid sphere heated within (Chandrasekhar 
1957 a) and has determined the critical Rayleigh number as a function of 
the Taylor number for a specific set of boundary conditions which corre- 
sponds neither to a rigid boundary nor to a free one at the surface of the 
sphere. The aim of the present paper is to re-examine the same problem 
with the correct boundary conditions. 

In the formulation of the problem, axial symmetry of the perturbed 
motions and of the modified temperature distribution is assumed. The 
assumed symmetry allows the expression of the velocity perturbation as 
the superposition of a poloidal and a toroidal field in terms of two azimuth- 
independent scalars. As we shall see in § 2, the equations connecting 
these two scalars and the function which characterizes the temperature 
perturbation are a set of three partial differential equations, two of order 
two and one of order four. In these equations the Taylor number, T, can 
be considered a parameter to be specified in advance, in which case the 
equations together with the boundary conditions determine a characteristic 
value problem for the Rayleigh number, R. The critical Rayleigh number 
for the onset of convection is to be found by evaluating the lowest 
characteristic value. An entirely equivalent procedure which we will 
have no need of considering is to consider the Rayleigh: number as the 
parameter and the Taylor number as the characteristic value. 

When one is confronted with a problem of the type described above, one 
might, in general, expect to proceed in either of the following two ways: 
Upon obtaining the most general solution of the differential equations and 
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requiring that the solution satisfy the boundary conditions, one obtains 
an equation to determine the set of characteristic values; or, upon 
obtaining the most general set of functions which satisfy the boundary 
conditions and requiring that the solution satisfy the differential equations, 
one again obtains an equation to determine the set of characteristic values. 
Because of the complexity of the problem we are considering, it is necessary 
to employ a method which contains elements of both the above mentioned 
procedures. We shall first construct a set of functions which satisfy 
all the boundary conditions and the two second order differential equations. 
By requiring that these functions satisfy the remaining fourth order 
equation (or the equivalent variational equation) we shall determine 
the critical Rayleigh number as a function of the parameter T. 

In constructing the set of trial functions for the variational equation 
we will find it necessary to expand the defining scalar of the poloidal field 
in a complete set of functions which satisfy two boundary conditions at 
r=1 and which are bounded at r=0. A discussion of similar expansions 
can be found in a paper by Chandrasekhar and Reid (1957). 


§ 2. EQUATIONS OF THE PROBLEM 


In Chandrasekhar’s paper (1957 a) the equations governing marginal 
stability in the case where rigid rotation is present are presented and 
discussed adequately. We shall be content with merely presenting them 
here. 

Let us consider a homogeneous fluid sphere of radius ry, rotating with 
a constant angular velocity 9, about the z-axis. We assume there is 
within the sphere a uniform distribution of heat source which, in the absence 
of conduction, will cause the temperature at every point to rise at the 
rate €. In addition we assume no external forces acting on the sphere, 
we ignore rotational flattening, we assume axial symmetry of any 
perturbations, and we take into account the variation of density due to 
thermal expansion in the Boussinesque approximation. 

The equilibrium temperature distribution is given by (Chandrasekhar 
1952, 1957 a) 

T=B(re?—17); p=5, re es ee ee 
where x is the coefficient of thermal conductivity. 

The expression of the perturbed velocity field as the superposition 
of a poloidal and a toroidal field in terms of the two azimuth -independent 
scalars, U and V, can be written as follows (Chandrasekhar 1956 a, b, 


1957 a, Liist and Schliiter 1954) : 
u=I_xrV(r, w)+eurl [I,xrU(7,p)].  . . . (2) 


In terms of spherical coordinates (r, = cos #, ¢), 


1—p,?) 0 
u=— = [(il—-p2)0]I,- Yee! ay eG Wehr — p?) VI. (3) 
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The equations governing marginal stability assume the following 


ie QaVv 1 06 

2 as 

Ti ai OS eee) Et em ae 

ay El vy oz vrou (4) 

Ai 4 2 Eo eee 
Ve Oe 

IM Per phe [Salem Reet 168 
x Ou 


where 6 is the deviation from the equilibrium temperature distribution, 
v is the coefficient of kinematic viscosity, y= 4nGap/3, G is the gravitation 
constant, « is the coefficient of volume expansion, p is the density, 


@ 490 (1-p2)e 4ua@ 


ae ees —_—- == SR te me ha 8 
A; a ror Oe Oe 7) 
and 
Cen? 0.) (Lp) OF ee cise 
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As are rar re Out re Ope (8) 


are the five- and three-dimensional Laplacian operators for axisymmetric 
functions. 
By the transformation to the dimensionless variables, 


n=] f, 
Ue, w= Ur, w)|ro _ 
Vey = VV (ro eer eee 5 a 


O'(r’, wu) =O(r, 1)/2Bro* 2, | 
the equations assume the more convenient forms, 
oV 1 00 


Nee Tee eee Re a a Ue ee 
2 dz r Om’ eg) 
Av+ SS =0, elt. RT ee ee Ley 
and 
fa) 
et alice. (1 =p") U)], ©. ee ee eco) 
where 
20Y oe. 4 92 
Re es 16s ne eS) 


and the primes have been suppressed. The number denoted R is the 
Rayleigh number; T is the Taylor number. 
We now make use of the easily verified identity, 


Ly Cy seen elie 


to rewrite eqn. (12) as 
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The boundary conditions which must be satisfied by the solutions of 
eqns. (10), (11) and (15) depend upon whether the bounding surface at r=1 
is assumed to be a rigid or a free boundary. In both cases we require the 
vanishing of the radial component of the velocity and of the deviation of 
the temperature from the equilibrium value. At a rigid boundary we 
require the vanishing of w, and w, and at a free boundary we require the 
vanishing of the tangential viscous stresses. Accordingly, the boundary 
conditions for a rigid boundary at r=1 can be written; 


0U 00 
U=0 —_—_ = ae ae Viana . 
aor on eV Oe ee ee ee ee LO) 
and for a free boundary at r=1; 
0” 00 OV 
U=0 are = —_— = = 
» 5a (rU) 0, Bn 0, a Oe AL) 


§ 3. THE VARIATIONAL PRINCIPLE 


Chandrasekhar has shown that eqns. (10), (11) and (12) together with 
either of the sets of boundary conditions (16) and (17) imply the following 
formula which provides the basis of a variational method for determining R 
in terms of the parameter T (Chandrasekhar 1957 b), 


1 1 vs oU 2 
R | dr | dy. [A0|2= —2 | (112) du (F) 
0 a4 = eg heme 


2 le ar [" 0-1) dy {aaswr+t] (5) s ae (5) |} 
- (18) 


For reasons of convenience we will use the following relation from 
which eqn. (18) is derived: 


| ee | (1—p2) du[UA,2U —TU(AV /éz)] 
0 —1 


aS co) 


1 1 
| 4dr | (1 — p2) duU(00/rdp.) 
0 —1 


The application of the Rayleigh—Ritz method consists in substituting 
an appropriate trial function, dependent on several arbitrary parameters, 
in eqn. (19) and minimizing with respect to the parameters. The require- 
ments on the functions substituted for U, V and @ in eqn. (19) are, (1) 
that they satisfy the correct boundary conditions and (2) that they be 
consistent with eqns. (11) and (15). Keeping these requirements in mind, 
we now proceed to develop an appropriate set of functions. 


§ 4. Tue Expansion oF THE PoLorpaL FIELD 


We begin the construction of the trial functions by expanding the 
poloidal scalar, U, in terms of the complete set of functions determined by 


the equation 
jaa) = atU . . . . . . . - (20) 
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where A, is again the five-dimensional Laplacian for axisymmetric 

functions. We find that by using this particular set of functions we are 

able to solve eqns. (11) and (15) with little effort and that, in addition, 

we are able to evaluate the integrals encountered in eqn. (19) analytically. 
The general solution of eqn. (20) which is bounded at the origin is 


T= [AF pranlor) + ByEnsai(or1Cq2%(u)fr%, . . (21) 
n=0 
where C,,3? are the Gegenbauer polynomials and J,,,3. and I,,,32 are 
the Bessel functions of half-odd-integral index for real and imaginary 
argument (Sommerfeld 1949, Erdelyi et al. 1953). 
The application of the boundary condition U(1, 1) =0 implies 


Tie S a | py -F — Lnssie peel | 0,8 32(p)[r32, =. (22) 


n=0 J nisiel) Lal ol 


On applying the second boundary condition at r=1, we obtain for the 
ease of a rigid boundary, 


Fnvsi2(%) T,+3/2(%) 
and for a free boundary, 


J" + 3/0() wx Lene 2(&) Fs : : i : . 4 (23) 


J’ +3/o() LD’ +3ie() 

eee ee meee ee C284) 

J n+3/2() L,43/2(%) ( 
Thus we may write 


ao: Ge dh Qn fT I hve 
U= > 3 Any| pera? — cpus caer, (25) 
n=09=1 J nisiol Gar) n+sio(& nj) 


where «,,; is the jth root of eqn. (23) when the boundary at r=1 is rigid 
and of eqn. (24) when the boundary is free. 


§ 5. THe TEMPERATURE PERTURBATION AND THE TOROIDAL 
FIELD 


We will find it convenient to introduce the following notation before 
solving eqns. (11) and (15): 


‘ Joa (aaty ee en 
fan=iile 7) = ee nee eG 
: i Jn 3/a(o aay) n+3/2(& Xnj) ( ) 


Ve 
J nasil 2(X%p;7') Lies ea) 


I(%nj") = Inin( Ong?) = Te ae y? (27) 
n+3/2\% n+3/2\ Ang 
Ti 3/0(Ong?) (251) 
Sf im(a = 1+:3/2\nj cS Lisia(o At43/2\Ong") ; ag, 1S ret 28 
a 5 EY (0 Xj) Dale) ( ) 
on 
Iilm\ nj?) = 7 -3/2(% pj") Lis sjo(x HH) (29) 


m+3/2 jel hmi) Le 3/2(Cmj) 
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The following relations are easily verified: 


Asl Frain (% ng”) Cyy?!2() [73 >] bart, Oni min (nj!) Cee is, (30) 
A5[9min(%ng?)Om?!2( 1) [79/2] on Det Fenda nit) Cay UE) tats (31) 
oe? 
EME [(1 = w*)C,,9?(u)]= — (n+ 1)(n+ 2)C,32(u), - « + (82) 
é - Cne!*(u) Xnj m+ 2 
Oz Exe ae | = z= aH (5) Sinatin Ong Oma 
m+ 1 
re Gn 5) Imatln(& Ons) | ;m#0, . . (33) 


g =a ) m+2 
Oz | Sn ro “| - S rene) Duin Cnt Cer iE) 


m+1 
“(és =) Tan Cmte) | MA 0, De TG (34) 


o 
5 | Fonlonsr) | = - x4 ) dunlonrICPu), + (85) 
C3 a5 
a2 =| dont Mngt) _ "] SEs eran: To) (630) 
[mC 2%(u) = (n+ 2)10,, Limetewe acu SoMeHs 
Soe u)l=0. me PPO Sn hot aged eee Poe 38) 


Equations (30 and (31) follow from the differential eqn. (20) and eqns. 
(32) through (38) are derived with the help of the recursion relations 
for the Bessel functions and the Gegenbauer polynomials and the further 
relation 


ie Yop ke 39 
dz ar r (OM 2) 
Armed with eqns. (30) through (38), we are now prepared to attack 
eqns. (11) and (15). 
We begin by expanding 06/rdy in terms of the Gegenbauer polynomials, 


7 Ob = by Tt) Cg ie): Sree pa en (40) 


Using eqns. (31) and (32) we obtain a particular solution of eqn. (15) to 
which must be added the bounded solution of 


1 00 
ead ee er ees, ae (LE 
a(75) Ue ri (41) 


The result is the following expression for 7’,,(r): 


S (n+ 1)(n +2) Ans 
LEA) > sam oa Sail a 
ree 


nj 


(Wont) + Br, pepe Soe (42) 
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The application of the boundary condition 00/04=0 at r=1 implies 
Byj= —9(%nj) = — 2. . . . . . ° (43) 
Equation (40) becomes finally 
MiGs meer (n+ 1)(n+ 2) 
= = A ;——_ 
if On 2, 2 a Day 
The determination of V is only slightly more complicated than the 
determination of 6. We begin by writing, with the help of eqns. (33) 
and an the eae for 0U/dz. 


= n+3 
= aps =o ~ oa n+l, j¢n-+1, J (5) Snin+i(n+1, i") 


i. Ao j~¢n-1, 5 (aa :) Inin—1(% n—1, |. (45) 


Using this expression in eqn. (11), we obtain for V: 


co © Ce 3/2 A, ; n+3 
V= ™ » p32 tH) opens (5) Inin+a(%n+1, j”) 


[G(onjr) — aman Cn UH) (44) 


n=0j=1 Ont, j Qn+5 
A =] ° nN 
* Pees & she :) Pape + B,r9? (46) 


where the constant B,,, is to be determined by the boundary condition on. 
V atr=1. Fora rigid boundary, one obtains 


Aveeno Anni 
Satay) ($5) Inin+a(%n-+4, 3) iv ea Ga 


n+, j 2n +5 X19 


) fain-a( On—1,;) + By =, 
(47) 


2n+1 


and for a free boundary, 


Aras ;{(nt+3 
oJ 
Ont, J ee Sie =) Lona, F ninta(%n+1, i) cal 39ninsa(%n+4, i] 
a 1,9 n 
— 1,7) 
On—1, 5 (53) on — ae liat nin— 1(o on 1,3) a ines 1( Xn -13)] 


+nB,;=0, (48) 


where we have made the following additional definitions: 


PF nlm(mj) = J n+3/2(%ms) = DP’ n+-3/2(%m3) oe (49) 


J m-3/2(%mj) Linesia(%ng) - 


Olin Song) = seit’ —- L +3/2(2mj) 


(50) 


m-+3/2(%mj) Lin+sie(%mj) , 


Equations (47) and (48), together with the recursion relations for Bessel 
functions, lead to the following expression for a rigid boundary ; 


~ x n+3 
n+1,79 
ein 2 [ow = i) [Graft (na, 1) = Gninei\naagr ed 


On, j 


>. meth n G: 3/2 
Seat (| alacalena f) “fant, emae]} (61) 
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and for a free boundary, 


be >, E ($)gox(4") —Boy®* | Co°"(u) 
f) 1j 


et Pe 
Sb sie n+3 
mtd, 
J = a = Cae ;) Ininsalo On+1, the Ininpal% n+1, gia a *| 
St een er 
2 Oe ees z) | frail, "+ (=) Frtn-aln-s, sv || 
C,37(u) 
Shale RS a eee seer ae git) 


In particular in deriving eqn. (52) we have made use of the following two 
relations 


On+1,59 ‘nlns(%n44, i) 39 nint1(%n+1, i) = 2 nin+t(%n44, i) ts (Op 41, j)> (53) 


On—1,5F nin—a(%n—1, 5) — 3Fnin—a(%n—a, 3) = — (+3) frin—1(%n—a, 3) +F(%n—2,5)- (54) 
Operating on eqns. (51) and (52) with 0/dz we obtain 


= oO n+4 n+3 
= 2 > { Ansa) (==)| ($55) Inin+2(%n-2, j") 


(n+ 3) : 
=o Cniinge(Cnea git iG 
On+2, 7 = 


el 2) 88) + Ca) 


pies n+3 
ere iG 5) Vain on Ageia 


C,,3?2 
Areas eS at (seri aol Inn-altn-a 0) wee een 5) 


p=—(n+1) for a rigid boundary | 


where 


fF (56) 
p=(n+4) for a free boundary J 


§ 6. APPLICATION OF THE RAYLEIGH—RiTz MetHop 
In order to make the minimization of the functional (19) with respect 
to the expansion coefficients A,,; more tractable, we introduce the 
following shorter notation: 


OSMAN EAC, nese 2 87) 
n=09=1 
ee NEST ne je ee = (88) 
POG pe 
co. 800 C,,3?? ) 
= = » MS [Aye pPnj(t) + AngQnglt) + Ante, ngl)] a : (59) 
n=0j=1 
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Upon differentiating eqn. (19) with respect to A,,, one obtains the 
following equation into which must be substituted the above expressions : 


: | Bays Ae “5, |} =° 60 
OA mr ‘J, ee tb ¢ ee qu | UAeO Le Oz RU Ou ( ) 


Using the orthogonality relation for the Gegenbauer polynomials 


2(n+1)(n+ 2) 5 


he (ie p”) DpC,,,3?(p)C 7,9 (12) = 
—1 


we obtain 


e { > > > A nA mint (nj |F?|nl) — RA, ;A,,<nj|FT |nl) 


OA nr =0 j=1 l=1 


— TA, {Ano (mg |FP|n— 2,1) + Ay (ng|FQ|nl) 
+ Ania (njlFR|n+2,1]} =O. oe eee ome Gd 


where 


1 
Cnjltfaly = PEE rdrP ysl \Falt), (63) 
2(n+1)(n+2) ft 
Torrie di 


ae) I, r drF ,;(7)P (7), (65) 


(nj | FT |nl rdrF,(r)T (7). (64) 


<nj|FP\n— 2,1) = 


dnjlFQ|nty = EE rarP (Qtr), (66) 


; 2(n+1)(n+2) (2 
AV == =e ‘ ; 
(nj|FR|n + 2,1) ame { 1dr ysl0) Blt) (67) 


Performing the differentiation indicated in eqn. (62) we obtain 


Ss Ang [ont (ny |? |mk ) + «,,;4mk| F2|mj Y] 
j=1 


— Ri (inj | FT mk) + (mk|FT mj] 

— TL Km |FQ|mk) + (mk|FQ|mj y]} 

— Anno, jT{(mk| FP |m — 2,5) + (m—2,j|F Rlmk y} 

— Anis, T{(mk|PR|m + 2,5) + (m+ 2,j|FP|mk)y} =0. (68) 


From the self-adjoint character of eqn. (10) or from the expressions to 


be derived in §7 of this paper one can demonstrate the following 
symmetries of the matrix elements 
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Ame (mj |F2|mk) =onj4{mk|F2\mj), . . 2. (69) 
(mj|FT|mk) =(mk|FT|lmj), . . . . . . (70) 
(mj|FQ|mk> =(mkiFOlm), . « . .. »« (71) 
(mk|FP|m—2,j5)= (m—2,9|FRimk). a et ea Ce) 


To determine the solution of the infinite set of eqns. (68), we set the 
determinant of the coefficients of A,,; equal to zero. The resulting 
polynomial of infinite degree in R is 


[Sranl nx onj|E2}mk) — Rj FT mk) —T (ng FQ|mk 
= 8 n—2, ne <mk|FP|nj > i Ones L (mk|FR\nj >| =0. (73) 
The critical Rayleigh number is now determined by approximating the 
above equation with a finite number of terms and computing the lowest 
root numerically. As may be easily demonstrated, we need only compute 


matrix elements with even values of m since functions with the subscript m 
odd do not lead to conservation of angular momentum. 


§ 7. EVALUATION OF THE Matrix ELEMENTS 


The first integral we wish to evaluate is 


n+1)(n+2) 


dnjfee|uty = 2+ Wie [ raiflonsr)floni) Sey 


Using the relations 


C,8%(u)]  C,2%(u)d? 1d (n+3)? 
As | Plan pal2 | == aR 773 a ae Zs = f(opj7) (75) 
and 
d? 1 d nN -+ 3. 2 
dr? os rar xe = |flew= = ng J (Onl), - (76) 


we can readily show 


2n+3 ne fore 
2(n + 1)(n +2) (Onn nj ) (nk|F [nj y 
~ On i LOnjh "(ong )G (On) a On td (%nj)9 (nx) | 


i Ong Omid (nn) I (%ng) va Ont (nn) (nj) I (77) 
Now we use the boundary conditiuns at r=1, 
F(cng)=9, J (Ons) =9: for a rigid boundary ; 


flonj)=0, f'(Ong)= —%ng for a free boundary ; 


to show 5 a 
(nk| F?|nj )=0, nj FA OXnn- ee a 1S) 
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In case opj=o,, we use l’Hospital’s rule to show that for a rigid 
boundary 
(n+1)(n +2) 
se i ee 
(nk|F*\nk) = (79) 
and for a free boundary 
(n+ 1)(n+2) 


2n+3 One I (nk) . Ms . (80) 


(nk|F?|nk) = — 


Using the same techniques, we obtain for either set of boundary 
conditions 


1 DE 2 P 
[ rdrf ona) 9 onj0) = FO [ong ng) Ane (nd (81) 
(0) Onk — nj 


3)2 
[rae oon) = HE On? + (01 AEP, (82) 
0 n 
i‘ ir Arf (On n7)I nim(%ms”) aa ee {206 42 Omg alm (mj) 
= On I (aK) Inlm(%mj) = in eng (oa) T rater Gn ales . - (83) 


With the help of the recursion relations for Bessel function we obtain 


[ rdrf(a,.r)re ts? = a pias — Lenk) | (84) 
~ 0 Xnk Xnk 
Equations (78) through (84) and further use of recursion relations 
enables us to write the rest of the matrix elements as 


n+1)?2(n+ 2)2 on +3 
(nk|PT nj) = et Dn 2)" { Seles 
2n+3 nk nj 
>) ! : , 
nk — nz Xnj Xnk 


dn PT ay = SOE LAF Cong) +9" Caad 


2 7 2Qg 
a Rc) ‘(a nail (86) 
nk nk 
<nk|PQ|nj ) = pe ee _ F' (nk)\(2H+3 9" (o%ns) 
(2n+ 3) nk Xn ke ngs On; 
2n? + 6n+1 Mae 
Qn+lan+5) lly), (87) 
(nk|FP\n—2,j) = 2@=Vnln+(n+2) 
(2n—1(2n +3) cing!—oyo 8 


x {inns nell ns i) 


2 ! One , 
a oni Gan + Leni) (20-1 _ 9" (on-2,) (88) 
Xnk Cli On—9, fa a : ? 
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{nk|FR|n+ 2,j)= 2(n+ L)(m+ 2)(m+3)(n+4) 


(2n + 3)(2n+7) One Sate, 4 


x { ~ Wy Kn +05 id (nn) S (nee 3) 


9 , , 
toga t( BHT + 9 Casas) (2048 : Tens) (39) 


2 
Xn+2, j An+2, 7 On ke On ke 
The evaluation of the coefficients in eqn. (73) for approximations up to 


nine by nine is given in the Appendix. In these tables a factor of two has 
been cancelled. 


§ 8. NUMERICAL RESULTS 


Because of the rather tedious calculations involved, it was decided to 
approximate eqn. (73) with determinants of order less than or equal to 
nine. As may be seen by examining the apparent convergence of the 
values of Rin table 1, eight parameters give a sufficient approximation for 
values of T less than 10° in the rigid boundary case and less than 105 for 
the free boundary. 

The values of n and k assumed in our approximations are: n=0, 2 and 
4 and k=1, 2 and 3. 


Table 1. The Value of R for the Onset of the Lowest Mode of Instability 
for Various Values of the Taylor Number 


For the case of a rigid boundary 


lst approx. | 2nd approx. | 4th approx. 6th approx. 9th approx. 


8-0611 x 10° 421 x 103 


8-0 =s 
9-4882 x 103 | 9-3855 x 108 | 9-3383 x 10° 
oul 


-7306 x 104 | 1-6991 x 10* — , 
== 0:56993 x 10° | 0-53344 x 10° 
— —s 0-32122 x 10® | 0-29105 x 108 


For the case of a free boundary 
i enne eo eee 


T | Ist approx. | 2nd approx. | 4th approx. 6th approx. 9th approx. 


0 | 3-0916 x 10% | 3-0912 x 10° pie ate - 
3) 9-4715 x 103 | 9-2493 x 10 -9870 x ~- - 

we ve oi as 2-6395 x 104 | 2-5517 x 104 = 

105 = = 0-96961 x 10° | 0-91229 x 105 


The results of the above tables are shown graphically in fig. 1. 


4¥ 


P.M. 
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Fig. 1 


logR 


9-5 
oh 


5.0 


4.0 
G8) 


3.0 
0 2 3 4 5 6 logT 


The dependence of the critical Rayleigh number for the onset of convective 
instability on the Taylor number in a rotating fluid sphere. The curves 
labelled (1) and (2) are for the cases of a rigid and a free bounding 
surface. 

Table 2. The Characteristic Vectors Describing the Lowest Mode of 

Instability for Various Values of the Taylor number 


For the case of a rigid boundary 


T=0 | =103 


0-99982 | 0-99859 


2 | 0-01917 | 0-02847 


0:04470 


0-00362 


=10° 


1: 


0-82340 
—0-33881 
0-42816 
0-03278 
0-14274 
0-04958 


0-85096 
0-37755 
0-02874 
0-33137 
— 0-00423 
0-00159 
0-13141 
0:05937 
0-03938 
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Table 2 (cont.) 
For the case of a free boundary 


T=0 = 108 =10! =10° 


Ay, | 1-00000} 0-99662} 0-80118| 0-61020 
Aog | —0-00151 | —0-04190 | —0-36978 | 0-63219 


An =m —0-00388 | —0-01094 | —0-11656 
Pr = 0-07050] 0-45896 | —0-39936 
As = a —0-03995 | 0-16196 
ao == = = 0-01481 
de = a 0-09497 | —0-16842 
ot = = an 0-00975 
ae = a = 0-00385 


Though these are too few values for any quantitative conclusions, we 


have graphed the results in figs. 2 and 3 to show qualitatively the 
behaviour of a few of the components. 


Fig. 2 


; isti tor describing the 
The dependence of the components of the characteristic vec bing th 
ire mode of instability on the Taylor number in a rotating fluid 
sphere with a rigid boundary. 


4X2 
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Fig. 3 


The dependence of the components of the characteristic vector describing the 
lowest mode of instability on the Taylor number in a rotating fluid 
sphere with a free boundary. 


§ 9. ConcLUDING REMARKS 


Because of the rather slow convergence of the successive approximations 
to the Rayleigh number for higher Taylor numbers, it is felt that further 
computations using the method described herein would not be worth- 
while for this problem. Since the only self-gravitating spheres we can 
observe are of astronomical dimensions and (because of the dependence 
of T on the fourth power of the radius) have rather large Taylor numbers 
(e.g. ~ 10° for the earth’s core) ; it is evident that an asymptotic relation 
giving the dependence of R on T when T is very large is needed. 

A second point which may be mentioned is that there is certainly no 
good reason for assuming that instability does not appear as overstability 
rather than convection. Indeed, it has been shown that in the case of the 
plane Benard problem with rotation, instability does generally set in as 
overstability in the case of liquid metals such as mercury (Chandrasekhar 
and Elbert 1955, Fultz and Nakagawa 1955). Before any more extensive 
calculations are attempted for the problem in a sphere, the overstable 
case should be thoroughly investigated. 
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APPENDIX 


Further Tables of Numerical Results 
Table 3(a). Roots of the Equation f(«,;)=0 


1 | 5-267568 | 7-748590 | 10-10830 
2 | 8-506951 | 11-19090 | 13-73380 
3 | 11-68768 | 14-47661 | 17-12845 


Table 3(b). Roots of the Equation f(a,j) = — %; 


4-180881 | 6-813002 | 9-232055 


2| 7-573795 | 10-309365 | 12-881869 
3 | 10-802354 | 13-618716 | 16-289137 
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Table 4. Numerical Values of (nk|F?|nk) 

n,k| Rigid boundary | Free boundary 

Ol 1-333333 6-330580 

02 1-333333 28-60704 

03 1-333333 63-92143 

2] 3°428571 54-65893 

22 3°428571 146-4251 

23 3:428571 271-3658 

41 5454545 175-8249 

42 5454545 378-4838 
5454545 632-3310 

Table 5. The Numerical Value of the Constant Appearing on the 


Principal Diagonal in the Basic Determinants 


n,k| Rigid boundary 


Free boundary 


The constant term appearing off the principal diagonal is identically zero. 


5:13309 x 10? 
3-49144 x 108 
1-24401 x 104 
Ol o7o sa lOs 
2-68870 x 104 
752923 x 104 
2-84735 x 104 
9-70268 x 104 
2-34748 x 10° 


9-67133 x 10? 
4-70649 x 104 
4-35201 x 10° 
5°88825 x 104 
8-27014 x 10° 
4-66735 x 10 
6-38621 x 10° 
5:21114x 108 
2-22590 x 107 
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Electrons in Polar Crystals} 


By G. L. Szweiu 
Department of Applied Mathematics, The University, Liverpool 


[Received July 16, 1958] 


ABSTRACT 


The theory of the properties of an electron in a polar crystal is investigated 
by means of a model which takes account of the atomicity of the lattice. The 
method used is essentially that of a tight-binding approximation in which the 
interaction of the electron with the lattice displacements is built into the 
theory from the start. This method is valid only in cases where the tight- 
binding approach is appropriate even when the effects due to the lattice dis- 
placements are ignored. : 

Thus, the eigenstates of the system are constructed as linear combinations 
of localized states. In each of the localized states, the electron is bound to one 
positive ion, and displaces the neighbouring ions. The transfer of the electron 
between ions is therefore accompanied by a transfer of the mean position of the 
neighbouring ions. The effective mass of the electron is consequently 
increased as a result of the displacements it carries with itself. It is found that 
this effective mass increases exponentially with the electron—phonon coupling 
and can thus be extremely large in the strong coupling case. It is also found 
that the effective mass increases with temperature, since the random thermal 
motion of the ions opposes the motion of the lattice displacements accom- 
panying the electron. 


§ 1. INTRODUCTION AND DiIscussION 


AN electron in a static periodic lattice behaves, in many respects, like a 
free particle, whose effective mass, m°, depends on the potential due to 
the lattice. In polar crystals, the interaction of an electron with the 
polarization field, generated by the lattice vibrations, increases its 
effective mass to m*, the polaron effective mass. 

The existing theories of the polaron are based on a model due to 
Frohlich (Frohlich 1954, Allcock 1956), in which the polarization field is 
treated as a continuum, and in which the electron is treated as a free 
particle, of mass m°. According to these theories, the polaron effective 
mass is not much larger than m° unless m? is itself large. In this latter 
case, however, the continuum model breaks down, since it leads to a 
picture in which the electron behaves like a wave-packet that is localized 
within a distance less than the lattice spacing. 

A different approach to the theory of polarons of very high effective 
mass has been discussed by Frohlich (1957), in connection with the Debye 
type of dielectric losses that occur in some ionic crystals, due to electrons 


ie wee a a 
+ Communicated by Professor H. Frohlich, F.R.S. 
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bound to lattice defect centres (e.g. colour centres). Fréhlich considers 
cases in which the original effective mass m® is large. These are precisely 
the cases where the Bloch tight-binding approximation would be 
applicable if the lattice displacements were neglected. The essence of the 
new approach is that, in such cases, one may also calculate the eigenstates 
of the electron-lattice system by means of a tight-binding approximation 
in which the interaction of the electron with the lattice displacements is 
built into the theory from the start. 

Thus, in determining the properties of an electron bound to a defect 
centre, one starts from a state in which the electron is localized on one of 
the positive ions near the centre and displaces the positions of the 
remaining ions. There will be a number of equivalently placed positive 
ions near the centre, and hence there will be a degenerate set of localized 
states of the above type. The energy required to localize the electron on 
any one of these ions is simply the overlap integral governing the transfer 
of the electron, together with its accompanying displacements, between 
neighbouring positive ions. This integral is given by 


W~Wexp(—A2/A,2), . .... . (1) 


where W® is the value the overlap integral would have if the effect of the 
lattice displacements were ignored, and the exponential factor arises 
from the transfer of the lattice displacements. In this latter factor, A is 
the displacement in the mean position of a positive ion, due to the presence 
of an electron on a neighbouring ion, and A, is the amplitude of the 
optical zero point vibrations of the positive ions. It is clear that the 
exponential factor may become extremely small in the case of strong 
coupling. In this case, the energy W, required to localize an electron on 
a positive ion near the defect, may be less than «7. If this is so, the 
electron will behave like a classical charged particle which is confined to 
the positive ions in the neighbourhood of the defect. On application of 
an electric field, the electron will therefore behave like a classical dipole, 
and will consequently lead to Debye dielectric losses. 

In the present paper, we shall employ the same tight-binding approach 
to determine the theory of polarons of very high effective mass, for the 
case when there are no lattice defects. Our theory should therefore be 
applicable, for example, to the experiments of Morin (1954), which have 
shown that the electronic mobility in some ionic crystals (e.g. NiO and 
a-Fe,O3) are extraordinarily low. Such low mobilities can presumably 
only be explained on the basis of an extremely high effective mass. 

Our theory will be based on a model in which an electron interacts with 
a polar lattice whose basic structure is that of two interlocking cubic 
sub-lattices (one for each set of ions). In order to determine the 
eigenstates of the system, we first obtain a localized state in which the 
electron is bound to one positive ion, whose motion it follows adiabatically. 
In this state, the presence of the electron on one ion will lead to displace- 
ments in the positions of the neighbouring ions. Moreover, there will be 
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a degenerate set of such states, one corresponding to each positive ion. 
By taking linear combinations of these localized states, we may remove 
the degeneracy and obtain the eigenstates of the system. The energy 
levels corresponding to these eigenstates will form a band, whose width 
is given by the overlap integral governing the transfer of the electron, 
and its accompanying displacements, between, neighbouring positive ions. 
This band-width will therefore be given by (1.1), provided that we now 
interpret W as the polaron band-width and W® as the corresponding 
electronic band-width when the lattice displacements are ignored. We 
see from this formula that the polaron band may be very narrow, and thus 
the polaron effective mass may be very large, if AS A,. As we shall 
show in §7, this criterion can be satisfied for reasonable values of the 
characteristic constants of the model. 

The polaron band-width is temperature dependent, and is given by the 
thermal average of the overlap integrals governing the transfer of the 
electron and its accompanying displacements. In fact it decreases with 
temperature, which means that the polaron effective mass increases with 
temperature. The essential reason for this is that the random thermal 
motion of the ions opposes the transfer of their mean positions. 

It should be pointed out that the electrons whose low mobilities were 
measured experimentally were apparently d-electrons. However, in 
order to simplify the mathematics, we have confined our theory to a 
polaron band based on electronic s-states. The present theory can 
therefore be regarded only as a qualitative, rather than a quantitative, 
approach to the problem of the properties of electrons in polar crystals 
of the type mentioned above. 

In a second paper, it will be shown that our theory leads to the 
prediction of a very low electronic mobility. 


§ 2. THE MopEL 


We shall employ a model consisting of an electron, of mass me}, Moving 
through a crystal built up of positive ions, of mass M, and negative ions, 
of mass M’. It is convenient to introduce the total mass and the reduced 
mass, for a pair of opposite ions, 1.e. 


iM = Mia, Minh GNF 2 <3 ieee Rois (2.1) 
and 
EY EE 


and also the reduced mass for an electron moving in the field of a positive 
ion, i.e. 
Mie hae ere TS) 


= 
M+m™ei 
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The mean positions of the positive and negative ions will be denoted 
by a,, a,’ respectively, where 


3 
a,=a> S$, U, = 4(81, 8g, 83), a te. Hie ee Sk ee 
i=1 


the numbers sj, 85, 8, being integers, and uy, u,, us forming a mutually 
orthogonal set of unit vectors ; and 


a, =a,+aa, a= (a, Oo, Xs), ° e . ° ° (2.5) 


the three components of « all lying between 0 and 1, and being independent 
of s. 

The position of the electron will be denoted by x, and the displacements, 
from their mean positions, of the sth positive and negative ions will be 
denoted by X,, X,’, respectively. The relative displacement of the sth 
ions, and the mass centre corresponding to the displacements of these 
ions are thus 


Y= Xi XS He ee ee 

and MX,+M'X,’ 
= s 8s DF 
[a (20 


It follows from these equations that we may express X,, X,’ in terms of 
Y;, Z, by the equations 


Ke Ze EO a) =) oe ee ee 
and X{=z;—(l=¢)y¥, eee ee ee 
M’ 
h = —_.. 2 
where g Wee (2.10) 


It will be assumed that the optical and acoustical modes are described 

by the y’s and z’s, respectively, though this is really the case only for the 
long waves (cf. Born and Huang 1954, p. 82). 

The Hamiltonian for the system is 
H=H(el)+H(lat)+ A(int), . . . . . (2.11) 
the three parts referring to the electron, the lattice and the electron— 
lattice interaction. Here, 

: Ce MLO 

H (el) = — —— < Ue eee ee 


2me1 Ox2- 


The interaction H(int) is the potential energy of the electron in the field 
due to the ions. Thus, 


A(int)=>[V(x—a,— X,)+ VAUX 8) eee Oat 


where V( x— X), V'(x— X’) denote respectively the potentials at x due 
to a positive ion at X and to a negative ion at X’. 


P We shall assume that the lattice Hamiltonian may be split up into the 
orm 


H(lat) =A. + He eet 
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the two contributions to H(lat) being due to the optical and acoustical 
vibrations, respectively. These contributions may be written as 


h2 


Donte y ate 

opt= 2 DM, dyz t Vovly), +. + + (215) 
h2 2 

and Hye =D sa 2 + Vacl2) toa 3 TORIC) 
8 a a s 


where V opt and Vac are functions of the y,’s and z,’s respectively. The 
formulae (2.14)-(2.16) rest on the hypothesis that the total potential] 
energy of the lattice can be written in the form Vopt(y)+Vac(z). It is 
easily seen that the lattice kinetic energy, as given by (2.14)-(2.16), is 


which, in view of (2.1), (2.2), (2.6) and (2.7), is equal to 
sf oe 
‘ aM OX ee OM OKT" | 6 
the correct formula for the total kinetic energy of the ions. 
The Hamiltonians Hop; and Ha. may now be simplified by writing the 
ys and z’s in terms of normal coordinates. We shall neglect the 


transverse modes, since these interact less strongly with the electron than 
do the longitudinal modes. Thus (cf. Wilson 1953, p. 251), 


> 1/2 
ye > 12 00s W ae sin Wa) Woe. me (cal) 
$ NM, a 
and 
2 1/2 ; = n 
ja (war) >» [Eé,3 cosw.a,t+é,,sinw.a,Jw. . (2.18) 


In these equations, the é’s are normal coordinates, N is the total number 
of pairs of ions, and the primes over the >’s denote summation over half 
of the occupied region in w-space. This region is given by 


T | 
= = Wj, Ws, W3< — 
ae et ee (2119) 
| 
J 


where (W 1, We, W3)=W.- 


The Hamiltonians H,,, and H,, may now be expressed in the forms 


A c= > Dal wg? + Pw Swi" J | 


j—t w 2 (2.20) 
and Lab 2 Pye 4 
NBL > ps Z| ows T Vwi Swi 1; J 
j=3w 


where nw; is the momentum conjugate to €,,, and v,; is the frequency 
for the mode w, j7. Hence, we may rewrite (2.14) as 


Tees eS Hee ee oT) 


g=1 ow 
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where 
Da bene, = $n; + Ver ewr le . . . . ‘ (2.22) 
and 
Nwj —tho/0Ew;- 4 A 5 A Fs (2.23) 
We shall assume that 
= j=1, 2 tical modes 
Vwi = V5 for ) , 2 (opti ) (2.24) 
and Vwi = cw, for j= 3, 4 (acoustical modes), 


where v (the infra-red absorption frequency) and c (the speed of sound) are 
constants. These equations for the frequencies are certainly valid for 
the long waves. 

The eigenstates and energy levels of H(lat) are given by 


4 
0,= I] LTO: Lon) gers ee (2.25) 
j=. ow 
and 
a (lait) = > D (hg +4) 9 mpm e pe ds) 
j=) Ww 


where O(n,,;, &w;) is the eigenstate of H,,; in which there are n,; quanta. 


§ 3. LocaxizEp STATES 
We shall now determine the wave functions for states in which the 
electron is bound to one positive ion, and displaces the positions of the 
remaining ions. (We shall ignore the polarization induced in the ions 
themselves.) For this purpose we split H into two parts, 


Heo, Wada lath ee ee 
where 
= V(x—a,) a X=a,/)"sa oh oh ee ose 


is the potential energy of He electron due to all the ions except the sth 
positive one, when they are in their 1 mean positions. Thus, by (2.11), 
(2.13), (3.1) and (3.2), 


H,=H(el)+ V(x—a,— X,)+ H(lat) + > IM X= ap ye al 
7a Koa Xiah eR): 


which means that H, is the Hamiltonian for a system in which an electron 
interacts with the sth positive ion and with the displacements of the 
remaining ions. 

The localized states we shall consider will be stationary states of ish 
given by the form 


p,=o(x—a,—X,)y,(lat), . . . . .. (3.4) 


where y, depends on the lattice coordinates € only. w, is thus a state in 
which the electron follows the sth positive ion adiabatically. 
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It follows from (2+1)-(2.3), (2.6), (2.7), (2.11)-(2.16), (3.3) and (3.4) that 


h* 2 
Ayb.=x, j ak 2m ox2 Ae V(x—a,— x) | o(x—a,— X,)— oe . oe 
+4] H(t) + & (Px a,— X)— Voxe—ay)}+ SIV xa — X/) 
FS t 
mV (ee a VTi | pee en a nea 


We shall assume ¢ to be an eigenstate of the electron, in, the potential due 
to the sth positive ion, thereby assuming that the lattice spacing is large 
enough for this atomic state not to be appreciably affected by the displace- 
ments of the other ions. Thus, we shall assume that 
2 2 

E — + V(x—a,— x.) | $(x—a,— X,)=e9h(x—a,— X,), (3.6) 
€) being the energy corresponding to the atomic state 4. Further, we 
shall stipulate that ¢ is an s-state, since the non-degeneracy of such atomic 
states leads to mathematical simplifications. ,, on the other hand, will 
be obtained from the condition that the energy of the system H, for the 
state 4, is stationary. ¢ and y, will be chosen so as to be normalized 
with respect to the electron and lattice coordinates respectively. Hence, 
it follows from (3.5) and (3.6) that, since 


$+Vd x= |ev¢ d’x=0 (¢ real), 


then the energy of the system H,, for the state y,, is 
(h,*| H, |b.) = €o+ (x.*| H,(lat) | x,); BP toe ee 


where 


H (lat) =H (lat) + {loo a,— X,)|2} > [Vx—a,—- %) — V(x—a)] 
t#s 
+ >[V'(x—a/ — X,')— V'(«- ai)! | Oo Xa eae) 


It follows from (3.7) that, if 4, is a stationary state of H,, then y, is an 
eigenstate of H,(lat), with energy, say H,’ so that 

H,(lat)y,=L,'x, hie yea Aan wee RE), 
and (b,+|H,|%,)=L, + &o- Pea ee (3210) 


We may simplify the formula (3.8) for H,(lat) by neglecting all powers 
of the X’s and X’’s higher than the first in the integral appearing in that 
formula. Hence, 


H (lat)=H(lat)— > X,. [ Ib(x—a,) PV V(x —a,) d?x 
t#s - 
—>X/ . | |p(x—a) PV (x— a’) Bx. 
; J 
This formula can be rewritten, in view of the lattice periodicity as 
H (lat) = H (lat) — x] X,. | [POOP VV — a+ a,) dex + XY | d(x) 2 
t ¢ d 


x VV'(x— a, + a,) x | + X,. |b(x) |? V V(x) d? x. 
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The last term in this equation is zero, since |¢(x) |? is an even function of 
x, while VV(x) is an odd function of x, since V(x) is a radial function. 


Hence, putting 
B(X)= [Ibo PV [gor — X)—(1—g) V(x — X + aa)] dx 
and 
C(X)= [BOOP VIP X)+V'(x— X +aa)] dx, 
it follows from (2.5), (2.8) and (2.9) that 
H (lat) = H (lat) — orale —a,)+z,. C(a,—a,)]. 


The last equation may be ae using (2.17), (2.18), (2.21), 
and (2.24), as 
4 
H (lat) a > Daler + V3 a So oa Ewis)1— 1 


j=lw 
where 


4 
1 2 
ap > aire wijs ? 
where the é’s are all real, ae nie 


= = @ 1/2 
6 prea) lal me . ; 
oe fe wes KH (war) a s B(a, a.) exp (iw * a,) 


and 
Oks as Vee = (cw)? (Gr : Pah C(a a.) exp (iw . a,). 

1.e. 

ewig t tEwas= N—U2B,, exp (iw a.) | 
and 

‘shes + 1b wa a tO exp (tw a Me | 
where 4 

9) /2 
b= (zr = “J >w . B(a,) exp (iw . a,) 


and 
: 5 fon eee 
snr 2)" dw. C(a,) exp (iw. a)). 


aC t 


By (2.21) and (2.23), we may rewrite (3.12) as 
A (lat) = TH (lat)T + — e,, 


. of = 
T =exp E 5 a Dass | 
fil 


(3.11) 


(2.22) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


is a unitary displacement operator. Consequently, the eigenstates and 
energy levels of H,(lat) are related to those of H(lat) by the equations 


Xs— Xsn = LEAS), 


(3.18) 
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and 
E, = E,(lat) — e,. 
Hence, by (3.4), (3.9) and (3.10), the localized wave functions are given, by 


b= Yon =P(X—a,— XX ons Pe enh) 
and the corresponding energy levels are . 
(bent |, |b.n) =H, (lat) + Ey — €. oe os ee (Bec) 


§4. Tor Banp SrructTURE 


For a given distribution of phonons (i.e. given @,,), there will be N - 
states y,,, one for each positive ion. These states form a degenerate set. 
By taking linear combinations of them we may obtain stationary states 
of the system. Thus, we shall seek stationary states of H, that of the 
form 


"Pie: > cabs be ae eee ek 
By (3.1) and (4.1), the energy E of the system H, for the state W, is given by 
dee (Pbon| Le vel. on U,| Bin) =0. 
st 
Choosing the c’s so that # is stationary, we obtain N secular equations 
Delben* | 1 Hi, me EF bin) =0. eer it Fen) 
Fi 
In order to solve this set of equations, we adopt a tight-binding approxima- 
tion (cf. Bloch 1928), and put 
(ben Lin) ra Sy: ss A . s (4.3) 
Thus, if we assume that the ¥,,,’s are eigenstates of H,, it follows from 
(3.22) that 
(bent | Hol tin) = (En (lat) + ep—ey)dy - - (44) 
The matrix elements (:/,,.7|U.|%,) depend on a,, a, only through their 
difference, in view of the lattice periodicity, so that we may put 
Chey | U,| bin) = W (a ca a,). x = = 3 . (4.5) 
By (4.3)-(4.5), we may rewrite (4.2) as 
Yel# + €y— €9 — H,(lat)) 3x a W (a aes a.) ] = 0. 
t 


The solution of this set of equations for # and the c’s is 
C,= exp (ik. a,), ape GAS ECD) 


where 
k=(Ky, has ha): ote dee eee (400) 


is an electron wave vector which lies in the region 
ffi TT P 
Se CAD al od eS ee ee (eS) 
a a 


cessary to consider other terms in this sum, since the atomic 


It is unne 
y egenerate, owing to the fact that ¢ has been chosen to be 


energy level ¢ is non-d 
an s-state. 


Z 
P.M. 4 
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and the energy H is given by 
B=«,—«,+L,(lat)+ >W,,(a) exp (vk . a). 
t 


It follows from this last equation, together with (4.1) and (4.6), that 
the wave functions and energy levels under consideration are 


Pah = > Oxp (tk. aslben a Ss 
and 
E=E£,,,=£,(lat) + €cn; , ee eae 
respectively, where 
Eun = €9 — 1 + > Wn (a) exp (tk . a,). eer ee 
t 


Since Z,,(lat) is the lattice energy, <,,, may be regarded as the polaron 
energy for the wave vector k, when the lattice state is 0,. 
The energies W,, may be expressed, using eqns. (3.2), (3.19) and (4.5), 
in the form 
W,,(a,— a.) -~ OG oa | W a | Xin)» oie px (4.12) 
where 


Wa= I $+(x—a,— X,)U,d(x—a,— X,) d?x. 


Assuming that the displacements X,, X, are much less in magnitude than 
the lattice spacing, we may neglect these terms in the above formula for 
W.,,. In this case, 


Wy= [Pea Ud —a) Px=W%a—a), . . (4.18) 


st’ 


a quantity depending on a,, a, only through their difference. Consequently 
by (3.18) and (4.11)-(4.13), 


Exn = (0,7 | Eop/( k) | @)); . : . ° ° (4.14). 
where the operator eop(k) is given by 
€op( k) = e)—€, + )W%a,)T tT, exp (ik . a). . . (4.15) 
t 


§ 5. TEMPERATURE DEPENDENCE OF POLARON BanpD 
It follows from (4.14) that when the lattice state is @,, the polaron 
energy is the expectation value of the operator e€op(k) for that lattice 
state. Consequently, the polaron energy for the wave vector k at 
temperature 7' is the thermal average, taken over all lattice states, of 
€op(k). It is therefore equal to 


e(k) = SPléon( k) exp (~ BH (lat) 


Splexp(—BH(lat))} — ” ae 
where 
Beil ee Pees SR. 


and « is Boltzmann’s constant. It will be shown in Appendix 3 that the 
mean square deviation from ¢(k) of the energies ¢,, is negligible for a 
macroscopic system. ¢«(k) may therefore be justifiably regarded as the 
polaron energy for the wave vector k. 
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It follows from (4.15) and (5.1) that 


e(k) = €)— €, + }W(a,) exp (ik . a), ee, Ooh 
where 


W(a,) = W°(a,)F(a,) care ued (Des) 
and 


F(a) = Sp[To*T,exp (—BH(lat))] 
Splexp (— BH (lat))] 
It will be shown in Appendix 2 that 
Wt(a)=W(a)=W(-a) ... . . (5.6) 
In view of this, eqn. (5.3) may be rewritten as 
(kK) = eg —€, + W(0)+25’W(a,) cos k . a, 
t 


(5.5) 


the prime over the } denoting summation over half of the lattice points 
a,~0. We shall now ignore all the W’s except those concerned with the 
transfer of the electron between nearest neighbouring positive ions. We 
shall therefore retain in our last equation for ¢(k) only the terms for which 


a,= Quy, aus or aus, 
the u’s being defined in (2.4). Hence, using (4.7), 


3 
e(k) = ey9—€, + W(0)+25'W, cos ka, Mas! Gee 
(=1 
where 
W,=W/PF;, Spats Seles Ah EE 
F,=F(au) 
and (5.9) 


W/?= W°(au). | 
It is easily seen from (5.7) that 
3 
e(k) = €(0) — SW ah? + 0(k4). 
l=1 


Consequently, the effective mass of the polaron, for propagation in the 
direction of u,, is 
h2 


oe 5.10) 
2W a? ae 


My = — 


§ 6. CALCULATION OF ™, 


It is seen from eqns. (4.13), (5.5), (5.8) and (5.9) that W, is the product 
of two factors. The first of these, W,°, is the integral corresponding to 
the transfer of an electron between nearest neighbouring positive ions in a 
rigid lattice of the same basic structure as the one under consideration. 
It is related to the effective mass m,°, of an electron moving in the direction 
of u,, when the lattice is rigid, by the equation (cf. (5.10)), 


hi eee res (61 


(0) = 
Mm, = . 
l 9 W a2 


4Z2 
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The second factor, F,, in the formula for W;, arises as a result of the 
transfer from ion to ion of the lattice displacements that accompany the 
electron. 
We may now express m, in terms of F, and m,°. For, by eqns. (5.10) 
and (6.1), 
m =m9/F,. tba OY se a Cae 


The effects of the lattice displacements on the polaron effective mass are 
therefore given by Ff. We shall now calculate this quantity. 
By (2.21), (2.22), (3.17) and (5.5), 
4 
F(a)= Ly LU fj oe ee ee 
j= w 
where fas Splexp {i(Ewso — Ew) 2w;/} exp (— BHy;)] (6.4) 
* Splexp (—PH,,;) | 
It will be shown in Appendix 1 the solution of (6.4) is 


2 


oe aaa 
F wit =eXP E a coth (1Btings | ; 2 2 MGA 
It now follows from (6.3) and (6.5) that 
Ps or ane eee 
B(a,)= exp] = 4h 2 > erie wines cw). coth ($Btin«) | ; 
j=1 w 
and therefore, by (2.24) and (3.14), 


ee : 
F(a, =exp | - > d'sin? $w.afv| By coth (4%) 


Nh j=1 w 


+cw|C\,|? coth (4Btiew)} |. 
Thus, by (5.2) and (5.9), 


F,=exp (—y,) ay eae vip Uo ee are OLE) 
where 
y= y (opt) + y,(ac), MP ore, we of ket 
and where os ve) a 
v h 
y(opt) = Wi; coth (x7) >| BaP sin? saw, | 
and 


(6.8) 
yi(ac) = 


: : hew 

2 sin? law 
w|C\,|? sin? Saw, coth a 
y(opt) and y(ac) are thus dimensionless quantities which increase with 
temperature, and which depend on the coupling of the electron to the 
optical and acoustical vibrations respectively. 
. The polaron effective mass, m,, and the energy W, may now be expressed 
in terms of y, For, by eqns. (5.8), (6.2) and (6.6), 


2, XD (4) (6.9) 


(6) 
Ni 


and 


Wi = WwW, exp (— y;)- . . . . . (6.10) 
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An order of magnitude approximation to this last equation, for the case 
7’ =0, in which the acoustical vibrations are ignored, is given by (1.1). 

Equation (6.9) shows that the polaron effective mass increases exponen - 
tially with y,, and may thus be extremely large if 


y>1. aC Se acne! ee (OSE Ly 
We shall show in the next section that this condition may be satisfied 
for reasonable values of the characteristic constants of the crystal. 


Moreover, by (6.8) and (6.9), m, increases with 7’; and it may be seen 
from those equations that 


m =m) exp (T/T), for eT shy, hcla, . . . (6.12) 
where 


2m 
(xT) t= Wie >’ {| Bw? + | C2} sin? daw,. ae Teme) 


The increase of m, with 7 might have been anticipated from eqns. (5.9) 
and (6.2)-(6.4), according to which m, is inversely proportional to the 
quantities f,,;,, which govern the transfer of the displacements in the 
various lattice modes. For /,,, is the thermal average of exp (iAn,;), 
where A is a constant. With rising temperature, the momentum distri- 
bution for the oscillator H,,; will tend increasingly to include high as well 
as low momentum components, and the thermal average of exp (tAy,,;) 
may therefore be expected to decrease with temperature. In fact, at 
T=co, when the momentum distribution will be completely random, 
this thermal average will clearly be zero. 

To sum up, we may say that the random thermal motion of the ions 
opposes the transfer of their mean positions (i.e. f,,;, decreases with 7), 
and, as aresult of this, the polaron effective mass increases with temperature. 


§ 7. MAGNITUDE OF y, 
By eqns. (6.7) and (6.8), 


V1 = ves ° . d . . ° : . (Tab) 
where oe Sis >’ | B, [2 sin? daw, ee ee ee ea) 
Ni 


is the value of y,opt) at 7'=0. The criterion (6.11) will therefore be 
satisfied. if 
yo >I. 
In order to ascertain whether this latter inequality is valid, we require 
only a calculation which yields a correct order of magnitude for 1°. For 
this purpose it suffices, using eqn. (7.2), to obtain formula for B, in terms 
_ of w for the long waves and then to assume that formula to be approximately 
correct even when w~am'. 
By (3.11) and (3.15), 


Bg=(3-) 1#[LOOE ES exp iw. a) Vd xa) 
: aa ~(1—g)V'(x—a,+.aa)] dx. 
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Since, for small w(<a-), 


PA exp (iw. a)=a-* | ...exp (iw. a,) da, 
it follows that we may rewrite the above equation for B, as 


B =(s,) 78" | lecoPw exp (iw. X) = [gV(x— X) 
—(1—g)V'(x— X+aa)]d®xd?X. . . . (7-3) 


Since, in the atomic state ¢, the electron is localized within a distance 
<a<w, we may put 

| d(x) |? =6(x) ; po eee 
in our calculation of B,. Also, for small w, we may calculate B,, on the 
basis of formulae for V and V’, which are valid when | x— X | is large. 
Such formulae are 


Lee 
Cy ee oe 
V(x ) [x] 
and 51 ies he Say 
Le," 
5a > 4 _ 0 ; 
V(x + aa) iexeena| 


where + Ze, are the changes per positive and negative ion, respectively, 
—é€) being the electronic charge. Inserting the expressions for |4|?, 
V and V’, given by (7.4) and (7.5), into (7.3), we obtain 


2 \"? 4rre?Zi a7: 
By= Gr) Ae [g+(1—g) exp (taw . a)]. 


For wa <1, this reduces to 
2 \'2 47e,?Zi 
By, = (— : ; av ate ds ce ore 
me (F -) varw 2) 
In accordance with the policy stated above, we shall assume this formula 


to be approximately correct for all the optical modes. Consequently, by 
(7.2) and (7.6), 


872Z2e,4 
il eee OE Well 
” NhM ,v?at : Rect 
where 
sin? Law 
= ee er ee RS 
Oe ap ve 


Since, by (2.19), |,aw,|<4a in the occupied region, it follows that 
(3aw,)?> sin? saw, > (aw,/7)? 


in that region. Consequently, it follows that a lower limit to S,, giving 
a correct order of magnitude to that quantity, will be obtained if we replace 
sin® $aw, by (aw,/7)? in (7.8). Thus, 
ow? - 
Si& >’ = 8,5 say | a eC) 


7 we 
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By symmetry, S,' is independent of J and is therefore equal to 


Hence, by (7.9), 


The summation in this equation is taken over half of the occupied region 
in w-space, and therefore 


Sim 2 
372 
Inserting this value for S, into (7.7), we obtain 
16Z7e,4 
or 0 
"= 35M ak en 
Typical values for M,, a, v, Z are 
M, S 10° Mey, 
a=5 x10 cm, 
(7.11) 


p=2x 10" sec}, 
Z=1. 


The lattice distance a, given here, is about ten Bohr radii, which is 
sufficiently large for the tight-binding approach to be justified, provided 
that the ionic radii <a. 

On using the values given by (7.11), eqn. (7.10) yields 


yo ~ 60. ies ere Be (7519) 


Hence, in view of (7.1), (6.11) can be satisfied for reasonable values of 
M,, a, v,Z. This confirms our assertion that the polaron effective mass, 
given by (6.9), can be extremely large in cases where the tight-binding 
approach is valid. 

In the above calculations, the ultra-violet polarization, due to the 
deformation of the ions, has been ignored. We can take this polarization 
into account by dividing the potentials V and V’, given by (7.5), by a 
dielectric constant K. XK is in fact the high-frequency dielectric constant, 
K.,, since the ultra-violet polarization can follow the field due to the 
electron adiabatically. It is easily seen that the division of V and V’ 
by K., leads to the division of our expression for y,° by K,,. Conse- 
quently, instead of (7.12) we have 

y° ~ 60/K ,,”. 
With K,,=2 (a typical value for a polar crystal), we then have y»,°~ 15. 
This value for y,° is still large enough for its exponential to be extra- 


ordinarily large. Therefore our main thesis, that the polaron effective 
mass may be extraordinarily large in cases where the tight-binding 
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approach is valid, is seen to still hold good, even when the screening effects 
due to the ultra-violet polarization are taken into account. 


APPENDIX 1 


We may rewrite eqn. (6.4) as 


fue SP LEAD Nidrgtieg) CFD Unie a) (A1.1) 
oi Splexp (—BA\,;)] ot) eae 
where 
Ang= (Guo Lop lille tate eee 


It may be seen from (2.22) that H,, is unchanged by the canonical trans- 
formation 
Nwj > Pe ats wh a ae Neil Yes 
Consequently, we may replace Hy,,, qwj, Ew; bY Hwy, YwiEwis — Nw5lYw9> 
respectively, in (Al.1). Hence 
=a Sp [exp (Aw j¥wjEws) exp ( —BH,,;)] 

Test Sploxp (—BHw,)I ; n(n een 
The eigenstates and eigenvalues of H,,, are O(n,;, €y;) and (nw;+ s)hvw;, 
respectively, where the 6’s are real, normalized functions. Therefore, if 
we drop the suffixes on 4, v, € and n, we may rewrite (A1.3) as 


fuu= | e€)exp (ieé)dé, . . . . (ALA) 
acre > [A(n, €) 2 exp (— nBhv) 
ple = a eee 
> exp (—nBhv) 


n=0 
In order to evaluate p(£), we make use of the well known relations (ef. 
Bohm 1954, p. 304), for the eigenstates of a harmonic oscillator, namely 


(1 =) 6(n, £) = [2h(n-+1)P0(n +1, é) 
and | (A1.6) 


(rns )e (n, €) = [2hn]"20(n —1, €). 


It now follows from (A1.5), (A1.6) that 


(8-5 nse \P =O > (n+p! “exp (—mPhv)O(n, €)0(n+1, €) — (A1.7) 
and 


(+5 —h ra p=0 > ms exp (—nBhv)O(n—1, €)0(n, €), (A1.8) 


where 


C= (enye] exp ce nBho) | 
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Equation (A1.8) may be rewritten as 


rs) co 

(« +H =) p=CSnt2 exp (—npliv)6(n—1, £)6(n, &), 
n=1 

since the summand is zero for n=0. On replacing (n—1) by n in this 

last equation, we obtain 


20 — 
(»¢+ Mh x) p=C > (n+1)!? exp [—(n+ 1)Bhv]O(n, £)0(n +1, €). 


n=0 


Consequently, by (A1.7), 


0 
| (+H 5) P= ex (—Bito)(ve— ah ae 
i.e. 
Op vv “ 
ae = ee tanh (SBhv)p, 
from which it follows that 
P=pp OXpP ls “tanh (26h») |, SS eee ee 


where py is independent of € and A. It now follows from (A1.4) and (A1.9) 
that 

fwy=A exp [—1hvA® coth (4ffiv)], . . . (A1.10) 
where 


A=p, [- _ exp | = ; tanh (yphione? | dé. 


A is therefore independent of A. By eqn. (A1.3), fy,=1 when A=0. 
It follows from (A1.10) that d=1. Hence, on reasserting the suffixes to 
v and A, we may rewrite (A1.10) as 

fugt=exp [— hry ;A,;? coth (FBhv,,;) ]. oy eee CA LST) 
Equation (6.5) follows immediately on substituting into (Al.11) the 
value for X,; given by (A1.1). 


APPENDIX 2 
In order to prove eqns. (5.6), we note that, by (5.3), 
W+(a,) =W(—a), fo ee eee 


since ¢(k) is real. 
By eqn. (4.13), W°(a,) is real, since the atomic wave-function ¢ is real, 


and, by (6.3) and (6.5), F(a,) is real. Hence, by (5.4), W(a,) is real, so that 
W*(a,) = W(a,). 
Equation (5.6) now follows from this last eqn. and (2.1). 


APPENDIX 3 
We wish to show that, for given k, the mean square deviation from 
«(k) of the energies ¢,,,, at temperature 7’, is negligible. This deviation 
is the difference between the thermal average, taken over lattice states, of 


1378 G. L. Sewell on 


€,,2, and the square of the average of «,,. It may be seen from (4.14) 
and (4.15) that this deviation is negligible if the mean square deviation of 
F,, is negligible, where 


F,,=(@,1|TotT;| On): pa, a. ae aie eae 
This latter mean square deviation is 
Asia ea eee 


where 
>F,,” exp (— BE, (lat)) 
Fr — a ee fi a 1 2, F e . A3.3 
= "Yexp(—FRE, (lat)? wae 
n 


so that, by (5.5), (A3.1) and (A3.3), 
F=F(a). ok Ree 


It remains for us to prove that A is negligibly small. 
By eqn. (A3.2), 
A Ff? 
FF 
and hence, since F?< 1, as may be seen from (6.3), (6.5) and (A3.4), it 
follows that 


|e 


F2 
<= -1; <6 gat ieey Seen 


In order to evaluate the r.h.s. of this inequality, we note that 
(O(n, €)| exp (iAn)| O(n, £)) = (8*(n, €)| cos An] O(n, €)), (A3.6) 


since the expectation value, in the state 6(n, €), of any odd power of the 
momentum 7 is zero. It now follows from (2.25), (3.17), (A1.2), (A3.1) 
and (A3.6) that 

F,=(0,* |G] ,,), Y ot ee i Sea) 


where 


4 
é! = a Il COs ded ihwh: . . . . . (A3.8) 
j=1 w 


The quantity 
(One| G?| ©) Ey (O74 OF) 


is the mean square deviation of G, for the state ©,,, and is therefore positive. 
Consequently, by (A3.7), 


He = (OF 67107). recs SB) 
It now follows from eqns. (2.21), (A3.3), (A3.7) and (A3.9) that 


4 
H= 11 [1 av(costis ye.) a) eee oa) 
and ‘ 


4 
He < |] OL 'av(cos? ie) ae TD 
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where, for any function g(y), 
S ' ee 
av (9(w;)) _ 'PL9(7w;) exp ( BH,,;)] : 


y A312 
Splexp(—PHw,)] iy 
Consequently, by (A3.5), (A3.10) and (A8.11), 

Le A< TY Ty tes Baste) (A3.13) 


eye 2av (cos AW ;7Nw;) 


Moreover, by eqns. (Al.1) and (A3.6), fw is the thermal average of 
cos (A,,;7~;), for the oscillator H,,;, and therefore, by (A3.12) 


> 


i wit = @V(cos Awe de 


Hence, by (A1.11), 


av(cos A, 7w;)=exp (—$o,,;), ee eC OsL eS) 
where 


Ow; = Hivy;d%y; coth (4Bhivy,;). Jel » eae 


It follows from (A3.15) that we may replace A, o by 2A, 4c, respectively, 
in (A3.14). Hence 


av(cos 2 dw Tw3) = exp (—2c,,;); 


and therefore, by (A3.13), (A3.14), 


4 
ce le COs Ov ee Cen eo. 10} 


j=l w 
In order to simplify this inequality, we note that, since ¢> tanh ¢, for t>0, 


0< | (¢—tanh t) dt=407—In cosh o; 
0 


and therefore, cosh o<exp (40”). Inserting this inequality into (A3.16), 


A<exp (A)-—1, 7 coer ee CAS La) 
where 


Tt follows from this last equation that, since the o’s are positive, 
RES le lage te eee eS 1) 


where 
Pea (oo Gu) + (Cyst Gway es ee (ane, 19) 


Since 7 is an even function of w, as may be seen from (3.14), (3.15), (A1.2), 
(A3.15) and (A1.19), we may replace >’ by 3> in (A3.18). Hence, 


A<}>7w- ee ere, (43.20) 


On expressing 7, in terms of B,, C,,, v and c, by means of (2.24), (3.14), 
(A1.2), (A3.15) and (A3.19), and on replacing 


Na 3 
> eee esas by Gap | Ow. 
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(A3.20) yields 


A< NI, ee eT 
where 
ee aml tl B,,|*v2 coth? (4Bhv) + |C,, |*c?® coth? (Bhew) | sin’ w . a, déw. 
(A3.22) 


The integration in this last equation is taken over the whole of the occupied 
region of w-space. It follows from (3.15) and (A3.22) that I is inde- 
pendent of N. Moreover, the contributions to the integral in (A3.22), 
due to the long waves, do not lead to a divergence. This may easily be 
verified on evaluating B, and C\, for the long waves, using the method 
employed in § 7 to evaluate B,. J is therefore finite. 

It follows from (A3.21) that 


INS) BI) INNS ge 


and therefore, by (A3.17), A+0as Noo. Consequently, A is negligibly 
small for a macroscopic body. 
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ABSTRACT 


The thermal expansion of a solid is due to anharmonicity of the interatomic 
forces. It may therefore be calculated, on the basis of a continuum model, 
from the observed elastic anharmonicity, i.e. the third-order elastic con- 
stants. Using Lazarus’ (1949) measurements on the pressure variation of 
the elastic moduli of several cubic crystals, the Griineisen y is calculated for 
the limiting cases of low and high temperatures. For KCl, NaCl the high 
temperature values agree with those obtained from thermal expansivities, and 
a decrease in y at low temperatures is predicted. However, Lazarus’ data on 
metals cannot be reconciled with the thermal expansions, which suggests that 
the continuum model is inapplicable to these cases. 


§ 1. INTRODUCTION 


THE assumptions of the Debye theory of specific heats are that the normal 
modes of vibration of the solid are independent, and that the frequencies 
of the modes, up to the limiting frequency, are those of a perfectly elastic, 
isotropic continuum. In this theory everything depends on one para- 
meter, the Debye temperature ©, which may be calculated from a know- 
ledge of the atomic volume and macroscopic elastic constants of the solid. 
This theory is inadequate to describe such phenomena as thermal 
expansion and the lattice thermal conductivity which depend on the 
anharmonic nature of the interatomic forces. Calculations of the lattice 
conductivity have been carried out (Leibfried and Schlémann 1954, 
Klemens 1956) and have relied upon a single phenomenological parameter 
as an average measure of the anharmonicity, namely the Griineisen constant 
y (Griineisen 1926). This is defined experimentally byy= a V/ XpCy, where 
a is the thermal expansivity, y, the isothermal compressibility and C,, the 
heat capacity at constant volume /. . 
In the Griineisen theory y is a measure of the change of the lattice 
frequencies with the volume of the crystal. This may be described in the 
continuum model by the addition of terms of the third-order in the strain 
components to the elastic strain energy. Since elastic anharmonicity has 
been studied experimentally it is interesting to see whether from such 
information we can predict the observed values of y obtained from thermal 


expansion measurements. 
+ Communicated by Dr. J. M. Ziman. 
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Bridgman’s measurements on the change of compressibility with 
pressure, were used in this way by Slater (1939), and fair agreement between 
the calculated and experimental values of y was obtained. However, 
more detailed information on elastic anharmonicity has been given by 
Lazarus (1949) who measured the pressure variation of the individual 
elastic moduli of several cubic crystals: in this paper we employ these 
data to calculate y on the basis of the anisotropic continuum model. It 
must be emphasized that the values obtained have only the limited validity 
of the continuum model. 

We shall first derive the theoretical formulae needed for the calculation. 


§2. APPROXIMATE THEORETICAL FORMULAE 


Consider only monatomic, cubic crystals: the normal modes of vibration 
may be labelled by a wave vector q and polarization index p; p=1, 2, 3 
corresponding to one longitudinal and two transverse branches of the 
frequency spectrum. By treating the crystal as an assembly of loosely 
coupled harmonic oscillators, whose frequencies depend on the volume of 
the crystal, and applying statistical mechanics, it is easy to calculate the 
free energy F from which the Griineisen parameter follows as a consequence 
of the thermodynamic relation 


oe | 
Xp aan) 


The result is (Slater 1939, Barron 1955): 


ae Veneer ys Cara . . . - . . (1) 
q,P q.Pp 
where 
Cup = Se ire ee Ves 
; kT /} {exp (hw, »/kT)— 1}? 
and 


W,,, is the frequency of a normal mode and the sums are over the 3N 
possible modes (N is the number of unit cells). Thus + is expressed as a 
weighted average of the individual y, ,,’s. The weighting factor C , is, of 
course, the contribution of the (q,p)th mode to the heat copocte 
Cy= > Cun: 


q,p 
The simplifying property of a continuum is the absence of dispersion of 
the elastic waves. The frequencies are given by w,,»=48,(9,¢) where 


q=|q| and s,,(@, 4) is the sound velocity of the pth branch in the direction 
specified by polar angles 0,4. Hence 


ye ding _ (Ce ”) 
T 


dln V ~ @nV 
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The wave vector q, according to periodic boundary conditions applied 
toa cube of side L, has components 27n,/L (integer values of n,); thus for 
a particular mode, the n; are fixed and qc V-8, We can write 


Vq.p =Yp(9, ), where 
_1_ /dlns,(9, 4) 
Yp(9, $)= 3 (Se), . fA on eae We) 


Equation (1) can be considerably simplified for the limiting cases of 
low and high temperatures. 


For low temperatures, 7 <@, converting the sums to integrals over 
q-space and using standard Debye theory, we obtain 


z ¥p(9, ) / | dQ 
ini | dQ VLR at ORS 
=Je,6 oR =) Rear se 
where dQ is an element of solid angle in the direction 8, ¢. 


For high temperatures, 7'> ©, the heat capacity per mode Cg p =* and 
hence 


which, in the continuum approximation becomes 


y= aad | rol? $40. he VAR are) 
p 


The experimental information we shall use refers to room temperatures, 
for which T~©. Writing C(x)=C,,,/k=2e*/(e” —1)?, and remembering 
the limiting frequency for the longitudinal branch is about twice that for 
a transverse branch, we have at room temperature 


Oat 


DOne et Cmax 0-5 for transverse branch, 
kT ‘a kO ~ \ 1-5 for longitudinal branch. 


Hence 1 > C(x) > ae that is, the weighting factor for the higher 


frequencies of either branch has not fallen significantly below the upper 
limit, unity. Equation (6) is therefore a good approximation at room 


temperatures. 
The difference between the low and high temperature eqns. (5) and (6) 


respectively, lies simply in the type of average taken over the y, (9,4). 
These limiting values of y will be denoted by y,, Yq. 


§ 3. ANHARMONICITY IN AN Enastic ConTINUUM 
Before explicitly evaluating y,(0,¢) we shall consider briefly how the 
anharmonicity is specified in terms of the third-order constants, and the 
available experimental information concerning them. 
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The free energy per unit (unstrained) volume may be written as an 
expansion in the Lagrangian strain components 7,, at a fixed temperature : 


f= : > Cog NpNgt a Opel airs ‘ 7 ‘ % : (7) 
psa pqs 

where p,q,7r take values 1, 2,3...6 (Birch 1947). The second-order 
coefficients are the elastic moduli of Voigt, and the C,,,, are the third-order 
or anharmonic coefficients. For cubic crystals belonging to the classes 
of maximum symmetry there are only three independent elastic constants 
C115 C19» Cqq and six independent third-order constants Cy4,, Cie, Cra Cres» 
C93, C56: 

Since cubic symmetry is preserved under uniform pressure the form of 
the relations between an additional infinitesimal stress and corresponding 
infinitesimal strain are also preserved, but the ‘effective elastic constants’ 
depend on the pressure. Birch (1947) showed that the elastic moduli 
under a finite pressure, corresponding to volume /’, are 


AV ‘ 
Cyy(V) = yy + > (244 + 2C yy + 604) + 4C 439), | 


NY ee | 
Cyo(V) = Cy + 37. (Chog + 4C 42 — C11 — C2); «cogs pees) 
OV 
Caa( V) = Cq4 + 37 (Cag + Cyq + 2Cy9 + 3C p44 + Cree), | 
‘0 


to the first order in AV=V—V,. At zero pressure, V = V4, these reduce 
tO Cy, C2, Cag, 1-€. the coefficients in eqn. (7). 

Lazarus (1949) has determined the adiabatic elastic constants of several 
cubic crystals by measuring the velocities of ultrasonic waves up to 
pressures of 10000 bars, and Hearmon (1953), by plotting the dependence 
upon volume, deduced values for the quantities 


Cg =6C y+ 4C ye, Cy=Crog +402, Ca= 4 ast Crean 


using the above relations. It is clear that in the derivation of the latter 
equations the c,,(V) etc. are taken to be the isothermal moduli and the 
volume change AV is carried out isothermally. Since Hearmon omitted 
to convert Lazarus’ adiabatic moduli to the corresponding isothermal 
moduli before applying eqns. (8), his values for the third-order constants 
are incorrect. The errors are not likely to be great; in the case of Cu the 
corrections to the values of C,, C,, C4 are estimated to be less than Loa 

However, for our calculation of the Griineisen ratio we require the 
dependence of sound velocities upon isothermal volume changes and this 
will be given correctly by substituting Lazarus’ values for the adiabatic 
moduli (distinguished by a superscript s) and Hearmon’s values for OC. Cee 
Cq (also distinguished by a superscript s) into the equations : 


a AV 
C1 °(V) =e 8+ BV, 


(2¢,,°+ 2c,9°+C,,), ete, . . . (10) 
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§4. Toe Maenirupe or y AT Room TEMPERATURE 
Previous estimations of y (Slater 1939) neglected crystal anisotropy 
and relied on special assumptions, for example that Poisson’s ratio is 
independent of volume. We can avoid these approximations by obtaining 
a general expression for y,(4,4) as follows. 
The secular equation giving the sound velocity s(0,¢) in a cubic con- 
tinuum may be written (Quimby and Sutton 1953): 


228+ BLT, (0,¢)2—B,T,(0,6)=0; . . . . (12) 


the variable z is related to s by 


z=ps?/K,—K, hooey a aw Pane are 7 areas MO.) 
and 
B,=1-K,?, B,=1—3K,?+2K,3, 
: é : (13) 
Ky=Cyy'/Ks, Ky=(6y.°+Cy4°)/Ky, Kg= yy? — Cy’, 


I’ (@, 6) =sin? @(sin? 6 sin? d cos? ¢ + cos? 6), 
[,(@, 6) =sin* 6 cos? @ sin? d cos? d. 


The elastic constants appearing here refer to the particular pressure at 
which we require the sound velocities; they are the c,,°(V) of eqn. (10). 
Differentiating eqns. (11) and (12) with respect to InV at constant 
temperature and eliminating (02/0 In V),, immediately gives (dIns/dIn V),, ; 
> € In s(@, 2) (6, 05,—268;, 15) NG. K 
T 


Pe eee le ee ee (1 
olnV (822-2+B,V,)(2+K,) z+Kk, K; ee 
where we have used the notation F’(V)= (AF (V)/aln V),,, for any function 
of V. 

The quantities K,’, B,/ can be expressed in terms of the elastic moduli 
and their variation with volume. For example: 


Cy = $3 (2C 41° + Wey +C,%), at V=Vo, 


Keg = (647819 — 41° Cyy® + Cre*Caa™ — Cra" Cag? + C11 °Cqq® — Cyy" Cag®)/K 5’, + (15) 
1 = —6K,K,.. | 
For special directions eqn. (11) can easily be solved for the sound velo- 
cities (de Launay 1956) whence eqn. (14) gives the values of y,, e.g. in 
the [100] direction 
84=(C5/p)? 3 S2,3= (Caa®/p)*? ; 
(16) 
dlns ae ; Tier bs 
Ga =4(e8 es +1); (Cag /Caa® +1), twice. | 
ih 
To calculate y for a given substance we take from experiment the values 
of p, 641°, C19°, Cag*, Cy°, Cy, Ca®- Equations (12), (13) and (15) at once 
give us Cy" - _. K,...B,, K,'....B,'. Fora particular direction 6, ¢ the 
secular equation is solved, the three roots <, give the sound velocities s,, 


5 A 
P.M, ve 
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and substituting in equs. (14) and (4) the values of y,, for this direction 
are obtained. Numerical integration according to low or high tem- 
perature equations yields y. 

The integrals were evaluated numerically on EDSAC 2, the computer at 
the University Mathematical Laboratory, a Gauss twelve-point inte- 
gration formula being used over the range 0,¢=0° to 90°. It was found 
that eqn. (14) is indeterminate for the longitudinal branch (p=1) in the 
[111] direction (it gives 0/0), but this does not invalidate the integration 
procedure as the points at which the integrand is evaluated do not approach 
closer than 6° to this special direction. 

Since in the course of the calculation a mean sound velocity s,, and 
Debye temperature © are evaluated, they are recorded below for their own 
intrinsic interest. Their definitions are 


dL Sell AS | pw age tA) 8 eek: 
Sn? Vaated ene) k V 
(n is the number of atoms per unit cell, V the number of unit cells.) 

To see the contributions from the individual polarization branches to 
the final average y and also the effect of anisotropy of the medium, the 
sound velocities s, and Griineisen constants y, are recorded for the [100], 
[110] and [111] directions using (16) and similar equations. 

The experimental values and data from which they were evaluated are 
included in table 1. 


Table 1 


we 
RSs 
Si 


Material ¢ ie x = 


Co p s 
(10~° °x-+) |(cal. mole? °K-1)| (em) |(cm? dyne 10-12) ta x5C 


Cu i673 5:848¢° ola! 0-716 2-03 
Al 2-368 5-817" 9-99 1-285 2°26 
KCl 3-67> 12-21 37-54 5-450! 1-48 
NaCl 3°98°¢ 11-868 27:04 4-076! 1-60 
CuZn 2-004 11-6]h 15-54 0-861! 2:22 


a, = coefficient of linear expansion | at 298°K 
Osa capacity per mole at constant pressure | and 
=volume per mole atmospheric 
Xs=adiabatic compressibility J pressure. 


Bijl and Pullan (1955). 

Dayal (1944). 

Buffington and Lattimer (1926). 
International Critical Tables. 
Giauque and Meads (1941). 
Berg and Morrison (1957) 
Clusius et al. (1949). 
Sykes and Wilkinson (1937), 
Lazarus (1949). 


OFS Cee Bee. 
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Table 2 
© ee A S,: LL ae i e) 
Material | Direction } 82: T) 83: Ty, : L| yy: T rl Vy Ya - (es 
100 | 4-37] 2-88] 2-88]1-69| 0-591 0-59 ; 
Cu 110 | 5-00 | 2-88] 1-62} 1-41} 0-59] 0-80 | 0-762 | 0-937 | 2-587 | 331-7 
lll 5-18] 2-12} 2-12] 1-40] 0-69] 0-69 
100 | 6-25 | 3-25] 3-25] 2-84] 257] 92-57 
Al 110 | 6-41 | 3-25] 3-93] 2-28] 2-57] 5-48| 3-77 | 3-26 | 3-430 | 400-2 
ne 6-57 | 2-96] 2-96] 2-12] 4.23] 4.93 “fl 
a 100 | 4-54] 1-78 | 1-78 | 2-73 | —0-87 | —0-87 
KCl 110 | 3-90 | 1-78 | 2-92 | 1-81 | —0-87| 3-03. | 0-525 | 1-37 | 2-377 | 225-4 
lll 3-68 | 2-60] 2-60/ 1-35] 2-41] 2-41 
100 | 4-80] 2-44] 2-44] 2-81] 0-14] 0-14 
NaCl 110 | 4-52 | 2-44] 2-92] 1-99] 0-14| 2-891 1-28 | 1-60 | 2-900 | 307-0 
lll | 4-40] 2-76 | 2-76] 1-70| 2-18| 92-18 
100 | 3-94] 3-15] 3-15]1-90] 1-69| 1-69 
CuZn 110 =| 4-93 | 3-15 | 3-42] 1-82| 1-69] 1-7511-74 | 1-76 | 2-201 | 280-3 
111 5-22 | 2-02] 2-02] 1-80] 1-70} 1-70 


L: longitudinal branch; 7’: transverse branch. 
{ In 10° cm sec}, 


The © values agree well with other approximate calculations (Blackman 
1955) but the values of y are significantly different. For comparison we 
have calculated y using Slater’s theory which, in addition to the approxi- 
mations already mentioned, neglects the difference between the adiabatic 
and isothermal compressibilities. The resulting formula is 

Pe OIE: bigest TS) 
where a, 6 are the coefficients in Bridgman’s expression 
—AV/V,=aP—bP?, T constant. 

The coefficient b represents some average measure of the anharmonicity 
of the solid. But by using directly the results of Lazarus the adiabatic 
modulus can be inserted into the theory; taking soc(y,p)*” gives 

nels (5) wall Soares (ose ee 19) 

6, 2, eln V ji 6 2\ 642+ 2075 

In table 3, y is computed from eqns. (18) and (19) using Lazarus’ data 
including the values of a and b he records for his five materials. 


Table 3 
: 2 b 1 ae 1 /dOlnys 
Material as =f a 6 olan Vie 
Cu 1-8: D2 
Al 1:86 223) 
KCl 1:58 Dail 
NaCl 2-00 3-07 
CuZn Poste) 2-25 


5A 
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§ 5, Discussion OF RESULTS 


When arranged according to increasing experimental y’s (table 4) it can 
be seen that the theoretical y,, yq also exhibit a general increase except 
for the case of Cu. As expected the actual numerical agreement is much 
better for y, than for y,, and for the alkali halides is well within the limits 
imposed by experimental error in the original data. But for Cu and Al 
some other explanation of the numerical discrepancy must be sought. 


Table 4 


NaCl 


Yexpt 1-60 


1-23 
1-60 


It appears from table 2 that it is the contributions from the transverse 
branches which cause the final values of y, (and y,) to differ from the 
observed y by so great an amount. Yet this does not offer an explanation 
of the discrepancies because the different values of y associated with the 
polarization branches are still unexplained. Indeed, the continuum 
model does not attempt to explain them ; it simply takes these values from 
direct observation of the low frequency end of the lattice spectrum. 

Assuming the basic theory leading to eqn. (1) is sound, there are two 
possible causes of the discrepancies: errors in the data or inadequacies of 
the model. The elastic moduli are themselves accurate to within a few 
per cent but errors of 10 or 20° may reasonably be allowed in the third- 
order constants. Granting this, a discrepancy of a factor of 2 for the 
case of Cu still cannot be explained. The more likely explanation lies in 
the basic assumption of our model which allows us to replace the y,, ,’s by 
the values at the long wavelength end of the spectrum. It would hardly 
be surprising if, in fact y, , were very different for short waves thus causing 
the average y to differ markedly from that predicted by the continuum 
model. 

The humiliating fact is that Slater’s calculation, based upon the crudest 
possible model, gives better overall agreement with experiment, especially 
when one considers the assumptions in its derivation. Indeed, if one of 
his false assumptions is corrected, giving the modified Slater values of 
eqn. (19), the agreement is made worse. It is not difficult to see why these 
modified values differ from y,,: they depend on the anharmonicity only 
through the volume variation of the adiabatic compressibility ; it is effect- 
ively the average value of y for compressive (longitudinal) waves that we are 
calculating as may be seen by comparing with the column y, in table 2. 

From our results some tentative predictions as to the temperature 
variation of y can be made. The Griineisen parameter is a function of 
temperature (at say, constant pressure, to conform with experimental 
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measurements), this dependence being due to two distinct causes. The 
Yq,p 8 themselves may depend on temperature, for the elastic moduli and 
third-order constants are temperature dependent, but this effect is likely 
to be small. Of greater importance is the temperature dependence of the 
weighting factors C,,,, of eqn. (1). This causes the weighted average of 
the y, ,'s to vary: at low temperatures, owing to the 1/s? factor, y is close 
to the value for the transverse branches, while at high temperatures the 
unweighted average is the appropriate one. If this is the main cause of 
the temperature variation the low temperature limit of y should be close 
to our y;. Thus, for KCl and NaCl y should decrease with decreasing 
temperature, the behaviour being more marked for KCl. Unfortunately, 
experiments are only available for the metals Cu and Al for which a slight 
decrease has been observed below 40°K (Bijl and Pullan 1955, Simmons 
and Balluffi 1957): our values of y, are too unreliable for any comparison 
to be made here. 

Finally we remark that the electron gas in a metal only contributes to 
the Griineisen parameter at temperatures low enough for the electronic 
and lattice specific heats to be comparable (Bijl and Pullan 1955): eqn. (1) 
need therefore only be corrected for temperatures below a few degrees 
Kelvin and the influence upon the room temperature y’s, with which this 
paper is concerned, is negligible. 
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Note added in proof.—Since this paper was written information on the 
pressure variation of the elastic constants of Cu, Ag and Au has been 
published (Daniels and Smith 1958). For Cu the pressure derivative of 
the shear modulus, ¢,,°, is larger by a factor of 3 than the earlier 


Material Yu Yexpt 
Cu 1-97 2-03 
Ag 2-40 2-48 
Au 3-02 2-90 


measurement and considerably alters the calculated Griineisen parameter. 
In fact the values of y,, derived from these data ave in good agreement 
with the experimental values for all three metals indicating that the 
continuum model is of wider applicability than suggested by Lazarus 


measurements. 
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ABSTRACT 

The magnetic moment of a multidomain grain at its blocking temperature 
is determined by the domain arrangement which gives the minimum magnetic 
energy; for a large grain the energy minimum is produced by a self- 
demagnetizing field equal and opposite to the external field. Below the 
blocking temperature the domain configuration is frozen in and the moment 
increases with the increasing spontaneous magnetization. This results in a 
TRM proportional to small inducing fields, as in the theory of single domain 
grains and as required for agreement with experimental observations. 
Artificial TRM’s of 17 rocks for which analyses of magnetic minerals are 
available show good agreement with calculated values. The theory also 
gives a formula for the maximum possible initial susceptibility of rocks 
which agrees well with experimental observations. 


$1. IyrRopUCTION 


In the process of cooling of igneous rocks, grains of magnetic mineral are 
separated, generally as magnetite or titanomagnetite. The ferro-magnetic 
character of these grains appears at their Curie points, several hundred 
degrees below the solidification temperatures of the rocks, but although 
spontaneous magnetization occurs at these temperatures, as is evident 
from saturation magnetization measurements, the property of remanence 
or ability to retain induced moments is not acquired until temperatures 
several tens of degrees lower still are reached (Nagata 1953). 

Thermal agitation is evidently responsible for preventing remanence 
above these lower temperatures, which have been termed ‘ blocking 
temperatures ’ (Néel 1955). The moment induced in a grain at its 
blocking temperature is ‘frozen in’ by further cooling and becomes 
thermo-remanent magnetization (TRM). This is important in rock 
magnetism because of its stability or resistance to demagnetization. 

The general properties of TRM, including its proportionality to the 
inducing field, and the occurrence of blocking temperatures are well 
explained by the theory of Néel (1955), in which single domain grains of 
magnetic mineral are considered. However, Néel recognized that most 
naturally occurring magnetic grains were much too large to be correctly 
considered as single domains. The purpose of the present paper is to 
present a simple theory of TRM in multidomain grains, which has a direct 
and quantitative application to the magnetism of igneous rocks. 


+ Communicated by the Author. 
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For Néel’s (1955) mechanism of thermal agitation of single domain 
erains to restrict the natural moments of igneous rocks to small fractions 
of their saturation moments we require 


emg] ae er es 


where H is the external field imposed on a grain which has a moment p: at 
its blocking temperature 7’, and k is Boltzmann's constant. Except in 
the special case of a grain for which 7’, is so close to the Curie temperature 
that the spontaneous magnetization has only a very small fraction of its 
value at ordinary temperatures, formula (1) can apply only to grains of 
the order 0-01 micron in diameter, whereas grains from one to 200 microns 
appear to be most important in the magnetism of igneous rocks (Nagata 
1953). If such large grains were single domains they would acquire 
saturated TRM in very small fields, and the fact that this is not observed 
is particular evidence that each grain contains a number of domains. 


§2. TRM ty LarRGEe GRAINS 


The total magnetic energy of a grain may be divided into the following 
six parts (Bozorth 1951, p. 811): 


(i) exchange energy ; 
(ii) erystal anisotropy energy ; 
(iii) magnetic strain energy ; 
(iv) domain wall energy ; 
(v) magnetostatic energy (self-demagnetization due to magnetic poles) ; 


(vi) mutual energy between magnetic moment of grain and external 
field. 


We are interested in small fields, in which the moment is well below 
saturation and the domains are magnetized in easy directions, all of which 
have equal anisotropy energy in the absence of strain (Bozorth 1951, 
p. 478), so that (i) and (ii) are independent of the state of magnetization 
of the grain and can be neglected for the present purpose. In a large grain 
in its unmagnetized state the domains are so arranged that there are 
virtually no surfaces of magnetic polarity, either on the surface of the 
grain or internally (Bozorth 1951, p. 834), so that (v) is minimized. There 
is a further requirement that (iii) plus (iv) must have the minimum value, 
but if the grain is much larger than the limiting size for single domains, 
(v) exercises a much stronger control on domain geometry than (iii) and 
(iv) because the volumes occupied by closure domains and domain walls 
are small fractions of the whole grain. Therefore only (v) and (vi) need 
be considered for a first approximation to the domain structure of the 
grain. 

In an external field, H,,, surfaces of magnetic polarity are induced on 
the surface of the grain, but not internally. This is equivalent to giving 
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the grain a moment J per unit volume, and results in a self -demagnetizing 
field ee NI), N being the demagnetizing factor, a geometrical quantity 
explicitly stated only for ellipsoidal grains. In this state the grain 
experiences an internal field H,, given by : 


Fiver moi Ge ne ee re) 
Contributions (v) and (vi) to the magnetic energy give a total value 
E=—H,I+4NIT, oe ee re), 


where the first term is (vi) and the second term is (v). If the grain. is 
above its blocking temperature, the domain readjustments in the field H z 
are not subject to internal restraints, so that J follows H,, reversibly and 
is self-determined to give # the minimum value. This is found by putting 


dE 
= —H, 7 Ni=—H;=0, 
whence 
H 
a PPO fo ee Ae 
That this corresponds to a minimum of His shown by a second differentiation 
fad Dy 
ie a 


which is necessarily positive. 

As an alternative derivation of eqn. (4) we may assume that, for small 
external fields, the permeability of the material of the grain is effectively 
infinite above its blocking temperature, which is equivalent to the assump - 
tion that there are no restraints to domain readjustment, and implies 
directly that the internal field H; is zero. 

As long as H, is small, so that J is well below saturation, the equilibrium 
described by eqn. (4) is maintained during the cooling of the grain until 
the blocking temperature is reached and the relaxation time for the 
decay of moments not in equilibrium with H,, becomes long compared with 
the time scale of cooling. Below this temperature the moment increases 
in the same way as the spontaneous magnetization, o, so that at 
laboratory temperatures the thermo-remanent moment is 

Hy oR 

Cre eae (5) 

where (o,/o,) represents the ratio of spontaneous magnetizations at room 
temperature and the blocking temperature. 

As an approximation to the TRM of the magnetic minerals in rocks, 
eqn. (5) gives a reasonably good fit to experimental results, but detailed 
comparison will be left until a refinement has been added in § 4. 


§ 3. SuscepTipitity of Rocks 


The susceptibility of rocks at laboratory temperatures will now be 


considered. This requires a re-examination of the term ‘ blocking 
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temperature ’ which has been adopted from Néel’s (1955) theory of TRM 
in single domain grains. 

The moment of a single domain grain can change direction under the 
influence of thermal agitation, but only in discrete jumps because it is 
constrained to lie in one of a few directions by shape and crystal aniso- 
tropies. The blocking temperature is the temperature at which thermal 
agitation provides just sufficient energy for the moment to cross the 
potential barriers between these directions. Below this temperature the 
probability that a change in domain direction will occur becomes so small 
that in effect the moment is completely frozen in. ‘Therefore a single 
domain grain below its blocking temperature cannot contribute to the 
low field susceptibility of a rock. 

The susceptibility of a multidomain grain at low temperatures is not 
zero because a small moment may be induced in the grain by a continuous 
and reversible process of domain boundary displacements and not in 
large, discontinuous jumps. This is the process used to explain the 
initial susceptibility of large ferromagnetic bodies (Bozorth 1951, p. 480). 
Its reversibility arises from internal restoring forces which are due to 
magnetic energy contributions (iii) and (iv) (§2). These are constraints to 
domain readjustment which prevent the material of the grain from having 
infinite susceptibility. At high temperatures thermal agitation overcomes 
these constraints and allows the domain geometry to readjust itself to 
give the minimum total magnetic energy. The blocking temperature of 
such a grain may then be defined as the temperature at which internal 
constraints become effective in limiting the adjustment of domain 
boundaries when the grain is subjected to small fields. 

Accurate estimates of the susceptibilities of rocks cannot be made 
because the internal forces limiting domain readjustment are not well 
understood. However, the maximum possible susceptibility of a rock 
may be estimated by the following argument. The condition for this 
maximum to be observed is that internal constraints have only a 
negligible effect in limiting domain wall movement, and the moment 
induced in a grain has that value which results in a minimum of the total 
magnetic energy, being given by eqn. (4). Putting N = (47) and consider- 
ing a rock with M% by weight of magnetic minerals of density 5-25 g/em?, 
we have for the moment per cm? of the rock : 


RS 
Hy 4m 100 5-25’ 
where 6 is its density. The maximum possible mass susceptibility of the 
rock is then given by 
MaX ¥y = —— = 0:-45x 103M, aes eS Cn) 
Nagata (1953, p.93) showed that the maximum initial susceptibility of 
igneous rocks is given approximately by 


max y)= 4:5 x 10-°0j,,, 
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where Nagata’s y, is the same as y,, in the present paper and Cy, is the 
fraction by weight of normative magnetite in a rock. The agreement with 
eqn. (6) is remarkably good if M is the percentage of normative magnetite 
and not the total iron—titanium oxide content. 

More useful is a comparison of eqn. (6) with the maximum susceptibilities 
of rocks observed at elevated temperatures, at which the internal con- 
straints to domain wall movement are least effective. This comparison is 
made in table 1, using data of Nagata (1953) taken from his table 3-IV. 
The calculated values of maximum susceptibility have been obtained, 
using normative magnetite contents for values of M in eqn. (6). It is not 
clear why this gives better agreement than the total Fe—-Ti oxide content 
which indicates a maximum initial mass susceptibility of the order 6 x 103 
for the rocks listed. 


Table 1. Calculated and Observed Maximum Susceptibilities of some 
Japanese Rocks. Specimen Numbers and Experimental Values 
are from Nagata (1953) 


eee (max yy) cale. x 10% (max yy) obs. x 10% 
ry 3-54 2-56 
18 Led 2-06 
19 2-40 2°75 
20 2-92 1-42 
21 orks 1-93 
2 1-77 1-08 
23 1-98 1-46 
24 ()-62 0-68 
25 1-15 2-05 
28 1-15 1-42 

1-98 1-60 
as 1:77 0-79 
33 0:83 O-15 
24 1-04 0-96 
35 1:35 0-42 
37 3-02 3°34 
38 1:76 2-06 
56 2-29 0-87 
9] 1-67 2-02 


; , OF [ BY - AGNETIZING FIELD 
§ 4. REDUCTION OF TRM sy THE SELF-DEMAG 


The self-demagnetizing field exerted on a grain by its moment L o per 
unit volume is (—NJI,), so that the value of Iepy calculated in §2 is 
reduced to J, by this field, where 


Toe —Lo=N Loxy 
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yy being the volume susceptibility of the material of the grain at 
laboratory temperatures. The observed TRM, I,, is thus given by 


vi 
aa) SEEM aR a ee 
1+Nxv 
Now we can estimate y, for the material of the grains in a rock from the 
observed susceptibility of the whole rock. The moment / induced in a 
grain by an external field H,, is given by 


0 


gee: 

Tee CE ee ees 
If x, is the mass susceptibility of the rock, density 6, with M%, by weight 
of magnetic minerals of density 6,, then 


fe xv (001M) 


MM 3Ha Lp Nyy iu 


Bh 


which gives 
1 


= SSS SSS eee’ 8 
os (0:01.M/x 8) —N Ae) 

Substituting (8) into (7) we obtain 
Jone ee te 


10= THOM x8) IN) 


The value of J, for the whole rock is obtained by multiplying eqn. (4) 
by (0-014 /5,,;) so that eqn. (9) becomes: 


(Hy/N)(ox/oR)0-0112 (8/0 yy) 


whence 
La ~*~ On 0-O1M | - 


Kquation (10) expresses the TRM per gram of rock per oersted in terms 
of the weight percentage M of magnetic mineral of density 5,,, and yy 
is the mass susceptibility of the rock. 


§ 5. COMPARISON WITH EXPERIMENT 


Equation (10) gives the required proportionality of TRM and field and 
also allows a reasonably accurate estimate to be made of the constant of 
proportionality. Using data of Nagata (1953), the average blocking 
ES aT TOS Rae of an assembly of grains is usually 60°C + 20°C below the 
ee ae so ce Bese (n/n) e 0-3 for a blocking 

» point of a mineral which appears to be 
nearly pure magnetite. This ratio will be subject to some variation but 
the value 0-3 is probably accurate within a factor of two for most igneous 
rocks. It is to be expected that the magnetic grains in igneous rocks 
have randomly oriented shape anisotropies, so that N=(47) as for 
spherical grains will be taken as an approximate effective value this is 
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reconsidered in §7. The value of 84, is taken as 5:25 g/cm, as for pure 
magnetite (Nagata 1953); the densities of other magnetic minerals do not 
depart very widely from that of magnetite. Substitution of these values 
into eqn. (10) reduces it to: 

og Bs 103M —: 

Hy, = ¥518% Soo va ee) 


5 
Values of (J)/SH,,) calculated from eqn. (11) are compared with 
experimental values in table 2, using data of Nagata (1953, tables 3-I, 
4-V, 4-VI), whose specimen numbers are given in the first column. 
(1,/SH,,) corresponds to Jy, and yy, to y) in Nagata’s tables. 


Table 2. Calculated and Observed TRM’s of some Japanese Rocks. 
Specimen Numbers and Experimental Values are from Nagata, (1953) 


Se % Magnetic 
men ae Xu X 108 |(19/SH_) cale. x 109\(L,/SH_) obs. x 108 
number 

17 13-45 (7-87) 2-25 12-9 (4:4) 12-1 
18 11-89 (3-94) 1-31 13-7 (1-6) 10-7 
19 15-54 (5-33) 1-82 16-5 (2-0) on 
20 13-17 (6-48) 1-28 15-7 (5-6) 10-5 
21 11-88 (6-95) 1-55 12-8 (5-4) 13-0 
22 12-44 (3-94) 0-70 16-5 (2-6) 41-7 
23 11-48 (4-40) 1-29 13-1 (2-4) 5:2 
24 6-84 (1-39) 0-40 9-0 (0-8) A) 
25 11-60 (2-55) 0-90 14-6 (0-9) 14-8 
27 11-19 (4-17) 1-28 12-7 (2-1) ont 
38 14-00 (3-91) 1-01 17-9 (2-6) 27:6 
52 10-76 (4-63) 1-09 12-7 (3-4) 2-8 
53 8:78 (3-94) 0-57 11-4 (4-1) 9-0 

-22 (3-47) 0-54 13-7 (3-5) 4-6 

‘64 (6-02) 0-78 13-5 (6-5) 4-8 

‘91 (6-25) 1-37 18-0 (4-9) aya 

24 (5-09) 0-73 16-1 (5-3) 4-6 


+ Figures are total Fe-Ti oxide content by weight; values obtained using 
normative magnetite are given in parentheses. 


It is not immediately obvious whether the total iron—titanium oxide 
content or the normative magnetite gives the more appropriate value of 
M. If the whole Fe-Ti content uccurs in strongly magnetic grains then 
the former is correct. However, in many rocks some of the Fe-Ti 
content occurs in non-magnetic minerals which do not contribute to the 
TRM. and for such rocks the normative magnetite content may give a 
better value of M. Columns 4 and 5 of table 2 indicate that in ten of the 
rocks listed the whole Fe—Ti content is magnetically active, but that in 
the other seven a separation of magnetite has occurred, 
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§6. Decay or TRM with TIME 

A matter of some geophysical interest arises from the theory which has 
been presented. The evident decay over geological time of the natural 
remanent moments of rocks has been regarded (Nagata 1953, p. 140) as a 
decay of natural TRM, which at ordinary temperatures has a very long 
relaxation time. It is imagined that thermal agitation causes a gradual 
disruption of TRM. If this were the case then it would be exactly the 
same process as that by which TRM is acquired, except that it takes a 
very much longer time, and therefore remanent moments with the stability 
of TRM would be gradually induced by the changing geomagnetic field 
at the same time as the original moments decayed. Rocks with decayed 
natural moments would then be quite unsuitable for palaeomagnetism. 

However, the TRM of multidomain grains causes self-demagnetizing 
fields within the grains which are several times stronger than the field 
in which the TRM is induced. These self-demagnetizing fields must be 
the principle cause of decay of TRM. They are exactly opposite to the 
moments of the grains and the decay of the moments therefore involves 
no changes in their directions. Since the influence of a changing geomag- 
netic field is only secondary, rocks with decayed natural moments are not 
necessarily unreliable for palaeomagnetic work. 


§ 7. Discussion 


Further consideration will be given to the following approximations or 
assumptions which have been made : 


(1) multidomain grains have blocking temperatures ; 

(2) the susceptibility of the material of a grain above its blocking 
temperature is so high that the self-demagnetizing field provides 
the only limit to the magnetic moment ; 

(3) an average demagnetizing factor (<7) can be used for all grains; 

(4) grain interactions can be neglected ; 

(5) grains of magnetic mineral are much larger than the limiting 
size for single domains ; 

(6) the ratio (o,/o,) is roughly constant for all rocks; 

(7) chemical analysis indicates the content of magnetic minerals. 


(1) A definition of the term ‘blocking temperature’ as applied to 
multidomain grains was suggested in $3, but it requires qualification, 
since it is not a single temperature but a temperature range. The crystal 
imperfections and internal strains which become effective in opposing 
changes in domain geometry at low temperatures are unlikely to be 
distributed so uniformly through a grain that movements of all domains 
are opposed equally. It is more probable that some domain boundaries 
become frozen in at higher temperatures than others within the same grain ; 
the maximum spread can be inferred from the sharpness of the blocking 
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temperature of a sample of artificial magnetite powder (Akimoto, quoted by 
Nagata 1953, p. 161) of uniform grain size, and appears to be of the order 
100°c. This is not inconsistent with an average blocking temperature 
60°C below the Curie point as assumed in §5, so that the estimated value 
of (,/0,,) is probably not seriously in error on this account. 

(2) This assumption has been discussed by Néel (1955) who concluded 
that “the apparent susceptibility of rocks with coercive forces of less than 
500 oersteds... must depend more on the demagnetizing field effects 
than on the intrinsic susceptibility of the material”: coercive forces are 
certainly less than 500 oersteds at elevated temperatures. 

(3) The value (47) for the average demagnetizing factor in the field 
direction of a random assembly of grains is only an approximation and in 
general a somewhat smaller value must be expected. This can be seen 
by considering two very elongated grains whose major axes make equal, 
but opposite angles, #, to the external field, H. Each grain acquires a 
TRM almost at an angle 6 and given by 


H cos@ oR 


sf ——— 
TRM N’ oR ? 


where N’ is the small demagnetizing factor of each grain. The net 
moment in the field direction due to the TRM of many such pairs of grains 
is then 


7/2, 3 
= | cos? 6 sin 6 dé, 
QO 


(sin @) being the weighting factor for the number of grains at an angle 0. 


Thus 
il ff OR 


Fomxt = 3 NY" Gy 


where NV’ may be much smaller than (37). 

Generally dimension ratios greater than 3 or 4 are not to be expected ; 
if all the grains in a rock are prolate ellipsoids with the dimension ratio 
3-5, then the TRM parallel and perpendicular to the major axes can be 
calculated as above and the total TRM is 50% higher than that indicated 

5 (bs 
ee e os magnetic interactions between grains have been neglected. 
These interactions were considered by Néel who showed that under certain 
conditions they could lead to a mechanism of self reversal of PRM. How- 
_ ever, one of the conditions was a higher concentration of magnetic grains 
than is normally found in rocks. Assuming that the magnetic grains are 
uniformly dispersed in a rock, interactions between them cause a reduction 
in the effective self-demagnetizing fields of the individual grains. For 


spherical grains this reduction can be expressed in terms of an additional 


magnetizing field H’ given by 
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where J’ is the moment per unit volume of the rock. If J) is the moment 
which the same grains would give if they were very widely separated then 


I’ =1,+5xyH', 
5 and yy being the density and mass susceptibility of the rock. Thus 
A 1 
Ly 1—(37)dyyq 
if yy,~ 10-%, as for the rocks listed in table 1, (I’/I,) ~ 1-01, so that unless 
the grains are arranged in a grossly non-uniform way this effect is negligible. 

(5) It must be noted that the theory presented in §§2—-5 applies only 
to large grains. Small grains, containing only a few domains, for which the 
absence of internal surfaces of magnetic polarity cannot necessarily be 
assumed or for which domain wall and magnetic strain energy cannot be 
neglected, must retain some of the behaviour of single domains, and there- 
fore become magnetized to an extent greater than is indicated by eqn. (11). 
This probably applies to grains smaller than one micron in diameter. 

(6) Assuming that the difference between the blocking temperature and 
the Curie point has been correctly estimated, the only variation in (c,/c,) 
arises from the differences between the shapes of spontaneous magnetiza- 
tion versus temperature curves. As an extreme case we may consider a 
grain with a blocking temperature of 520°c and spontaneous magnetization 
decreasing linearly between 20°c and the Curie point at 580°c. Such 
a grain would have (cg/o,)~ 9-3 compared with the value 1-9 estimated 
from the curve for nickel. These represent extremes and it is more 
reasonable to expect 


(7) Magnetic analysis would be more appropriate than chemical analysis 
since both the magnetic minerals and their concentrations can be estimated 
from saturation magnetization versus temperature curves; also (c/o) 
could be found separately for each constituent. Nagata (1953, p. 107) has 
found that the saturation magnetization of rocks is roughly proportional 
to their normative magnetite, but that this does not represent the actual 
mineral content is evident from the constant of proportionality, which 
is only half the saturation magnetization of pure magnetite. Therefore 
the separation of magnetite referred to in §5 must be a separation of 
titanomagnetite with saturation magnetization roughly half that of pure 
magnetite. The saturation magnetization data would be satisfied by a 
separation of pure magnetite to the extent of half the normative magnetite 
but according to eqn. (11) this would give a TRM much smaller than is 
observed. Although no conclusion can be drawn, the results of table 2 are 
consistent with a separation of titanomagnetite having about 50 mol% 
Tike,O, particularly in specimens 27, 52, 93 to 96. Nagata (1953) ae 
that numbers 27, 52, 53 and 93 behaved irreversibly when heated, probably 
on account of a separation process, but that numbers 17 to 25 did not. 
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ABSTRACT 


The increase in adhesion obtained for aluminium films by the pre-deposition 
of films of chromium is investigated, using Heavens’ method of adhesion 
testing. It is found that the adhesion varies with the thickness of the 
chromium film, and with time. The variation with time suggests that the 
nature of the increase in adhesion is primarily one of alloying at the metallic 
interface, and subsequent age-hardening. 


§ 1. INTRODUCTION 


VACUUM-DEPOSITED films of aluminium on glass are extensively used as 
front surface mirrors for optical purposes and it has been known for some 
years that the durability of these mirrors could be considerably increased 
by depositing a thin layer of chromium on the glass immediately before 
the aluminium (Holland 1956), both evaporations being performed in the 
same vacuum chamber without admitting air. The reflectivity of the 
aluminium film, which is relatively thick, is unaffected ; but the adhesion 
to the glass is improved and it has been suggested that the chromium 
forms a barrier, preventing attack of the aluminium by the alkali in the 
glass (Holland 1956) with consequent loosening of the bonding between 
the aluminium. Such a screening effect would only prevent deterioration 
of the reflecting film and cannot be regarded as a complete explanation 
since an increase in adhesion is actually observed. A preliminary 
investigation into the nature of this increase in adhesion by Heavens (1950) 
showed that large increases in adhesion could be obtained when the 
thickness of the chromium deposit was greater than about 3004. It was 
later suggested (Heavens and Collins 1952) that this was due to an 
oriented overgrowth mechanism but it appears from the present work 
that the observed increases are largely due to a diffusion—precipitation 
phenomenon at the interface between the aluminium and the chromium. 

It is difficult to measure the adhesion of thin films and almost all 
accepted methods fail when the adhesion is strong. The results of 
abrasion testing (Townsley 1945) are not reliable since a burnishing action 
is involved, which affects results. Consequently, in this work, the 
adhesion of the evaporated films has been measured using the method 
adopted by Heavens. A chrome-steel point having a smoothly rounded 
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tip was drawn across the film and the load on the point was graduall 
increased until a critical load was reached at which tl : z 
S the film was completely 

removed from the glass substrate leaving a clear channel in the film. 
This stage was readily determined by viewing the scratches under a low 
power microscope using transmitted light. At the critical load a broad 
evenly illuminated channel with sharply defined edges could be observed 
whereas lower loads did not remove the metal completely and a scratch 
was produced which appeared bright at the centre, tapering to complete 
opacity at the edges, indicating a channel of more circular cross section. 
The results obtained by this method were readily repeatable and an 
accuracy of +5°% could be obtained in determining the end-point. 

Heavens did not give any analysis of his method but the application of 
the method to single metal films on various substrates has since been 
investigated by Weaver and Benjamin (1958) who showed that it gives an 
effective measure of adhesion. They have shown that the shearing force 
at the interface is a direct function of the vertical load on the point and 
that calculated results are in good agreement with theoretical estimates. 
The validity of the method for double films such as aluminium upon 
chromium has not, however, been established and will be discussed later 
in the light of experimental results. 


§ 2. EXPERIMENTAL 

Specimens were prepared on glass microscope slides by evaporating 
chromium and aluminium in succession so that each metal extended over 
two-thirds of the length of the slide, overlapping in the centre. Thus the 
centre section of the slide was covered first with a chromium film and then 
by an aluminium film while the two end sections were covered with pure 
chromium and pure aluminium respectively. The evaporations were 
performed in a glass bell-jar evacuated by an oil diffusion pump, which 
produced pressures of 2-5x 10->mm Hg as measured by an ionization 
gauge. The slides were supported in a rotating jig, about 20cm above the 
heaters which were made of tungsten wire, 1mm diameter, wound in the 
form of spirals. The slides were cleaned carefully before placing them in 
the chamber, being scrubbed first with Teepol solution to remove dirt and 
grease and then rinsed under running water before drying on a linen cloth 
and finishing with lens tissue. One position of the rotating jig allowed the 
full surface of the slide to be exposed to a cleaning discharge in the bell -jar 
for about 15 minutes during the pumping process. When a sufficiently 
low pressure had been reached the jig was rotated so as to bring the slide 
first over the chromium heater and later, over the aluminium heater. 
Suitable masks were used so that only two-thirds of the slide near one 
end were coated at a time. The rates of evaporation were approximately 
104/sec for the chromium and 20A/sec for the aluminium, the two 
evaporations being performed in rapid succession. 7” 

The thicknesses of the chromium films were monitored during deposition 
and measured in vacuo by using a parallel] beam of white light and 
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measuring the transmission. The light beam was projected into the 
vacuum chamber and reflected by a glass prism, screened from the 
evaporated material, so that it passed through the specimen slide and 
struck a barrier layer photocell connected to a sensitive galvanometer. A 
calibration curve of film thickness against transmission was obtained by 
weighing a number of specimen slides, covered with chromium only, on a 
microbalance immediately after deposition. The thickness was calculated 
from the bulk density so that the measurements gave an effective weight 
thickness. The aluminium films were evaporated to be just opaque, 
corresponding to a thickness of about 600-800 A. 

Immediately after preparation of each specimen it was taken from 
the vacuum chamber and measurements of adhesion were made. Values 
were obtained for aluminium, chromium and for the double film of 
chromium overlaid by aluminium on the same slide and prepared under 
the same conditions. From the relative values the effect of the chromium 
underiayer could be seen. It was known that commercially produced 
films were aged after deposition and it was noticed during the investigation 
that the adhesion values obtained changed with time. Measurements 
were therefore made periodically until steady readings were obtained and 
ageing was complete. For films kept at room temperature this ageing 
process took about three weeks. In accordance with commercial practice 
some of the specimens were aged in ovens at a temperature of about 120°C 
and, in this case, the ageing process was completed in about three days. 


§ 3. RESULTS 


Typical ageing curves obtained for one specimen which had been aged 
at room temperature are shown in fig. 1. The initial adhesion figures 
obtained both for the pure metals and for the double film were all low. 
The adhesion of the aluminium film increased gradually with time but 
quickly reached a steady but low value. The final value obtained did not 
vary noticeably with thickness or from one specimen to another. The 
adhesion of the pure chromium also increased steadily but more rapidly 
at first and over a somewhat longer time, finally reaching a steady value. 
This value varied slightly with the thickness of the specimen, being 
somewhat greater for thicker films. The significant feature of these 
ageing curves is, however, the difference between the chromium curve and 
the curve obtained for the double film of aluminium overlying chromium. 
Chromium undoubtedly sticks to glass better than aluminium and this 
could partially explain the increased adhesion of the double film, as 
compared with aluminium, but the difference between the two ageing 
curves must still be explained. The difference is even more pronounced 
for specimens aged at 120°C and typical curves are shown in figs. 2 (a) 
and 2(b). The final steady values obtained for pure aluminium and pure 
chromium are comparable with the values obtained when ageing at room 
temperature. ‘The final values reached by the double films are, however, 
much greater and the shape of the ageing curve has changed, 
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The final steady value of the adhesion of a double film after ageing 
varies considerably with the thickness of the chromium underlayer ; this 
is shown in fig. 3, where results obtained from a large number of specimens 
are given. This graph does not compare very well with the results 
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obtained by Ifeavens, as shown in fig. 4. His results compare more 
favourably with fig. 5 in which we have plotted the adhesion values 
obtained immediately after evaporation, against the thickness of the 
chromium underlayer. This suggests that Heavens did not take ageing 
effects into account but measured his specimens immediately after 
removal from the vacuum chamber. 
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§ 4. Discussion 


The increased strength of the double metal film cannot be explained in 
terms of any change in the basic adhesion forces between the chromium 
and the glass. The same chromium films act as a single metal film over 
one section of the slide and as an underlayer for another section. This 
implies that the difference between the measurements on the two sections 
is due to the aluminium overlayer and for the aluminium to affect the 
interface between the chromium and the glass, aluminium atoms would 
have to diffuse through the chromium film. This would be a slow process 
compared with the rapid initial rise of the ageing curves, even if we 
assumed that the aluminium atoms at the interface could increase the 
adhesion. This seems most unlikely since our measurements show that 
chromium adheres to glass better than aluminium. 

The ageing process for the double films at room temperature normally 
appears to consist of two stages occurring at different times. Similar 
processes have been observed in the age hardening of alloys and our 
curves correspond closely with those of Gayler (1937) who attributed the 
two stages of ageing to an interdiffusion of the two metals followed by a 
partial precipitation of an intermetallic compound. According to the 
theories of age hardening (Smith 1949), the interdiffusion would cause a 
hardening or increased resistance to deformation with no change in lattice 
parameters and this would be followed by a release of strain and 
precipitation of an intermetallic compound. At higher temperatures the 
processes would be accelerated and could overlap, giving one continuous 
change such as is observed in our graphs of high temperature ageing. 

There appears to be some doubt about the phase diagram of the 
aluminium—chromium system (Metals Handbook 1948), and information 
is not complete but it has been claimed that there are nine distinct 
intermetallic compounds (Bradley and Lu 1937). Chromium does, 
however, have strong similarities to nickel and cobalt and complete phase 
diagrams are available for aluminium—nickel and aluminium—cobalt. In 
each case the phase diagram has a pronounced peak indicating formation 
of a stable intermetallic compound of high melting-point. These 
compounds have a high energy of formation and are known to be hard 
and strong. It would therefore not be unreasonable to expect that at 
least one of the chromium—aluminium compounds will have similar 
properties. This leads us to suggest that, when a film of aluminium is 
deposited on top of a film of chromium, an intermetallic reaction occurs 
at the interface and a compound precipitates forming an intermediate 
layer. Such a process could explain the results which we have obtained 
since precipitation of a hard intermetallic compound would explain the 
second step in the room temperature ageing curves. In the method of 
adhesion measurement which has been used the effect of a strong 
intermediate layer between the aluminium and the chromium would be 
to spread the concentrated load upon the point over a wider area so that 
a greater force would be required to produce the same shearing effect at 
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the interface between the glass and the chromium. At the same time 
the bonding between the aluminium and the chromium would be improved. 
This implies, of course, that although this method of measuring adhesion 
has been shown to be suitable for single metal films and capable of 
theoretical interpretation, this is no longer true for double films such as 
are considered here. Valuable information may nevertheless be drawn 
from the results provided care is taken in the interpretation. 


$5. Execrron DirrrRacTION INVESTIGATION 


An electron diffraction investigation was carried out on some specimen 
films. These were prepared in the same apparatus as the adhesion slides, 
but the metals were condensed on to collodion films supported on copper 
grids. The table of the apparatus was designed to give films of the three 
types already investigated, i.e. both the pure metals on separate grids and 
a composite film on a third. An electron microscope was used to obtain 
the transmission electron diffraction patterns. Before the films were 
examined they were allowed to age for about 300 hr at room temperature. 

The first set of grids examined consisted of a chromium film approxi- 
mately 3004 thick, and an aluminium film 5004 thick. Both of the single 
metal films showed the characteristic diffraction pattern of the bulk 
materials, fig. 6(a,b), Pl. 98. The compound film, fig. 6(c), Pl. 98, gave 
a pattern that was simply a summation of those obtained for the single 
metals. No trace of intermetallic compounding was detected. If, 
however, the diffusion rates were low, and the degree of penetration small 
then diffraction effects from such a layer might be obscured by the 
relatively thicker metallic layers. 

Another set of films was thus prepared using much thinner films, the 
thickness of the chromium layer being 304, and of the aluminium 50-754. 
After the ageing period, the single aluminium film showed only the 
presence of an amorphous substance, fig. 7 (a), Pl. 98, probably aluminium 
oxide. No trace of the metal itself was obtained. This is in agreement 
with the theories of oxide growth, since it has been estimated that the 
thickness of an oxide layer on aluminium would achieve the value of 
approximately 354 in 300 hours at room temperature (Mott 1947). The 
chromium pattern, fig. 7(b), Pl. 98, also showed the presence of oxidation, 
but some of the metal has remained unoxidized. 

For this specimen, the compound metal film gave a pattern completely 
uncharacteristic of any of the previous specimens, fig. 7(c), Pl. 98. The 
rings are broad but distinct, representing a face-centred lattice of unit 
cube cell and spacing 2-924. This corresponds to the « phase of the 
aluminium—chromium alloy system (Bradley and Lu 1937) and appears 
for a concentration of approximately 87°, of chromium by atomic volume. 
This is not the most prominent phase in the system, but when it is 
considered that the structures of the other phases are based upon the 
stacking of three of the a phase unit cells on top of each other, and due 
to the direction of interdiffusion that the electron beam would more 
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likely traverse this system parallel to the ‘c’ axis, then we can see that 
it is doubtful whether more information could have been obtained. 


§ 6. ConcLuUsIon 


The evidence deduced by electron diffraction appears to give reasonable 
confirmation of the general validity of our interpretation of the ageing 
curves which have been obtained for double films of aluminium on 
chromium. Such values as we have measured for the adhesion of these 
films immediately after deposition are in general agreement with the 
results given by Heavens in that films of chromium thicker than about 
300A cause a marked increase in the adhesion of aluminium films. But 
even thinner chromium films cause some increase in adhesion according 
to our measurements although not so great as that produced by thicker 
chromium films. Heavens found that thin chromium films did not have 
any appreciable effect and we can only suggest that this may be due to 
differences in cleaning methods or evaporation conditions whereby our 
chromium films had a greater basic adhesion than those prepared by 
Heavens. 

Heavens also indicated that in a case where air was admitted to the 
vacuum chamber after evaporating the chromium and before evacuating 
again to deposit the aluminium, no increase in adhesion was obtained. 
Such a result is to be expected if our theory is correct, since an oxide 
film would form rapidly on the chromium when exposed to air and such 
an oxide film would at least hinder if not prevent an intermetallic 
reaction such as we have suggested. 

The existence of a critical thickness for these chromium underlayers, 
cannot readily be explained from the results. Previous work on the 
optical properties of thin films of chromium (Hill and Weaver 1958) has, 
however, shown that thin films of chromium are composed of aggregates 
which may be regarded as oblate spheroids with the minor axis 
perpendicular to the substrate. Each aggregate is surrounded by a very 
thin layer of oxide when examined in air, and the size of the aggregates 
changes with increasing film thickness, the axial ratio becoming greater 
so that the aggregates tend to become platelets. This type of structure 
exists up to weight thicknesses of about 100-1204, above which there is 
a transition region and the structure changes to a more continuous type 
corresponding to metallic chromium. It appears quite feasible that films 
below about 200-3004 have not reached the stage of having a continuous 
metallic structure even in their upper layers, whereas it has been shown 
that the upper layers of films of about 506A weight thickness have the 
same refractive index as bulk chromium metal. In thin films the 
screening of the aggregates by a layer of oxide would prevent an inter- 
metallic reaction with a superimposed aluminium film such as has been 
suggested and, consequently, these films would have little effect in 
increasing adhesion. . cn . 

Finally, we ought to point out that since obtaining the results given it 
has been found that the basic adhesion of the chromium films themselves 
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can be increased by using a shutter in front of the source and moving it 
only when the source is at operating temperature so that only high 
velocity chromium particles strike the substrate to form a deposit. Under 
these conditions the overall adhesion of a double film is increased and the 
effect of the intermetallic reaction becomes less noticeable but it 
nevertheless still persists and plays some part as can be seen from the 
ageing curves of fig. 8. 


Fig. 8 
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ABSTRACT 
Observations have been made of slip, grain-boundary migration and grain- 


boundary extrusion during the fatigue of high purity lead in vacuo (approx. 
5x 10-3 mm Hg). 


§ 1. INTRODUCTION 


Recentiy, Forsyth (1957 a, b) has reported slip-band extrusion during 
fatigue of an aluminium +7-5% zine +2-5% magnesium alloy, and 
Forsyth and Stubbington (1957) have observed grain-boundary extrusion 
during elevated temperature (250°c) fatigue of an aluminium + 4°, 
copper alloy. In addition to slip-band extrusion, Cottrell and Hull (1957) 
have detected slip-band intrusion in copper fatigued at room and low 
temperatures. Extrusion and intrusion effects have been suggested, in 
theories by Cottrell and Hull (1957) and by Mott (1958), to be intimately 
associated with the start of fatigue cracks. 

The following observations describe grain-boundary extrusions and 
deformation during the fatigue of lead in vacuo at room temperature. 


§ 2. EXPERIMENTAL 
Flat cantilever specimens with tapered sides were prepared from lead 
strip of 99-999% purity. The specimens were annealed at 100°c, 
chemically polished and tested at a small cyclic strain (estimated to 
be 0:05%) in reverse bending fatigue in a vacuum of the order of 
5x 10-3 mm Hg at a rate of 500 c.p.m. 


§ 3. RESULTS 
Observations on the polished surface showed that the initial effect of 
cyclic straining caused slip and some boundaries to migrate. After 
about 104 to 10° cycles, the rate of migration slowed down and almost 
ceased. The initial and actual positions of a migrated boundary are 
shown in fig. 1¢ and are indicated by A and B respectively. After the 
rate of migration had diminished, specimens were lightly repolished to 
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reveal more clearly the subsequent changes when given further cyclic 
straining. Figure 2 shows the same area as in fig. 1 after repolishing at 
7:5 x 105 cycles with 1-25 x 108 additional cycles, i.e. a total of 2-0 x 10° 
cycles. It is seen that the grain boundaries appear darkened and show 
little sign of further migration. 

After repolishing and further cyclic straining slip tended to be denser 
in the central region of grains than in zones along either side of boundaries, 
fig. 2. 

Examination of the grain-boundary zones at higher magnification 
showed that two types of extrusions with a metallic lustre had occurred 
along some boundaries. A whisker-like extrusion is shown in fig. 3 at 
8 x 10° cycles and fig. 4 indicates the subsequent development of further 
whiskers, marked by arrows, in the same area after 2 x 10° cycles. The 
size of the whiskers varied and some were found up to 50, long, with 
a diameter of approximately 1. Occasional ones had a fine structure 
—appearing twisted or kinked along their length. 

A more continuous type of boundary extrusion is illustrated in fig. 5. 
This type often occurred together with the whisker-like extrusions. 

In addition, some boundaries showed what appeared to be localized 
slip patterns which had geometrical forms, fig. 6. These were found on 
one or both sides of boundaries and their sides were frequently parallel 
to traces of slip systems in other parts of the same grain. 

Further work is in progress. 
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ABSTRACT 


A description is given of tests in which a } in. or 3 in. diameter weighbar, 
seated vertically on a thin steel plate, was impacted by a finite mass. The 
force applied to the plate and the principal stresses produced at the bottom of 
the plate have been derived from the output of resistance strain gauges. 
The stress-time curves at the centre of the bottom surface of the plate, 
obtained at a number of impact velocities, have been compared with stress— 
time curves derived from theoretical considerations. The agreement 
between the theoretical and experimental results is found to be exceedingly 
good. 

In such tests the plate showed no sign of having yielded even though the 
stresses reached were much greater than the static upper yield stress of the 
steel in the plate. This is attributed to the fact that strain rates of the order 
of 100 see- were attained in the plate in these tests. 

The manner in which the disturbance created in this type of test was pro- 
pagated radially has been determined from the way the principal stresses, at 
the bottom of the plate, varied with radius and time. It is found that there 
was an annular area of compression which propagated radially in advance of a 
central area of tension. 

Tests have also been made in order to determine the manner in which a 
disturbance caused by an impulsive point load propagates radially in a plate. 
A dural and a steel plate were impacted by a high-tensile steel sphere and the 
principal stresses produced at the bottom surface of the plates were once 
again determined from the output of strain gauges. It is found that the 
disturbances produced in these tests propagated in a similar manner to the 
way the disturbances propagated in the tests in which the force was applied 
over a finite area. 


§ 1. IyrRopvUcTION 


investigation has been made in order to obtain information 


concerning the magnitude of the elastic stresses produced in a thin plate 
by a normal impulsive force. Tests are described in which a zin. or 2in, 
diameter weighbar, seated vertically on a thin steel plate, was impacted 
by a finite mass, thus applying an impulsive load to the plate. The 
force-time curve applied to the plate in such a test has been obtained 
from the output of resistance strain gauges mounted on. the bottom of the 
weighbar. Also the radial and circumferential stresses on the bottom of 
the plate have been derived from the output of strain gauges fixed to the 


plate surface. 


T Communicated by the Author, 
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The stress—time curves at the centre of the bottom surface of the plate, 
one for each impact velocity concerned, have been compared with 
stress-time curves derived from a theory presented by Sneddon (1945). 
Sneddon considered the bending of a plate caused by the application of a 
constant force uniformly distributed over a circle of radius a and acting 
for a finite time 7’. He derived expressions for the principal curvatures 
of the centre plane of the plate from which the stresses throughout the 
plate could be calculated. The results of the present tests have been 
compared with the theory by assuming that a particular force-time curve 
applied to the plate can be approximated by a number of superimposed 
square waves. 

The manner in which the disturbance, created in such a test, has 
propagated radially has been investigated by plotting curves of radial 
and circumferential stresses against radius at a number of instants. Also, 
it was considered that a determination of the way in which a disturbance 
caused by an impulsive point load was propagated radially would be of 
interest. Hence tests have been made in which a dural and a steel plate 
have been impacted by a high-tensile steel sphere. The manner in which 
this type of disturbance has propagated in these plates has again been 
determined from the output of strain gauges fixed to the bottom surface 
of each plate. 


§ 2. APPARATUS 


Figure 1 is a diagrammatic representation of the impact apparatus 
used in the present investigation. The cylindrical high-tensile steel mass, 
which impacted the weighbar, was free to slide in the head of a hammer 
arrangement. The length of the hammer, measured from the pivot to the 
centre of the steel insert, was 54in. The dural hammer arm was of 
channel section and the head was machined from a cylindrical aluminium 
casting. The hammer was released from an electromagnet and arrested, 
after the steel insert had impacted the weighbar, by a rubber pad fixed to 
an aluminium angle which was supported at its ends so as to be above the 
plate. At the moment of impact the steel insert was horizontal and 
concentric with the weighbar. 

The plate used in such tests was 20-5in. in diameter and 0-285 in. thick, 
and was bolted round the circumference to a levelled base plate. Two 
3gin. long weighbars were used, one }in. and the other 2in. in diameter ; 
also the impact end of each weighbar was machined to a spherical shape 
of radius 10in. For the steel used in the weighbars and plate, the elastic 
modulus was taken as 30 x 106lb/in?, the Poisson ratio as 0-30 and the 
density as 490 Ib/ft?. 

Three wire resistance strain gauges were mounted on each weighbar. 
These gauges had a gauge length of }in. and were centred Lin. from the 
plate end of each weighbar. Also Sie gauges were fixed to the bottom 
surface of the plate at the centre and at radii of 1, 2, 3, 4, 5 and 6in. from 
the centre respectively. At each of these radii there were two gauges, 
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one in a radial direction and the other in a circumferential direction. Only 
one gauge was required at the centre of the plate, for at any given instant 
the radial and circumferential strains were equal at this point. 

The electronic apparatus used to record the outputs from these gauges 
was almost the same as that used in an earlier investigation (Campbell 
and Maiden 1957). In these previous tests the cathode ray oscilloscope 
spot was triggered and brightened by the magnetostrictive output of a 
coil surrounding a weighbar. In the present tests the spot was triggered 
and brightened by the gauge output itself. 


Fig. 1 


Electromagnet 
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Mechanical impact apparatus. 


§ 3. THEORETICAL SOLUTION FOR A PLATE IMPACTED BY A CONSTANT 
FORCE 

Let p, # and o refer to the density, elastic modulus and Poisson ratio 
respectively of the material in a plate of thickness 2h. Also let 
D=2Eh?/3(1 —o?) denote the flexural rigidity of the plate. 

Boussinesq (1885) has shown that if a normal force is applied at a point 
on a large elastic plate, then at any instant the velocity v of that point 
is proportional to the force f acting on the plate. He found that the 
relationship between these two quantities was given by 

if = Gphbt meee a9 ee ieee ce eeare) 
where b?= D/2ph. 

Sneddon has used this relationship in deriving, frem the differential 
equation governing the transverse displacement of a thin plate, 
expressions for the principal curvature at the centre of a plate when 
loaded by a constant force F applied over a circle of radius a and acting 
for atime 7. The expressions for the curvature « at a time ¢ are given by 

«= —[log,(a?/4bt)+y—1]F/16pahb® for O<t<T . . . (2) 
and k= —[log,(1—T/é)|F/16prhb? fort>T  . . . s+ (3) 
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where y is Euler’s constant. It is to be noted that this formula is accurate 


only if a?/4bt is small. 
Also, if M is defined as the bending moment per unit length at the centre 


of the plate then 
M = Dk(1+o) ee ee ee ec i, 


and the tensile stress p at the bottom surface of the plate can be 
determined from the formula 


p=1:5M |h?. oi AW d 2 eee eae 
Hence from eqns. (2), (3), (4) and (5) the stress p, at any time ¢, is given by 
p= — (log, (a?/4bt)+y—1]F(1+0)/5-33nh2 for O0<t<T . . (6) 

and 
p= —[log, (1—T/t)|F(1+0)/5-337h? for t> 7. oa bg ee 


This solution only applies when a plate is subjected to a force of 
constant magnitude, a condition which would not often be met in practice. 
In the next section it will be shown how the theory can be used when a 
force-time curve of general shape is applied to a plate. 


§ 4, EXPERIMENTAL RESULTS 


4.1. The Magnitude of the Tensile Stress Produced 
Immediately Below the Load 


The in. diameter weighbar, seated vertically on the plate, was 
impacted at four different hammer head velocities. At each impact 
velocity four strain-time records were taken from the weighbar gauges 
and four from the gauge at the centre of the plate; the weighbar records 
were found to be repeatable to within +3 and the records from the 
plate gauge to within + 2%. 

A typical plate strain-time oscillogram, with calibration lines super- 
imposed during enlarging, is presented in fig. 2(a), Pl. 101. All the traces 
were interrupted every two microseconds in order to give an accurate 
time base. 

The force-time curve applied to the plate, at a particular impact 
velocity, was determined by multiplying the relevant average weighbar 
strain-time record by za?H, where HF is the elastic modulus for the 
weighbar steel. The derived force-time curves, one for each of the four 
impact velocities concerned, are shown in fig. 3 (a). 

The stresstime curve at the centre of the bottom surface of the plate, 
for each impact velocity, was determined from the relevant average 
strain—time record using the formula relating the stress p and strain « at 
the centre of a plate. This formula is given by 


p—Hel(l—o). rn 


The stress-time curves obtained in this manner are shown in fig. 3(b), 
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Fig. 3 
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Force-time and stress—time curves obtained at various heights of drop using a 
4 in. diameter weighbar. 

(a) Force-time curves. Curves A, B, C and D correspond to heights of drop of 
37-4, 28-0, 22-8 and 10-2 in. respectively. 

(b) Stress-time curves. Curves A, B,C and D refer to experimental results 
and curves E, F, G and H to theoretical results at height of drop of 
37-4, 28-0, 22-8 and 10-2 in. respectively. 
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The force-time curves applied to the plate in these tests were by no 
means square waves. Hence, in order to compare Sneddon’s theory with 
the curves shown in fig. 3(b), each force-time curve has been considered. 
to be equivalent to a number of superimposed square waves. The 
magnitude of all these square waves, excluding the top increment, has 
been taken as 200lb. Curve D in fig. 3 (a) shows such an approximation. 
For each force-time curve eqns. (6) and (7) have been used to determine 
the effect of each increment of load on the stress at the centre of the plate. 
Thus, by adding all the increments of stress at corresponding times, the 
theoretical stress-time curve has been found at each impact velocity. 
The four theoretically derived stresstime curves are presented in fig. 3 (b) 
(broken curves), and it is seen that they compare very favourably with the 
experimental stress-time curves. It is to be noted that eqns. (6) and (7) 
only apply to a plate of infinite radius, whereas the plate used in the present 
investigation had a radius of 10-25in. In the next section it will be shown 
that, for the time interval over which the theory has been compared with 
the experimental results, the centre of the plate acted as though the plate 
was infinitely large. 

It is of interest to note that the thickness of the plate was chosen so 
that the plate should theoretically act as an extension of the weighbar. 
This has been arranged in the following manner. Provided there is no 
wave reflection at the weighbar-plate interface then the force f applied 
to the plate at any instant is given by pcvra?, where c is the velocity of 
sound in the weighbar. This force is also given by 16phbv, hence equating 
these two values for f we obtain 


16hb =cza? 
ie. h=a23sa(1—o2)"4/4-. 2 ww... (9) 


which for a $in. diameter weighbar yields h=0-1425in. That the actual 
reflected waves were small is seen from the face that the oscillations on 
the force-time curves are slight ; and also from the fact that the force-time 
curves die off approximately exponentially, this being a theoretical 
condition for a finite mass striking a semi-infinite rod. 

Tests were also made using a 3in. diameter weighbar and the results 
obtained have been compared with the theory. Only one impact velocity 
was considered, this corresponding to a 28in. height of drop. The 
average force-time curve applied to the plate is shown in fig. 4(a) and 
the average stress—time curve, at the centre of the bottom surface of the 
plate, is presented in fig. 4(b). Once again the force-time curve was 
considered to be composed of a number of square waves in order to derive, 
from eqns. (6) and (7), the theoretical stress-time curve at the centre of 
the plate. This theoretical curve is shown in fig. 4(b) and the agreement 
with the experimental curve is again seen to be excellent. 

In all the above tests the resistance, measured to five figures, of the 
strain gauge at the centre of the plate, was found to be the same after 
test as before test. It was thus assumed that in all tests the plate 
remained elastic. 
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Fig. 4 
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The force—time and stress—-time curves obtained using a 28 in. height of drop 
and a 2 in. diameter weighbar. 


(a) Force—-time curve. 
(b) Stress—time curves. Curves A and B refer to the experimental and theo- 
retical results respectively. 
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4.2. Propagation of a Disturbance Radially 


4.2.1. Impact over a finite area 

It was described in §2 how gauges had been placed in both radial and 
circumferential directions at radii of 1, 2, 3, 4, 5 and 6in. from the centre 
of the plate. In the present tests the 4in. diameter weighbar was impacted 
at a velocity corresponding to a 28in. height of drop and oscillograms 
were obtained giving the output from the central gauge plus that from 
each of the above twelve gauges in turn. The known output of the central 
gauge was subtracted from these dual outputs in order to obtain the 
strain-time curves, both radial and circumferential, at each of the six 
radii considered. In this way the time base for all the strain—time curves 
was the same. Figure 2 (6), Pl. 101, shows a typical strain—time oscillogram 
obtained in the above manner. 

At each radius the radial and circumferential stress-time curves have 
been found from the radial and circumferential strain—-time curves by 
using the formulae relating the principal stresses (p, and p,) to the 
principal strains («, and e«,) at a point on a plate. These formulae are 
given by 

p,r=(e,+c€,)H/(1—o?) oo Ee 
and 
po=(€o+ ce,)H/(1—o?) eo 


where the subscripts r and c refer to the radial and circumferential 
directions respectively. From the derived curves of radial stress—time 
and circumferential stress—time, the radial and circumferential stresses 
have been plotted against radius at times of 20 and 60 microseconds 
respectively. These curves are shown in fig. 5 (a). 

It is to be observed that up to 60 microseconds the disturbance had not 
significantly effected points at any radius greater than about 6in. Hence 
it is considered that, for the time interval over which the theory has been 
compared with the experimental results presented in the previous section, 
the centre of the plate behaved as though the plate was infinitely large, 


4.2.2, Impact at a point 


It was decided to investigate the way in which a disturbance caused by 
an impulsive point load was propagated radially in a plate. A reasonable 
approach to this type of loading arises when a plate is impacted by an 
elastic sphere ; hence tests have been made in which two plates, one dural 
and the other steel, were impacted by a lin. diameter high-tensile steel 
sphere at a velocity of 243in/sec. The dural and steel plates were 
0-133in, and 0-330in. in thickness respectively. In order to prevent 
localized yielding at the point of impact in such tests two methods, one 
for each plate, were used. For the dural plate the sphere impacted a 
small shim of dural, lin. in diameter and J in. thick, placed on the centre 
of the plate; and for the steel plate impact occurred at the centre of a 
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high-tensile steel insert, lin. in diameter and }in. thick, which had been 
pressed into a similar sized cavity machined out of the centre of the upper 
surface of the plate. For each plate the radial and circumferential stress— 
time curves at the centre and at radii of 1, 2, 3, 4, 5 and 6in. were derived 
from the output of strain gauges in a similar manner to that described in 
the previous section. In using eqns. (10) and (11) to derive the curves for 
the dural plate, H was taken as 10 x 10°lb/in? and o as 0:34. For both 
plates the radial and circumferential stresses have again been plotted 
against radius at times of 20 and 60 microseconds respectively; these 
curves are shown in figs. 5(b) and 5(c). A typical strain—time oscillogram 
used in the derivation of these curves is presented in fig. 2(c), Pl. 101. 


§ 5. Discussion AND CONCLUSIONS 


The agreement between Sneddon’s theory and the results of the tests 
using the hammer impact apparatus is exceedingly good. With reference 
to figs. 3(b) and 4(b) the discrepancy, over the first few microseconds, 
between each theoretical and experimental curve can be partially 
explained by the fact that a?/4bt must be small for eqns. (6) and (7) to be 
accurate. For the 2in. and tin. diameter weighbars the magnitudes of 
a?/4bt at a time of 10 microseconds are 0-0512 and 0-0912 respectively ; 
hence it is considered that only for times less than this would eqns. (6) 
and (7) be in significant error. Three other factors must also give rise to 
differences between the theoretical and experimental results. The first 
factor is that the applied force-time curves had to be approximated to 
in order to calculate the theoretical curves; the second is that the effects 
of ‘ rotary inertia ’ and shear forces are neglected in the theory ; and the 
third factor is that the effects of local elastic deformation at the weighbar- 
plate interface have not been taken into account in deriving the theoretical 
stress—time curves. However, the difference between the theoretical and 
the experimental results is considered to be small enough to justify the 
use of the theory in the design of plates to sustain transverse impact loads. 
This could be done provided the impact velocity, the area over which the 
load is applied and the duration of the impact are known. The maximum 
tensile stress in the plate could then be calculated from eqn. (6) if it is 
assumed that the force-time curve applied to the plate is a perfect square 
wave. 

Also it is to be noted that in the tests made, using both the din. and 
gin. diameter weighbars, the plate showed no sign of having yielded. The 
maximum tensile stress reached during these tests was 66 5001b/in? which 
is well above the static upper yield stress (35 0001b/in?) found for the steel 
used in the plate. This increase in strength is to be expected due to the 
fact that during these tests strain-rates of the order of 100sec-! were 
attained in the plate. Impact compression tests (see Campbell and 
Duby 1957) were made on specimens of the same steel as used in the plate 
in order to give an indication of the increase in strength of the steel at 
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such high strain rates. It was found that for a pre-yield strain-rate of 
the order of 100sec™! the upper yield stress was raised to approximately 
85 000 Ib/in?. 

It is of interest to calculate the tensile stresses that would be reached at 
the centre of the bottom surface of the plate if it were loaded statically 
by each of the maximum forces given in fig. 3(a), and to compare these 
stresses with the actual maximum tensile stresses produced in the plate. 
This has been done as follows. Woinowsky-Krieger (1933) has shown 
that for a plate of radius R, clamped round the edge, the tensile stress p 
immediately below the load P is given by 

p=P(1+<)[0-485 loge (R/2h)+0-52]/4h2? . . . . (12) 
where p is independent of the radius a of the area of loading provided 
a<R. Now fig. 3 (a) shows that, for the four impact velocities considered, 
the maximum loads applied to the plate were 4540, 4000, 3160 and 
2240lb. Substituting these values for P in eqn. (12) the corresponding 
tensile stresses are found to be 165000, 145000, 114000 and 81 0001b/in? 
respectively. A comparison of these values with the actual maximum 
tensile stresses produced at the four impact velocities, i.e. 66500, 56000, 
47 500 and 33 500 1b/in?, show that eqn. (12) is quite inadequate for dealing 
with impulsive loads. 

Figures 5 (a)-(c) show that, for the steel plate (0-285in.) impacted over 
a finite area and for the steel (0-330in.) and dural plates impacted at a 
point, the radial and circumferential tensile stresses produced at the 
bottom surface of each plate were still very localized even after 60 micro- 
seconds. Also it is to be observed that in all three tests there was an 
area of compression which propagated radially in advance of the central 
area of tension. These results are similar to those that would be predicted 
from the curves of plate deflection against time derived by Sneddon for 
the particular case of a square loading wave applied to a plate. 
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ABSTRACT 


An examination of the K satellite emission bands observed in lithium and 
beryllium has indicated that the presence of the hole which occurs in the 
core-level during the soft x-ray emission process is unlikely appreciably to 
affect the intensity distribution of the spectra. 


§ 1. INTRODUCTION 


THE method of soft x-ray spectroscopy has been used for a number of 
years as a means of investigating the electron energy distribution of the 
valence bands of solids. Reviews of the subject have been given by 
Skinner (1946) and Tomboulian (1957). Much of the experimental work 
has been carried out on metals and alloys and it is with the soft x-ray 
bands of metals that this paper is concerned. 

Until recently it was accepted that the soft x-ray emission spectra gave 
a general indication of the distribution with energy of the various types of 
electron states in the valence band. That is to say, K spectra gave a 
measure of the density of states of p symmetry, and L spectra the density 
of states of s+d symmetry. However, recent calculations of the band 
structure of lithium (Schiff 1954, Parmenter 1952, Brown and Krumhansl 
1958, Glasser and Callaway 1958) which show a predominance of p states 
at the Fermi surface, are at variance with the observed K spectrum which 
shows a maximum before the high energy end of the band is reached 
(fig. 1). Jones and Schiff (1954) have pointed out that since x-ray 
emission takes place at an atom which has already lost one electron the 
states in the neighbourhood may be strongly perturbed. Skinner et al. 
(1954) do not agree with this view since Schiff’s (1954) calculated band- 
width is considerably different from the experimental value. In the 
present paper an attempt is made to assess the magnitude of the effect of 
the loss of an electron on the band shapes and widths by considering 
these properties in the satellite spectra. 

Soft x-ray emission bands are produced when an electron is ejected 
from an inner shell of an atom and the resulting hole is filled by a transition 
from the band of valence electrons. The intensity is given by the 
expression [(#)=v*. F(Z). N(H) where v is the frequency, F(E) the 
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transition probability, and N (HZ) the density of states. However, in 
addition satellite bands are frequently observed at higher energies, arising 
from double ionization of the inner shells. In this case the radiation 
corresponds to a final state in which there is a hole in the valence band 
and a hole also remaining in an inner shell. If the effect of a positive hole 
is large one would expect a considerable difference between the shapes 
and widths of main bands and satellite bands, and a comparison between 
the two enables an estimate to be made of the effect of the hole on the 
product F(H).N(£). In comparing these bands it is necessary to allow 
for the broadening due to the shorter lifetime of the doubly-ionized state. 


Fig. 1 


Experiment _ Calculated 
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It is the K spectra which are of primary importance for comparison of 
the bands. The unexcited K level has a spectroscopic designation 185. 
Removal of one electron converts it to the state 7S,/., whilst removal of a 
second electron results once more in the state 1S). Owing to the selection 
rules the K spectra represent the p fraction of the valence (or conduction) 
band. The main emission corresponds to transitions down to the 78,). 
state, and the satellite emission to transitions to the 1S, state. Since in 
both types of transition the lower level is a singlet a comparison between 
the two types of emission may be made directly. L satellite spectra are 
much more difficult to interpret. Removal of one electron from the 2s 
state gives the term °S,., whilst removal of one from the 2p state gives 
the terms 2P,, and ?P3.. However, since emission bands corresponding 
to transitions down to the 2S,). level are not observed, due to an Auger 
process which fills the 2s hole, we may disregard ionization in the 2s shell 
and consider only the removal of 2p electrons. Removal of a second 2p 
electron gives the terms *P,, *P;, 3P,, 1D, and 189. The L satellite spectra 
consist therefore of transitions between the p fraction of the valence 
band and the levels 1D, and 1S), and the s+d fraction and the levels 
3P,, 8P, and °P,. Calculations by Kennard and Ramberg (1934) of the 
energy separation between the above five levels show that the two types 
of satellite band must overlap. Consequently interpretation of al 
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satellites is quantitatively very difficult, although an attempt at a 
qualitative interpretation of the bands obtained by Skinner has been 
made here, in order to show that the general features are those that would 


be expected. 


§ 2, EXPERIMENTAL TECHNIQUE 


The spectrograph in use at the National Physical Laboratory and the 
technique used in its operation have been described elsewhere (Catterall 
et al. 1958). 


§ 3. K SaTevxites or Li anp Be 


The K emission main band of lithium has been measured by Skinner 
(1946), by Bedo and Tomboulian (1958) and by Catterall and Trotter 
(1958). The general features of the band are in agreement. In addition 
Tomboulian has measured the satellite band of lithium, with which we 
also agree. The intensity of the satellite is much lower than that of the 
main band, and occurs at 88ev, whereas the main band is at 55ev. In 
fig. 2 are shown the main and satellite bands superimposed to match at 
the peak intensities. It will be seen that the two bands are very similar 
except for an overall broadening of the satellite. We have also measured 
the K satellite band of beryllium which occurs at 134ev, and it is shown 
in fig. 3. 

Fig. 2 
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& Calculated Satellite Points 


The emission intensity distribution of both main and satellite bands 
depends not only upon the electron energy distribution in the valence 
band of the solid but also upon the ‘ natural width’ of the core level. 
Theoretically the line shape of the core level should be Lorentzian, 

= Lo 
a+ (HE) 
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with a width at half-maximum intensity of 2a. Unfortunately the half 
widths of levels for elements with atomic numbers as low as lithium and 
beryllium are not known. However, from a study of the width of the K 
absorption limit in lithium Bedo and Tomboulian (1958) conclude that 
the K level width may be about 0-lev. Using the empirical rule that the 


Fig. 3 
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half-width of a line increases with the square of the excitation energy, the 
half-width of the doubly-ionized state should be about 0-3ev in the case 
of lithium and 0-7 ev in the case of beryllium. It is possible to calculate 
the expected shape of the satellite band from a knowledge of the main 
emission band and the natural width of the doubly-ionized level. The 
calculated shape will only be approximate since the main band will itself 
be broadened by the natural width of the singly-ionized level. However, 
in figs. 2 and 3 are the points calculated by a ‘ folding ’ operation, which 
involves determining the integral 
+0 ’ 1 ; 
L@)= | 10). a 2 


where I,(Z) is the intensity of the satellite band, and [ ,(Z’) the intensity 
of the main band. The agreement between the observed and calculated 
satellite points is quite good. The fit at the low intensity parts of the 
beryllium satellite is only fair, probably because the beryllium satellite is 
more difficult to measure since its intensity is very low. 


§ 4. L SATELLITES OF Na, Mg, Al anp Si 
The bands discussed in this section are those obtained by Skinner (1946) 
and shown in fig. 4. No detailed analysis of the shape of the L satellite 
bands is possible because of the splitting of the L level with the removal 
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of the second electron. The separation in electron volts between the level 
8P, and the others, calculated by Kennard and Ramberg for isolated 
atoms, is shown for the four elements in the table. The 1S level has the 
highest energy. The relative intensity of the emission from the valence 
band (V) to these levels should be given by their quantum weight, and 
should be in the ratio 3P,,) :4D,:1S,=9:5:1. The intensity of the 
satellite bands is usually only 5% or so of that of the main bands, and in 
view of the greater experimental difficulties in the measurements we may 
neglect the level 18, and consider the satellites to be made up from the 
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transitions V->"D, and V-*P,,). The transition V>1D, has the highest 
energy and will be responsible for the high energy edge of the band, whilst 
the transition V-+*P,,9 will overlap at lower energies. The selection rules 
indicate that the former transition will arise from the p fraction of the 
valence band, and the latter from the s+d fraction. The intensity 
distribution given by the p fraction is the same as in the K spectra of 
main emission bands, and the s+d fraction the same as in the L spectra 
Accordingly, a qualitative analysis of the shape of the satellite bends 
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should show that the high energy end has the same features as the main 
K spectra, with a somewhat broader edge, and with a more intense band 
resembling the main L spectra overlapping it at an energy lower by an 
amount equal to the °P,,,-'D, separation. The two overlapping bands 
may not have their intensities in the precise ratio 9:5 since the rule 
relating the relative intensities requires a common initial level before the 
transition, whereas in the present case the initial levels in the valence 
band are different. The most that can be expected is that the satellite 
bands should show a rise in intensity at the lower energies compared with 
the low energy end of the K spectra. 


4.1. Sodiwm 


There are no published data concerning the shape of the sodium 
spectrum, but the theoretical calculations of the valence band show that 
there is a predominance of s states at the Fermi surface (Schiff 1954) so 
that the K spectrum should decrease before the high energy edge is reached. 
Consequently no sharp features such as a peak will be expected at the 
high energy edge of the satellite band (although there should be some 
broadening of the edge) and none are observed. The satellite shows the 
expected broad hump at about 2-5-3-0ev from the edge. The °P,—'D, 
separation in sodium is 3-lev. 


4.2. Magnesium 


The K spectrum reported by Farineau (1937) has a sharp peak at the 
high energy edge (not unlike that in the L spectrum). A peak is also 
present at the satellite band edge, and a broad hump occurs at an energy 
lower by about 4-5ev. The °P,-'D, separation in magnesium is 3:87 ev. 


4.3. Aluminium 


The K spectrum obtained by Farineau (1936) is illustrated by Skinner 
(1946), and like magnesium shows a sharp peak at the high energy edge. 
The satellite also shows a peak, with a very broad hump at about 7ev 
lower. The ®?P,—!D, separation is 4:5ev. 


4.4, Silicon 


The K spectrum of silicon obtained by Farineau (1937) has no sharp 
emission edge but decreases from a maximum to zero over a range of 4ev. 
The satellite band edge shows the same general feature. Skinner reports 
the satellite band as extending to 100-04, and there is a broad hump at 
a lower energy some 10ev from this point. The band V->*Po seems to 
begin at the point 104-14 where Skinner records a change of slope in the 


emission curve, and this is 4-9 ev from the high energy end. The °?P,—'D, 


separation in silicon is 5-lev. we 
The interpretation of the silicon satellite is further complicated by the 


fact that in insulators and semiconductors, unlike metals, an excess 
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positive charge on an ion will definitely produce a bound state lying below 
the valence band (Koster and Slater 1954). In silicon this will be a 3s 
state, and transitions from this level to the lower levels will give a broad 
line on the low frequency side of the spectrum. The energy difference 
between this line and the high energy end of the satellite spectrum will be 
approximately the same as the (3s)?+(3s)(3p) transition in the element 
of next highest atomic number P4*. This is 8-4ev for parallel spin of the 
3p screening electron. The transition accordingly may account for the 
broad hump appearing some 10ev from the high energy end of the 
satellite. 


§ 5. DiscussIon AND CONCLUSIONS 


Examination of the shape of the satellite bands of lithium and beryllium 
shows them to be identical with those of the main band, except for a 
broadening which can be explained by the natural width of the core level. 
In addition, a simple qualitative interpretation of the shapes of the L 
satellite bands seems possible by considering the splitting of the L level 
and the shapes of the K and L main band spectra. 

We may consider separately the possible effect of the extra positive 
hole on (a) the electron energy distribution N(/) and (b) on the transition 
probability F(/) respectively. 

(a) The excess positive charge due to the hole remaining in the inner 
level will be screened either by a level subtracted from the valence band 
in the manner of Friedel (1952) (although Slater (1955) considers this 
unlikely with only one excess positive charge) or by the electrons in the 
valence band itself. In the first case a broad line may be visible, if the 
transition is allowed, on the low frequency side of the band, and in the 
second case the band itself may be somewhat altered in shape. 


(0) The transition probability /(#) for the main emission bands 
depends upon the wave functions of the initial and final states. The 
extension of the wave functions in space will be changed by the increased 
effective nuclear charge due to the presence of the positive hole, and 
consequently we may expect a difference between F(Z) in satellites and 
main bands. a 

Since the extra positive hole, however, has no apparent effect on either 
F(E) or N(E) we may reasonably conclude that the hole in the core level, 
which appears during the main band emission process, does not 
appreciably alter either of these terms. We may also conclude, that the 
suggestion of Jones and Schiff, that the presence of the positive hole could 
account for the large discrepancy between the observed and calculated 
band forms of lithium, is extremely improbable. 
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ABSTRACT 


High explosives have been employed to generate intense shock waves in 
Perspex and from the measurements made of the shock and free surface 
velocities the ‘ Hugoniot ’ equation of state has been determined within the 
pressure range from 3 to 13x10! dynes em™. The results fit an equation 
of an inverted Bridgman type. Graphs are given of the variation of pressure, 
shock and particle velocity with Perspex thickness. 


§ 1. INTRODUCTION 


Tue detonation of high explosives in contact with a solid results in the 
transmission of a steep fronted pressure pulse into the material. The 
pressures obtained by this means are greater than the highest pressures 
reached in static experiments such as carried out by Bridgman (1945 a, b, ¢, 
1949). The propagation of these high intensity shock waves in metals is 
a subject which has received considerable theoretical attention in the 
last few years (Duval and Zwolinski 1955, Walsh and Christian 1955, 
Goranson et al. 1955) and it has been shown that an equation of state 
can be derived from simultaneous measurements of the velocity of the 
shock wave and the particle velocity (mass velocity) in the wave. The 
pressure volume relations (often called Hugoniots) obtained from 
experiments of this nature are in a pressure region so far unattainable 
by other methods. These dynamic experiments also have the advantage 
that they avoid the use of cumbersome apparatus, with its inherent 
difficulties of unknown relative distortions which lead to large 
uncertainties in measurements. The results from the two types of 
experiments are not immediately comparable since in one case the 
pressure is applied very slowly and under isothermal conditions, while 
in the other the pressure is applied extremely rapidly (< 10-6 sec) with 
an associated rise in temperature. Nevertheless the difference between 
the two cases is not large and it is possible to convert from one to the 
other by the use of thermodynamic relations (Goranson et al. 1955). 

Much experimental work has been carried out along these lines since 
the pioneer experiments of Pack et al. (1948). Their original technique 
has been developed (Minshal! 1955, Bancroft e¢ al. 1956), optical 
techniques have been used (Walsh and Christian 1955, Walsh and Rice 
1957) and as a result a large amount of data is available on the 
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compressibility of metals (Walsh and Christian 1955, Minshall 1955, 
Bancroft et al. 1956, Walsh ef al. 1957). Some work has been carried 
out on liquids (Walsh and Rice 1957, Rice and Walsh 1957) and air 
(Deal 1957), but very little has been reported about the transmission 
of these shock waves through non-metallic solids (Lawton and Skidmore 
1956). Due to the increasing utilization of plastics it was considered 
desirable to obtain quantitative information about their behaviour. 
Perspex was chosen since it is widely used, readily available and easily 
machined, which facilitated the preparation of the test assemblies, 


§ 2. THEORY 


The equations based on the conservation of mass and momentum 
across a plane shock wave were derived by Rankine and Hugoniot and 
are respectively, 


and 
dS gic NBD haa Pd ee, Sete a. (G2) 


where P, p and wu are the pressure, density and particle velocity in the 
shock wave which is travelling with velocity D relative to the material 
at rest, with initial density py at pressure Py. PP, is usually taken as 
zero when P> P,. 

The shock velocity, D, can be measured directly but the particle 
velocity, wu, in a shock wave must be deduced from the measurement 
of the free surface velocity which results from the normal reflection 
of the shock wave at the free surface. There must be zero pressure 
difference across the free surface boundary, and this condition requires 
that a wave equal in amplitude but opposite in sign must be reflected 
from the surface. Under these conditions the surface will acquire a 


velocity 
v=ut+o 7 Cube Wea, are eae) 


where u is the particle velocity in the stress wave, and a is the excess 
of the free surface velocity over the particle velocity uv. If it is assumed 
that the unloading takes place isentropically then o is the Riemann 
function and is given by 


P 

P dv\ 12 
x BS ade ptremmumred i ali wrt LEY, 
li a 


where subscript s denotes constant entropy, and c is the velocity of 


sound. 

The approximate relation 
i= CL ee ee) 
when combined with eqns. (1), (2) and (3) enables a (P, p) point to be 
evaluated for each set of pp, u, D data, 


D 
P.M, ? 
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A more exact relation between wu and o has been given by Goranson 
et al. (1955, eqn. (10)), but the thermodynamic data necessary to evaluate 
this equation are not available for Perspex. However, Walsh and 
Christian (1955) have given relations setting upper and lower limits to 
the ratio o/u. These are 


(lar) log. Cole erate ort! 


where B=C,,p)/z, and C,, is the specific heat at constant pressure and 
x the thermal coefficient of volume expansion of Perspex, and 


Gb ALC a) 1S oe eae 


w [P(% —v)]*? 


where the subscript H denotes integration along the pressure-volume 
curve. An approximate pressure-volume curve can be derived using 
eqn. (5) and used to evaluate eqn. (7). This can then be used to construct 
a corrected pressure-volume curve and the process repeated until self- 
consistency is obtained. This method was not employed since any error 
introduced by the approximation w=o is less than the inherent 
experimental errors (see § 5.3). 


§ 3. EXPERIMENTAL METHOD 
3.1. Aim 


The experiments were designed to measure the velocity of propagation 
of a shock wave passing through Perspex, and also the velocity of the free 
surface caused by the normal reflection of the shock wave at the 
Perspex—air boundary. The shock wave was produced by detonating 
an explosive charge in contact with Perspex, the maximum pressure 
generated in the Perspex being of the order of 101! dynes/em?. The 
Perspex was in the form of blocks of various thicknesses. 


3.2. Hxplosive 

The explosive charge consisted of a mild steel tube measuring 3 in. 
long by 3in. internal diameter, with in. walls, filled with plastic 
explosive (P.E. No. 2) (see fig. 1 (a)). Initiation was achieved by means 
of a detonator, electric, No. 8 Briska, placed in contact with a tetryl 
pellet, } in. dia. by } in. long, which was countersunk flush into the top 
surface of a detonation wave-shaper. The latter generated a plane 
detonation wave in the plastic explosive charge. 


3.3. Target Plates 
The test plates were 4in. x 4in. and cut from Perspex sheet up to thick- 
nesses of lin. Greater thicknesses were built up by carefully placing 
sheets together and taping with Sellotape, The two large surfaces were 
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machined to be parallel and were flat to within 0-0015 in. over a diagonal. 
The average roughness was about 15 micro-inches. 


3.4. Probes 


Three types of double probe were used in the experiments to determine 
shock and free surface velocities. The probes acted as high speed switches 
in the event circuit shown in fig. 1 (c) and consisted essentially of two 
poles, at a potential difference of 300 volts, separated by about 0-002 in. 
The poles were joined electrically when a shock front produced intense 
ionization across them, thus discharging the event circuit condenser into 
a matched coaxial cable which transmitted a pulse directly to the 
deflector plates of a high speed oscilloscope. 

A probe made from two nylon-covered copper wires (32 s.w.g.), glued 
together as a flat twin arrangement, was placed in the P.E. surface where 
it met the Perspex surface. The arrival of the detonation front at the 
probe caused the two wires to short together and thus operate the event 
circuit. The twin-wire probe used for the shock velocity measurements 
was made from two nylon-covered copper wires (32 s.w.g.) twisted 
together. The probe end of the wires was cut level and carefully rubbed 
flat on abrasive paper. The insulation of the probe was then checked 
at 500 volts. 

For the determination of the free surface velocity, a probe of more 
rigid construction than the twisted wire type, but still of small end-area, 
was required. A double probe was developed of concentric form which 
consisted of a copper tube at earth potential, 0-02 in. bore and 0-032 in. 
outside diameter, containing a nylon-covered wire glued in the bore so 
that the base of the wire was inset 0-002 in. with respect to the end of 
the tube : the end of the tube and wire had been carefully rubbed flat. 
The probe was soldered into a hole drilled along the axis of a 6 B.a. 
brass screw. This enabled the probe to be easily set in its carrier to the 
required distance above the free surface. A small piece of Sellotape 
covered the end of the probe thus enclosing a bubble of air which would 
ionize readily when compressed. This was necessary as, just before 
firing the charge, the probe was enclosed in an atmosphere of propane 
which does not ionize easily. 

An account of experiments with these and other forms of probe is 


given in § 5.2. 


3.5. Shock Wave Velocity 
The shock wave velocity was determined by measuring the time of 
arrival of the shock front at measured distances apart (fig. 1 (a)). For 
these experiments, the explosive was placed vertically above the test 
plate. Eight twisted wire probes were placed into flat-bottomed holes 
drilled on a circle of 0:6 in. dia. at the centre of a | in. thick Perspex 
plate. The holes differed in length by }in. and were grouped so that 
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Fig. 1 
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distortion of the shock wave produced when the shock reached a probe 
did: not affect those parts of the wave which would encounter the following 
probes. In the first series, the probes were in contact with the bases 
of the holes, and in the second they were set up 0-002 in. from the bottom 
of the holes. Each probe was connected into the standard event system 
(fig. 1 (c)), and functioned as described in § 3.4 above. 


3.6. Free Surface Velocity 


The free surface velocity was determined by placing concentric probes 
at known distances above the surface (0-005 in. to O-110 in.) fig. 1 (5), 
and obtaining from oscilloscope records the time of arrival of the surface 
at these probes. Eight of the probe screw assemblies described in § 3.4 
were fitted into a brass plate with tripod legs, and this rested on the test 
surface. The probes were equi-spaced along the circumference of a 
circle of 0-6in. diam. at the centre of the Perspex. Since the probe 
carrier rested on the test surface under its own weight, it was necessary 
to place the charge vertically below the target plate. The probe carrier 
served as a common earth and the centre wire of each probe was connected 
to an expendable event system mounted above the probes. 


3.7. Recording 


Three high velocity linear oscilloscopes and one spiral oscilloscope 
were used to record the time between event pulses. Initially Mark IIB 
linear oscilloscopes were used, which had been designed in the A.R.E. 
Instrument Section, but subsequently improved versions, Type 58, 
developed by Mr. Barnsley, were used. The Type 58 oscilloscope contains 
a new cathode-ray tube, 4EP1, which gives smaller spot size, greater 
deflection sensitivity and improved frequency response in comparison 
with the cathode-ray tube used in the Mark IIB. No published figures 
for the performance of the two types of oscilloscope are available. The 
following are approximate values under the conditions of the experiment, 
Type 58 oscilloscope quoted first : deflection sensitivity, 60 volts/em : 
120-140 volts/em: frequency response, 400 Mc/s; 200 Mc/s: spot size, 
0-55mm;0-9mm. The velocity of the spot across the cathode-ray tube 
may be varied so that the transit time will be 0-1 or 0-2 or 0-5 increasing 
by a factor of 2” to 128 psec. The Type 58 oscilloscope functions 
adequately on all these ranges whereas the Mark IIB oscilloscope is not 
satisfactory on sweep times of less than 2 psec duration. For the 
experiments on Perspex, the time ranges used were of 0-5, | or 2 usec 
duration, depending on the total time that event pulses were expected 
to occupy. The spiral oscilloscope was of the blanking type with a 
collapsing spiral, one revolution of the spot being crystal-controlled to 
take 2 psec. The event pulses were fed into the cathode of the cathode-ray 
tube, positive pulses causing gaps in the trace and negative ones 
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brightening the trace. A one shot camera containing a Wray f1-0 lens, 
giving a reduction in image-size of 4: 1, was fitted to each oscilloscope 
to photograph the trace. The time resolution attained was +5 musec 
on a microsecond trace in the case of the Type 58 oscilloscopes, and 
+ 10 mpsec for the spiral oscilloscope. 

A block diagram, fig. 2, outlines the electronic equipment used and the 
method of connection. Coaxial cable, of 100 ohms impedance, was used 
for all connections in the pulse circuits. A twisted wire probe inserted 
into the charge triggered the D.P.S. unit (Delay Pulse Synchronizer) 
about 3 psec before the detonation wave would have reached the Perspex 
surface (fig. 1 (a)). The function of the D.P.S. was to emit positive 
pulses, either without delay (Channel A) or after selected delays 
(Channels B, C and D) with respect to the trigger pulse. The spiral 
oscilloscope was tripped from Channel A of the D.P.S. and recorded the 
complete sequence of event pulses, all of which occurred within 30 psec. 


Fig. 2 
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The linear oscilloscopes were triggered from Channel B after a set delay 
which ensured that all the event pulses with the exception of that due 
to the surface probe would be displayed. The linear oscilloscopes were 
connected in cascade, the start of the traces on the second and third 
being delayed with respect to the first and second, respectively, using 
the beam delay control on the oscilloscope. A crystal-controlled time- 
marker unit fed pulses on to one Y-plate of each oscilloscope for the 
purpose of calibrating the traces : the other Y-plate received the event 
pulses. The three traces were linked by displaying two marker pulses, 
each one being common to consecutive pairs of traces with the minimum 
of trace overlap. These pulses were obtained from Channels C and D 
of the D.P.S. via the time-marker unit. Thus the time resolution was 
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increased in comparison with the simpler arrangement. employing one 


oscilloscope. The spiral oscilloscope, in addition to recording the surface 
probe pulse, enabled a satisfactory result to be obtained if the linear 
oscilloscopes failed, due to faulty delay setting, to record the event 
pulses. 

The oscilloscope traces were photographed on Ilford H.P.3 and 
developed in I.D.33. Typical records are shown in fig. 3, PIP102: 


§ 4. ExpERIMENTAL RESULTS 
The results of experiments on shock velocity measurement are shown 
in graphical form, fig. 4. The common zero for each experiment, 
irrespective of block thickness, was the tripping of the surface probe. 
This was assumed to trigger immediately the detonation wave reached it, 
so that a correction was applied for the time taken for the shock to reach 
the Perspex. No difference was observed between the results obtained 


Fig. 4 


K 


Z 2 ee 1 
- = 7 > © PROBES IN CONTACT WITH PERSPEX 


x PROBES SET UP 0.002" FROM PERSPEX 


SI DISTANCE (MM) 
” 


| 
| 
| | 
| 


| | 
5 75 10 
TIME (4 SECONDS) 


Distance versus time plot of shock front in Perspex. 


from probes in contact with Perspex and those set up 0-002in. In the 
latter case, a correction—varying from 0-01 to 0-03 psec—was made to 
the times recorded. This correction represented the time taken by the 
air shock to travel from the measured end of the hole to the probe. 

A quadratic equation, calculated by the method of least squares, was 
found to fit the exprimental points : 


= 0183s + 0-0005s? Se I oe ok) 


where ¢ is in psec, s in millimetres. The velocity at a point s within 
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the range covered by the experiments (s<75 mm) was obtained by 
differentiating eqn. (8) : 
1 


ee (9) 
0-183 + 0-001s 


shock velocity, D 
where D is in mm/psec. The equation gives a maximum shock velocity 


of 5-47 mm/sec. Figure 5, curve 1, shows the variation of shock wave 
velocity with Perspex thickness over the range covered by the experiments. 


Fig. 5 
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Curve |: Shock wave velocity—calculated from fig. 4. 
Curve 2: Free surface velocity—calculated from table 1. 
Shock wave velocity and free surface velocity versus Perspex thickness. 


Table 1 gives the results for the free surface velocity experiments on 
thicknesses ranging from } in. to 3in. The average velocity results were 
plotted for each thickness, and a best fitting straight line calculated, 
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This is shown in fig. 5, curve 2. The equation of the line is : 


free surface velocity, V =4:03 —0-035s. ess gece 
Using eqns. (1), (2), (3), (5), (9) and (10), the dynamic compression 
(Hugoniot) curve, fig. 6, was constructed. A value of 1-2 g/em?® was 
used for the density of Perspex. 


Table 1. Perspex : Results of Free Surface Velocity Experiments 
I ee eee 


Perspex thickness Free surface Standard 
Oe Standard pi 
————_—_—_—_—_————| __ velocity Gheenen Mean deviation 
(Nominal) (Actual) (mm/sec) a. of mean 
in. mm 
3°18 3-92 + 0-09 
$ 3:27 4-06 + 0-16 4-33 — 
3°20 5-03 + 0-28 
5:97 3°85 + 0-25 
6-17 3°28 +011 
4 6-22 3°95 +0-20 3°63 + 0-13 
6-22 3°67 + 0-29 
6:22 3°38 +0:12 
12-62 3°31 +0-18 
12-39 3°34 +0-06 as. 
4 12-22 3°55 + 0-06 3°28 +0:10 
; 12-55 2-92 +0-11 
12-34 3-28 + 0-07 
24.93 2-98 + 0-06 
1 24-93 3°05 + 0-06 3°19 = 
24-77 3-53 + 0-08 
37-5 3.08 + 0-05 2 eas 
15 37-4 2-66 #013 ef 
49-6 2-24 + 0-06 2.294 == 
2 49-7 2-23 + 0-06 
61-6 1-98 +0-i3 1-89 2 
25 63-1 1-80 +0-02 
75-0 1:39 + 0-02 1-40 = 
3 73.8 1-41 +0-06 


After detonating the explosive charge in contact with the Perspex 
plates, a few small fragments of Perspex were recovered in the case of 
the 3 7 thick target plates only. An odour, characteristic of vaporized 


Perspex, could be detected in the firing chamber. 
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§ 5. ERRORS 
The errors inherent in these experiments are, in the main, random 
and are discussed below together with a possible more serious systematic 
error due to the gaseous shock which precedes, and is generated by, the 
moving free surface. 


Fig. 6 
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Calculated from eqns. (1), (2), (3), (5), (9) and (10). 
Dynamic compression curve for Perspex. 


5.1. Random Errors 

It is implicitly assumed in the measurements of both shock and free 
surface velocity that the shock wave refracted into the Perspex block 
is plane and parallel to the surfaces of the block. Any deviation from 
this condition would have resulted in some of the probes being struck 
earlier or later than should have been the case. Since the probes were 
arranged to increase in distance in both directions from the first probe, 
the usual effect would have been an increased scatter in the times recorded. 
A study of the results obtained indicates that the error from the skewness 
of, and irregularities in, the shock wave was not great. 

The shock velocity results for blocks of different thicknesses were 
combined by taking the probe event at the explosive Perspex interface 
as the time zero. Due to variations in the surface flatness of the Perspex 
block, this probe could have varied in its distance from the other probes 
by +0-0015in. (+0-04mm). This is equivalent to a time error of 
+ 0-007 usec at the maximum velocity and since the times were measured 
to +0-01 psec the resultant error in the time zero should be less than 
0-02 psec. The twisted wire probes inside the Perspex block were set in 
holes whose depths were measured to an accuracy of +0-0002 in. 
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(+0-005 mm). The position of these probes might have varied due to 
the way the probe was pushed into the hole and glued in position, or 
handling and temperature effects might have caused movements in the 
probes after setting. The error in the shock velocity determination was 
estimated from the combined distance—time results (fig. 4) as less than 
+ 2% and the major source of error was in the combining of the results 
from different experiments. 

In the case of the free surface velocity the probes were set in the probe 
carrier to the required distance on a flat table fitted with a dial gauge. 
This enabled probe distances to be measured to +0-0001 in. The probe 
carrier was then transferred to the Perspex surface which was not as 
flat as the measuring table. The error in the distance from this operation 
was minimized by arranging that the probes were set on a circle at the 
centre of the test plate. The design was such that temperature effects 
were self-cancelling but handling might have caused some variation. 
Another possible source of error was in the probe. The centre wire 
might not have been exactly 0-002 in. inset and this would give an error 
in the time of an event. The centre wire is unlikely to have been 
misplaced by more than +0-002in. and the corresponding variation 
in time is +0-01 usec. All the above errors would lead to increased 
scatter of the experimental points about the ‘ true’ line. 

The experimental values obtained for the free surface velocity are 
given in table 1. The accuracy is generally better than +5°% for each 
experiment though the variation between nominally identical rounds 
is much larger in some cases. (This is discussed later.) The scatter of 
the experimental points is on the average omy slightly larger than that 
expected from the accuracy of reading the oscilloscope traces ( + 0-005 psec 
corresponding to +0-0006 in. in distance) indicating that the effect of 
all the other sources of error is about the same size as the time resolution. 


5.2. Systematic Errors 


The movement of the free surface sets up a shock wave in the gas ahead 
of the surface. This shock wave moves with a greater velocity than the 
surface and, if the probe system responds to the shock, an enhanced 
velocity will be recorded. The relation between the free surface velocity 
and the shock velocity in air and in propane has been calculated (Deas, 


private communication). ; “a 2 

Various experiments with probes of differing anticipated sensitivity to 
shock were carried out, and the results are given in table 2. Four types 
of probes were used ; a brass screw making contact with a conducting 
surface, a concentric probe with Sellotape over the end, a concentric 
probe without Sellotape and a twisted wire probe. _The results using 
tin. thick aluminium plate indicate the relative merits of the systems. 
The brass screw array in propane gives the true free surface velocity, 
while the concentric probes with Sellotape in propane or in air give mean 


results which are not significantly different (one result is high and would 
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agree with the propane velocity). However, the concentric probes without 
Sellotape in air give a significantly higher result which corresponds to the 
expected air shock velocity. 

The results using the 2in. thick Perspex show that the free surface 
velocity is approximately the same as that of the 4 in. thick aluminium 
sheet, and thus we would expect the concentric probe with Sellotape to 
give the correct velocity. However, there is a tendency to get a higher 
result when the probes are surrounded by air rather than propane. 


Table 2. Measurement of Free Surface Velocity. Experimental Results 
from Various Probe Systems in Air and Propane 


locit 
Test material Gas Probe system Peet Bi 
Propane Brass screw 2:34+0-06 
do. Concentric, Sellotape £ iis 20-8 
eat 2-56+0-07 
eee Brass screw . 2-20+0-10 | ft 
din. thick do. ; 2 ne 
oe Concentric, Sellotape 2-28+0-07 f 
Air Concentric, Sellotape 2-36+0-07 
do. Concentric, no Sellotape 2-73 +0-07 
(2-244 0-06 
Propane Concentric, Sellotape 2-23 + 0-06 
Pp 1-95+0-14 
Deedee tee 
ees Air Concentric, Sellotape 2:23 + 0-06 
2-47+0-10 
do. Concentric, no Sellotape 2-47+0-05 
3°85 + 0-25 
3:28+0-11 
Propane Concentric, Sellotape 3-95 + 0-20 
Perspex 3-67 + 0:29 
tin. thick 3°38 + 0-12 
do. Twisted twin wire, no Sellotape| 3-83 + 0-08 
Air Concentric, Sellotape 4-10+0:10 
do. Concentric, no Sellotape 4-07 + 0-03 
s No events. 
erspex en aiek ees etna 4 : Insufficient 
ee Rae Propane Twisted twin wire, no Sellotape EERE ci 
trigger probes 


+ Probes mixed, four of each type, same target plate. 


The third group of firings using }in. thick Perspex gave a large 
variation in the results from experiments of the same type. One 
interpretation is that in some cases we were recording the movement 
of the shock front and in others the true free surface velocity. The 
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value obtained for the air shock velocity using the concentric probes 
without Sellotape supports this interpretation, as does the high velocity 
figure obtained in propane using an uncovered twisted wire double probe. 
However, the experimental errors in these firings are of the same 
magnitude as the difference between propane shock and surface velocity 
and a definite decision must await improvement of the probe system. 
Any error introduced by this will only affect the high pressure end of 
the pressure-volume curve and will be no larger than that due to the 
scatter in results. To keep this uncertainty as small as possible the 
probes were immersed in propane rather than air, since there is a larger 
difference between the free surface velocity and the resulting air shock 
velocity than between the free surface velocity and the resulting shock 
velocity in propane. 
5.3. Observations 


In order to obtain an insight into the effect on the calculated Hugoniot 
of the experimental errors in the shock velocity and free surface velocity 
results, the following calculations were made. A nominal thickness of 
Perspex was considered and the shock velocity calculated from eqn. (9) 
increased or decreased by 2%. These values were then combined with 
the maximum and minimum values for the free surface velocity at the 
selected thickness (from table 1) im such a way as to give the greatest 
range in values for the compression ratio. The results obtained at 
various thicknesses are given in table 3 together with the values for 
the pressure at the corresponding thicknesses calculated from the shock 
velocity (eqn. (9)) and the mean free surface velocity. 


Table 3 


Pressure Compression ratio range 
(dynes/em) (p/p0) 


4-5 x 1010 1:28 to 1-33 
9-0 107° 1-44 to 1-60 
12-010" 1:46 to 
12-6 x 101° 1-56 to 


The particle velocity was taken as half the free surface velocity. The 
possible maximum and minimum values for the particle velocity at 
various pressures were calculated using eqns. (6) and (7). Te apo 
heat of Perspex at constant pressure was found to be 0:33 cal g—! deg c+. 
The compression ratio obtained from these values varied from 1:47 to 
1:52 at a pressure of 101! dynes per cm?, the discrepancy decreasing 
rapidly as the pressure reduced. The error introduced by taking 
u=4 free surface velocity may thus be neglected in comparison with 


the experimental errors, 
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§ 6. DiscussIon 


These experiments have enabled compression data for Perspex to be 
obtained in a region not covered by other experiments. The results 
can be fitted by an expression of the form 


p=A (#0) +B (pore) as eRe 

Po Po 
where A =6-64 x 10! dynes/cm? and B= 26-2 x 10!° dynes/cm?. This fits 
the results well over the range covered by the experiment but extra- 
polation to greater pressures must be used with great care. The only 
other published information in this region is the single point which 
can be obtained from the results of Lawton and Skidmore (1956). Their 
point (P=1-77 x 104 dynes/cm?, p/p) = 1-63) agrees well with the values 
of p/p» =1:71 given by eqn. (11) for the same pressure. Pressures up to 
~ 10x 10!° dynes/em? have been reached statically but no results have 
been published for Perspex under these conditions. It would be very 
interesting to be able to compare the values of A and B obtained in a 
slow isothermal experiment with the values obtained in the present 
experiment. 

It is possible to compute the velocity of propagation of elastic 
dilatational waves from eqn. (11) (i.e. as P+0) but this involves an 
extrapolation of the experimental results and no great reliance can be 
placed on the value obtained. The velocity thus obtained is 2-4 mm/ysec. 
Kolsky (1956) has published results which enable a velocity of 2-3 mm/sec 
to be deduced for the velocity of longitudinal waves along rods of 
Perspex. Kolsky, private communication, suggests that this value of 
the longitudinal velocity corresponds to a dilatational wave velocity 
of 2-5-2-6 mm/sec, which agrees reasonably well with the prediction of 
eqn. (11). Equation (9) is an empirical relation between shock velocity 
and thickness of Perspex and is valid only over the range covered by 
experiment. It is obviously incorrect at large thicknesses of Perspex 
as it indicates a continuing decrease in shock velocity with increasing 
Perspex thickness instead of an asymptotic approach to the sonic velocity 
as given by eqn. (11). 

A feature of the experimental results for the velocity of the free surface 
(tables 1 and 2) is that there was a greater round to round variation than 
might be expected from the standard deviation of each experiment. 
The source of this round to round variation is not clear ; it could be 
due to variability in the explosive charge, in the refraction of the shock 
wave into the Perspex or in the Perspex or a combination of all three. 
The effect of the first two could be eliminated if it had been possible to 
measure both shock velocity and free surface velocity in the one 
experiment, but this was not possible due to the limitation on the size 
of explosive charge which could be used. The other possible source of 
the round to round variation is that in some of the free-surface velocity 
experiments some of the probes were responding to the shock in the 
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Fig. 7 
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Pressure versus Perspex thickness. 
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propane gas as mentioned in § 5.2. This interpretation is supported by 
the results which show that there is a smaller round to round variation 
in the shock velocity experiments where this effect is absent. To obtain 
greater accuracy in the experimental results it will be necessary to improve 
the probe system used to detect the movement of the free surface, 
particularly at the higher velocities. It is believed that with some 
development this difficulty should be overcome and more consistent 
results obtained. 

In addition to the graph of P vs p/p, the results have been plotted in 
the form of pressure against thickness of Perspex (fig. 7) which gives the 
effective attenuation of a plane wave in Perspex. This is useful in the 
design of systems where high amplitude shock loading occurs. The 
slope of the graph of P vs u (fig. 8) is the shock impedance (p D) of Perspex 
which is seen to increase steadily with increasing pressure. This means 
that the shape of a high amplitude stress pulse will be changed when it 
is refracted at a Perspex boundary since the amplitude of the refracted 
wave depends on pD. Figure 8 will also be useful in calculations on 
high speed impact where pressures of this order are met. 
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ABSTRACT 


The preparation and basic dielectric properties of the metaniobates and 
metatantalates of iron, manganese, lead, tin, calcium, strontium and barium 
(all in the divalent state), and trivalent aluminium are described. Manganese 
metatantalate, barium metaniobate, and aluminium metaniobate have been 
found to be ferroelectric, and lead metatantalate antiferroelectric. These 
are in addition to ferroelectric lead metaniobate discovered by Goodman 
in 1953. Strontium and barium metatantalates may have transition points 
below —183°c. 


§ 1. InTRopDvCTION 


Surveys of the properties of perovskite type materials of general formula 
A” B®” O,-, together with a review of the whole field of ferroelectrics, can 
be found in works by Shirane, et al. (1955) and Kanzig (1957). In 1953 
Goodman prepared lead metaniobate, PbNb.O,, and found it to be ferro- 
electric. Ifthe general perovskite formula for n = 1 is written as A4*+(BO,)!-, 
then the formula which includes lead metaniobate is A?*(BO,)!-,, which can 
be regarded as the second example in a sequence where the valency of the 
A ion increases by unity each time, and B remains an ion from group 5. 

Metaniobates and tantalates of divalent metals occur naturally as the 
two dimorphous mineral groups columbite and tapiolite, which are each 
finely graded series having the general composition iron-manganese 
metaniobate-metatantalate. Many specimens of these minerals were 
obtained from the Bristol Museum, British Museum, and Universitat 
Museum of Oslo, to whom the authors are indebted for their cooperation. 
However, despite good or excellent habit the crystals were mesomorphic, 
and the axes of the constituent crystallites were poorly defined ; some had 
one axis common to all crystallites, while others had no common axis at 
all. Moreover, the electrical conductivity was so high as to preclude any 
reliable electrical information being obtained so that all specimens had to be 
produced in the laboratory. 

The research to be described surveys the dielectric properties of the 
metaniobates and metatantalates of iron, magnanese, lead, tin, calcium, 
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strontium and barium, all in the divalent state. In order to make a pre- 
liminary investigation into the third formula type, A**(BOs;)'3, the 
metaniobate and tantalate of aluminium were also prepared. 

As ceramic specimens only were investigated, and no single crystal elec- 
trical measurements made, the results are not presented as a complete 
description of the electrical properties of these materials, but rather for the 
guidance of those concerned in dielectrics research. 


§ 2. SPECIMEN PREPARATION 


The starting materials were the pentoxide of niobium or tantalum and the 
sulphate, carbonate, or oxide of the A metal. Whenever possible com- 
pounds of the A metal were used which decompose at a high temperature 
to give the required oxide, for these are far more reactive (e.g. Cohn 
1948, Goodman 1953). In most cases alternative starting materials 
were tried for the same final product to ensure that the desired form was 
obtained. 

Approximately stoichiometric amounts of the two starting materials 
were ground together under acetone, and, when almost dry, pressed into 
discs of 0-5in. diameter by 0-2in. thick. The prefiring temperatures were 
about 50°C above the decomposition temperature of the sulphate or car- 
bonate used, and this was maintained for 2 hours, achieving the correct 
value from 400°¢ in 1} hr. 

After regrinding under acetone and drying off to a putty-like consistency 
the final discs were pressed at 10 tons/sq. in. ina stainless steel die, and fired 
for 24 hr. at about 50°c below the fusion point, provided this did not 
exceed 1500°c. All firing was performed with an oxygen atmosphere 
which, together with a 5% excess of the B,O; constituent in the original 
mass, gave specimens of lowest conductivity. X-ray powder diffraction 
photographs were taken at each stage to check and follow production of the 
required material. 

The disc faces were ground plane parallel, the dise porosity measured 
by water absorption, and finally aluminium electrodes were evaporated on. 
to the faces. 


§ 3. Meruops or MEASUREMENT 


For all electrical measurements the disc was held in a combined specimen. 
holder and heating muffle. Below room temperature measurements were 
made at two fixed points—the temperature of boiling liquid oxygen, and 
that of vaporising carbon dioxide. Here, the disc was first maintained 
at 150°c for half an hour to ensure its being thoroughly dry, and then 
plunged rapidly beneath the coolant. Sufficient time was allowed for the 
assembly to reach the required temperature before any measurements were 
taken. 

Dielectric constants were measured in two ways. First, a modified 
Sawyer Tower (1930) display apparatus was used to give a visual indication 
of the displacement, (D), to field, (E), relation at 50 c/s. An indication of 
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the presence of ferroelectricity is given by saturation of the displacement at: 
the tips ofthe loop. The peak to peak dielectric constant at the given applied 

field can be calculated from the geometry of the loop and from a knowledge 
of the dise electrode area and disc thickness. The variation of dielectric 
constant with temperature for a fixed field value, and the variation of 
dielectric constant with applied field at constant temperature were measured 
for each material investigated. 

The dielectric constant values so obtained are the moduli of the complex 
dielectric constants, and therefore each includes a contribution from the 
loss component. When the conductivity of the specimen becomes too 
large, generally for temperatures above about 200-250°o, the real com- 
ponent of the modulus is swamped, and an inaccurate value of 
dielectric constant results. 

The second method of dielectric constant measurement was by means of a 
@-meter. This gives the real part of the complex quantity only, and there- 
fore in general will be smaller than the value obtained on the display 
apparatus. Moreover, since fields of only a few volts/cm instead of many 
kv/cm are used, the incremental value is obtained which is again usually 
smaller than the peak to peak value for non-linear dielectric materials. 
All discs were measured for dielectric constant and loss variation with 
temperature up to 500°C, at frequencies of 0-1, 1 and 5 Me/s. 

All quoted values of dielectric constant have been corrected for specimen 
porosity by the formula due to Van Santen (1946): 


ke, keq = 1/(1— 3p) 
where 4k, is the real dielectric constant for no voids, 


k,, is the measured dielectric constant, 


p is the specimen porosity. 


§ 4. SuRvEY or RESULTS 


The main results for each material are given in the accompanying table. 
Additional details are given below for the most interesting cases, and full 
details of tan § vs 7’, 1/kvs T and kvs E for each specimen tabulated may be 


found elsewhere (Coates 1956). 


4.1. Lead Metaniobate, PobNb,O, 


The ceramics prepared by the authors gave clear saturation trends 
on their D/E loops, with 50 c/s dielectric constant value of 26 at —183 c 
rising to 90 at 200°C. No transition temperature was observed, as this. 
is at 570°c which is higher than the range examined. The room 
75, is half that quoted by Goodman, which can 
existence of a non-ferroelectric rhombohedral 
polymorph (Francombe 1956) produced for firing temperatures below 
1200°c. Since the disc investigated here was fired first below 1200°c 
and second at 1230°c it is a mixture of both forms, as is confirmed by 


temperature value of &, 
be explained from the 
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its x-ray powder diffraction pattern. The dielectric constant value is 
thus reduced by the amount of rhombohedral phase present. 


4.2. Lead Metatantalate, PbTa,O, 


Figure 1 shows the k vs 7’ graph for this material. It is immediately 
obvious that some important transition occurs at about 70°c ; the 
curve of dielectric constant against temperature falls steeply above and 


Fig. 1 


©:50 c/s; E= 31 kV/cm. 
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Relation between dielectric constant and temperature 
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below this point for all frequencies. Despite this peak, the D/E trace 
at 50 c/s was a thin ellipse, tending to a straight line *: — 183°o, with 
no sign of any saturation effects. This suggests, therefore, that lead 
metatantalate is not ferroelectric, but might be anti-ferroelectric, Tf 
this is so, then at sufficiently high fields the trace should split into two 
small hysteresis loops, joined together by a linear portion where the field 
is insufficiently strong to cause a parallel arrangement of the dipoles. 
Cross (1954), working on sodium niobate, needed a field of 90 kv/em to 
achieve this splitting, whereas the highest field applied in this instance 
was only 43-6kv/em. Similarly, Shirane and Pepinsky (1953) were 
unable to split the trace of lead hafnate with a field of 40 kv/cm, despite 
a dielectric constant peak of 540 at 215°c. 

A plot of 1/k against 7 for the three sets of rf. readings shows that 
above 70°c the points lie on three straight lines, confirming that this 
region obeys a Curie-Weiss law of the form k=k,+C/(7—6), with a 
fairly small contribution from kj. From the straight lines, average 
values of the constants in this equation are: kj=4; C=85 000°K ; 
6= — 80°C. 

Lead metatantalate is probably anti-ferroelectric with a transition 
temperature at about 70°c. 


4.3. Barium Metaniobate, BaNb,O, 


The normal firing technique as described previously was not suitable 
for barium niobate : although it fuses at 1350°c, approximately, from 
about 1220°c upwards a ceramic disc slowly crystallizes into a matted 
mass of needle-shaped crystals, giving a very porous and fragile specimen 
quite unsuitable for investigation. The discs were therefore fired at 
only 1200°c, and for 48 hours at that temperature since the normal 
24 hours was insufficiently long to give good sintering. The growth of 
microcrystals at 1200°c was quite slow, and only a thin surface covering 
was produced, which was removed during polishing. 

The k/T relation for this material is given in fig. 2. It will be noticed 
that all curves peak at about 70°c. The D/E relation shows a definite 
saturation effect, both ends of the loop flattening towards the horizontal. 
The degree of saturation is less than that found in aluminium niobate 
(see below), but much greater than that found in manganese tantalate. 

The &/7 curves suggest the transition temperature is about 70°c. 
The peak position is difficult to observe at 50 c/s: above 70°c the 
increasing power-factor is causing an increase in the apparent value of &, 
thus obscuring the fall in true value. Saturation can still be observed 
in the D/E loop at 120°c, for an applied field of 77 kv/cm. It may be 
that at this temperature a sufficiently strong field will pull barium 
niobate into the ferroelectric state. However, the loss at 120°c for this 
ceramic specimen would be sufficient to mask the central linear portion 
if it exists, leaving only the saturation effects to show. 
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For frequencies up to 200 c/s the peak was very flat and ill-determined, 
as at 50 c/s. At 250 c/s the peak was clearly below 100°c, and at 1000 c/s 
was fixed as in the region of 70°c. It is thus likely that the position 
of the transition temperature is independent of the applied frequency. 

The graphs of 1/k vs 7’ are curved, showing that the k) contribution 
to the whole dielectric constant is large. A reasonable correlation can 
be obtained with the measured values of & and J at a frequency of 
0-1 Mc/s for the following constants : ky=86 ; C=6060°K ; 0= —82°c. 

Barium metaniobate is thus ferroelectric, with a transition temperature 
about 70°C. 


Fig. 2 
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4.4. Aluminium Metaniobate, AINb,O, 
The 14 materials so far ri 
Is gs described have been of the o 
ice ae e general formula 
f a bs » A 18S a metal in a divalent state, and B is niobium or 
tantalum. Previous niobates and tantalates examined have been those 
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of the perovskite system A(BO3),, where A is a monovalent metal ion. 
It would thus seem logical to continue the sequence to materials of 
general formula A(BOs)s, where A is now a trivalent metal ion. Ferric 
iron would be an obvious choice in view of the ferrous materials so far 
investigated, but the results obtained from ferrous niobate and tantalate 
suggested that the corresponding ferric compounds would be equally 
conducting, and hence give no definite information. The most stable 
trivalent metal readily available is aluminium, and therefore the niobate 
and tantalate of this were prepared. 


Fig. 3 
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The D/E relation for the niobate was a thin loop curving gently towards 
the horizontal at each end and resembled that given by barium titanate 
ceramic when just below its Curie point. Due to increasing ee ee 
increasing temperature no transition temperature could be found from 

; measurements. 
“5 eal nf. curves, together with the 50 c/s curve, are phe an ae 
and show quite clearly the transition to occur at about a 0; . or 
the r.f. values of k, graphs of 1/k vs T have been plotted. “aS the ne of 
is at 220°c there are not many points above this ee aa 1rough 
which to draw the straight lines indicative of a Curie-Weiss ee : ee, 
the last four points of each set do lie on straight lines, and these have been 
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extrapolated to cut the temperature axis. Average values of the constants 
in the relation k=k)+C/(7'—6) for the lines are: ky=3, C=50 000°K ; 
?= — 80°C. 

Aluminium metaniobate is ferroelectric with a transition temperature 
about 220°c. 


§ 5. CONCLUSION 


Fifteen new metaniobates and metatantalates have been prepared 
and their dielectric properties examined in the ceramic form. Possibly 
three new ferroelectric materials—manganese metatantalate (very 
weakly), barium metaniobate, and aluminium metaniobate—have been 
discovered, and one new antiferroelectric—lead metatantalate. Strontium 
and barium metatantalates may be of interest at temperatures below 
— 183°C. 

The research described here was carried out between October 1953 
and September 1956. Since this latter date several Russian papers 
have appeared (e.g. Isupov 1957) dealing with metaniobates and meta- 
tantalates of divalent metals. Generally results are in good agreement 
with those given here, but there are two cases of severe disagreement. 


(i) PbTa,O, is reported as ferroelectric, with a transition temperature 
at 240°c ; k=310—400 at 1 ke/s and 25°c. 

(ii) BaNb,O, is reported as non-ferroelectric, with k=40 at 25°c and 
1 ke/s, and having a negative k/7T' slope. 


It must be emphasized that all results given in this paper are for 
materials prepared under the conditions quoted ; polymorphism, with 
the various forms having quite different dielectric properties may well 
exist in several cases if production conditions are varied, as it does for 
BaTiO, and PbNb,O,. An independently produced supply of barium 
metaniobate ceramic discs gave properties and x-ray diffraction patterns 
exactly similar to those of the authors’ own discs. 

No electrical or optical results are available for single crystal specimens 
due to the difficulty of growth of large enough crystals. Such work is 
needed especially on barium metaniobate. 

Little structural work has been done on any of these materials, again 
due to lack of single crystals and to the cell complexity. The iron and 
manganese compounds correspond to the mineral structures, all being 
the orthorhombic columbite cell except for iron tantalate which is the 
tetragonal trirutile dimorph. A complete structure analysis of any of 
these compounds is extremely difficult because of the low symmetry 
and large number of atoms per unit cell: Goodman (1953) describes lead 
niobate as having an orthorhombic cell of approximately 25 A x 2547 A, 
with 40 molecules per cell. Francombe (1956) reports for the orthorhombic 
polymorph 17-51 A x 18-81 Ax 7-73 4, with 20 molecules per unit cell ; 
the large sides are both approximately 1/./2 those measured by Goodman, 
thus giving a cell of half the volume. 
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Single crystals of barium metaniobate of size suitable for x-ray analysis 
were prepared. The unit cell is orthorhombic (space group either 
Ci2-A2,am, or C}$-Ama 2), of size 10-274x12:164 7-784, having 
eight molecules per cell. A possible structure has been advanced from 
packing considerations (Coates 1956). 
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Decay of Lithium-7 Hypernucleus 


By P. H. FowLer 
H. H. Wills Physical Laboratory, University of Bristol 


[Received August 7, 1958] 


HYPERFRAGMENTS are so infrequently produced in nuclear reactions, 
that to the present time few cases of those with Z>3 which offer an 
unambiguous interpretation have been found. This paper describes 
an event found in photographic emulsion which by good fortune can be 
interpreted with considerable confidence as the decay of ,Li’. It gives a 
value of 5-2 +0-5 Mev for the binding energy of the A°. In addition, it 
is almost certain that the A° hyperon was produced directly from the 
capture of the K~-meson by O1. 

The event was found in a stack of 600 » Ilford K5 emulsion exposed to 
a separated beam of K~ particles at the Bevatron, Berkeley. A photo- 
micrograph is given. The K~ adsorption at rest gave rise to three short 
black tracks (two of which are superimposed in the photograph, Pl. 103), 
one track of a lightly ionizing particle, and the hyperfragment, which 
decayed apparently at rest after travelling 32-54. The hyperfragment 
decay resulted in a negative 7-meson and two black tracks, which were 
coplanar to within the limits of error; the plane of decay was inclined at 
only 14° to the emulsion plane, so that accurate measurements could be 
made. 

To obtain the best data on the angles of dip with respect to the emulsion 
plane, the dip measurements were made with a long working distance 
objective of numerical aperture equal to 0-95, and with the emulsion 
swollen with water so that it was approximately 1000 yu thick. 

The table gives the data obtained for the hyperfragment decay. 


Energy Momentum 
aie Angle in as 8. 
Track | Range pu sorimonsplane ae : aie rae ee 
Tle 
7 10060 0+0-3° 24-1 + 0-45 85-5 + | 
A Het 33 144-541 6:9 Oz) 227 + 1:5 
B He? 23 307+1 4:94 0-1 165 +2 


The directions of emission of the three tracks were coplanar to within 
the limits of error determined by multiple scattering and grain size; this 
suggests analysis as decay at rest. Using the measured range of the 


Correspondence 1461 


7m -meson, one finds that tracks A and B should have momenta of 226 + 16 
and 164 + 16 Mev/c respectively. If track A is due to He! its range should 
therefore lie between 27 and 39 por if He®, 48 and 74 1; any other assign- 
ment gives an even poorer fit. We conclude therefore that track A is 
very probably produced by an «-particle, and that its momentum is 
therefore 227+ 1-5 Mev/c. Using this estimate of momentum for A we 
obtained a new estimate for the momentum of track B, of 165 + 2 mev/c, 
so that if it is a He? nucleus, its range should lie between 22 and 24 pts 
which it does; no other particle has a range close to 23 » with this value 
of momentum. The whole event is thus consistent with the decay of a 
ali’ hyperfragment as given below: 


ali’ > a + Het + He® + 35-9 + 0-5 Mev. 


The best estimate for the Q value for the reaction is obtained from the 
measured ranges, and the error is due almost entirely to the straggling 
of the 7-meson. This Q value gives a figure of 5-2+0-5 mev for the 
binding energy of the A° in ,Li’, taking 37-2 Mev for the Q value for decay 
of a few free A°. 

It is of interest to note that the other examples of , Li? hyperfragments 
decay according to the schemes: 

(a) ,Li’+>727-+Be’; 
or 

(b) —>a-+H!+ Li’. 
Identification of the first mode of decay is somewhat uncertain. For 
example, in the E.F.I.N.S—N.U. collaboration (private communication) 
eight examples of collinear decays were found that could be attributed to 
this mode. The range of the recoil averaging 1-8 1, however, does not 
permit an unambiguous assignment. A further important point here is 
that in an appreciable fraction of such decays, Be’ might be expected to be 
produced in its first excited state at 430 kev; decay in this manner would 
not be resolvable from decay to the ground state, owing to straggling of 
the z-meson. A number of examples of the second mode of decay have 
been put forward (Ammar ef al. 1958 and Limentani et al. 1958). 

But there is difficulty in distinguishing this mode of decay from that 
of ,Li8>n-+H1+Li’, the range of the recoil of Li’ being ~85% 
of Li® for the same value of momentum—at values of momentum 
~ 200 mev/c, and an even greater fraction at lower values of momentum, 
where separation by range-momentum is not possible. There is no 
assistance to be obtained from the @ values of the two reactions. If any 
event gives an acceptable vaiue for the energy release when analysed 
according to (b), it will give a value ~ 0-5 Mev lower if the recoil is taken as 
Li? in a typical case; this will also be acceptable. In decays of type (6), 
once again the possibility of the Li® being produced in an excited state 

ne in mind. cts 
Tn aaa described in this paper is significant in that it 1s a clear 
example of the decay of , Li’, and it gives an accurate value for the binding 
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energy of the A°, free from doubt as to identification of the products and. 
the possibility of the fragments being formed in any excited state. The 
value of 5:2 + 0-5 Mev obtained is in fact in excellent agreement with the 
estimates from those decaying according to either reactions (a) and (6) 
and of an accuracy comparable to the average of all of them. Thus it is 
almost certain that the collinear decays have been correctly interpreted 
as decays of ,Li?, but the question as to what fraction of the Be’ recoils. 
are produced in the ground state is left open. 


Tur K- ADSORPTION 


Assuming the nature of the particles producing the three short black 
tracks it is possible to interpret the event as a K~ capture by O18 involving 
no neutral particles, and one obtains an excellent balance of energy and 
momentum. 

The scheme proposed is: 

K- +016 -> 7- + 2He*+ H1+4 ,Li’. 

The resultant momentum of the four heavily ionizing tracks has a 
magnitude of 230 + 10 Mev/c, and its direction lies very close to the observed 
light track; in fact the transverse momentum is only 8+ 10 Mev/c. It is 
therefore very plausible that this track is that of the only other particle 
involved in the reaction, and that the reaction has been correctly inter- 
preted. The final check that can be made is from the total energy release, 
which is 499 + 9 Mev using this scheme, in excellent agreement with the K— 
rest mass—494 Mev. 

This excellent balance of momentum and energy obtained when one 
uses the reaction above, gives one confidence that it was very probably 
correct. Other schemes involving disintegrations of oxygen or heavy 
elements with emission of neutrons are unlikely to give such a good fit 
when analysed according to this reaction. 

If one accepts this reaction, two further points can be made. First, 
we have further evidence of the correct identity of the hyperfragment, for 
no other hyperfragment having a momentum of 416 + 6 Mev/c, as required 
from the dynamics of the parent star, has a range of 32:5 «. Second, we 
can make a better estimate of the kinetic energy of the 7, using the K- 
mass; it is 124+0-7 Mev. This value is very high, and means that either 
the A° was made directly in the interaction of the A~ meson with a nucleon, 
or that if indeed a & hyperon was formed, it was formed in a bound state 
of ~30 Mev. The former supposition seems the more likely. 
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Heat Transfer between Copper and Liquid Helium II 


By N. J. Brow} and D. V. Osporne 
Department of Natural Philosophy, University of St. Andrews, Fife 


[Received August 6, 1958] 


WHEN heat flows from a solid into liquid helium II, there exists at the 
interface a temperature discontinuity, the magnitude of which is 
proportional to the heat flow density. The constant of proportionality 
(usually expressed in deg watt-!em?) has been called the boundary 
resistance or Kapitza resistance, K, after its discoverer (Kapitza 1941), 
and it has been measured subsequently by White ef al. (1953), Fairbank 
and Wilks (1955) at temperatures below 1°K, and Dransfeld and Wilks 
(1957) at high pressures. All these measurements have been carried out 
with the aid of a constant heat current resulting in a steady temperature 
difference ; we have thought it of interest to determine the same quantity 
using an alternating heat flow. 


Copper foil 


The coupled resonators. 


The apparatus (see fig. 1) consists of two accurately similar second -sound 
resonators separated by a thin copper foil. Each resonator consists of a 
hollow cylinder of Perspex, 11mm long and 6mm inside diameter, 
closed at one end by a glass plate and at the other end by the copper foil. 
Second sound waves are generated in one resonator by means of an 
aquadag heater painted on the glass plate, and the second sound transmitted 
through the foil can be detected by an aquadag resistance thermometer 
painted on the other glass plate. The copper foil itself is so rconty. on 
(0-05 mm) and of sufficiently high thermal diffusivity (~ 10*cm*sec™) to 


: . 
+ Now at Chapel-cross Works, United Kingdom Atomic Energy Authority, 
Annan, Dumfriesshire. 
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offer no obstacle to the passage of thermal waves at the frequencies used 
(500 c/s to 2ke/s). The Kapitza resistances on the two sides of the foil 
are, however, so great that the coupling between the two resonators 1s In 
fact relatively weak, and a measurement of the signal transmitted to the 
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second resonator serves to determine the value of K. In a blank 
experiment, using a Perspex disc 9:5 mm thick in place of the copper foil, 
it has been established that the accidental coupling between the 
resonators is negligible compared with that due to the transmission of 
heat through the foil. 

The determination of K in this way requires measurements of the phase 
and amplitude of the received signal as a function of frequency in the 
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vicinity of resonance. Since, moreover, it is not possible to measure 
absolute temperature amplitudes very reliably with aquadag thermometers, 
it is necessary to repeat the experiment with the same thermometer and 
heater but with the copper foil removed; such an experiment serves 
effectively to calibrate the thermometer at the frequencies used. The 
final results presented here have been derived from a series of four runs, 
two with the foil and two without, arranged so as to verify that the 
properties of the aquadag remained reasonably constant throughout. 

The signal received in the second resonator depends not only upon 
the coupling, but also upon the attenuation present in each resonator and 
upon the difference in length, Al, between the two resonators. Al was 
small, about 545, and could not easily be measured accurately. The 
attenuation was therefore deliberately made appreciable by the use of a 
fairly narrow tube so that Al became a negligible quantity in the final 
calculations. 

The signals are small, and a phase-sensitive detector is therefore used. 
Such a device facilitates the measurement of phase, and at the same time 
has a very low noise level on account of its narrow bandwidth; a signal 
corresponding to 4udeg can be detected. The amplitude has been kept 
sufficiently small to avoid non-linear effects; the received signal is in all 
cases sinusoidal and its amplitude is proportional to the input power. 
Measurements have been made for two different modes, the fundamental 
for which each resonator is half a wavelength, and the second harmonic 
for which each is a whole wavelength. 

The results, together with those of previous workers, are shown in 
fig. 2. The accuracy of the present values of K is about + 20%, and there 
is satisfactory agreement between values derived from measurements on 
the two different modes. The dependence upon temperature agrees with 
earlier work, but the actual magnitude of 1/K is somewhat higher. The 
theoretical curve of Khalatnikov (1952) is also shown for comparison. 
The disagreement between the earlier experiments themselves discourages 
us from laying too much stress on the absolute magnitude of any of the 
results, but if it should prove that the resistance measured with an 
alternating heat flow is genuinely different from the value obtained in a 
steady state experiment, the discrepancy may be associated with the 
turbulence present in a steady heat flow, but probably absent from second 
sound. This would imply that the mechanism of the Kapitza resistance 
is to be sought in the liquid itself as well as at the interface, and that the 
theory of Khalatnikov is therefore incomplete. The matter can be 
approached experimentally by studying second sound in the presence of a 
steady heat flow, though the apparatus would take on a more complicated 
form than that shown in fig. I. . ’ 

This work is published at this stage because It has had to be discontinued 
for the time being with no immediate prospect oF its resumption. One of 
us (N.J.B.) is indebted to the Department of Scientific and Industrial 
Research for a grant during the tenure of which this work was carried out. 

5F 

P.M, 
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Electron Distribution in Transition Metals 


By W. Hume-RoTHEery 
Department of Metallurgy, Oxford University 


P. J. Brown, J. B. Forsytu and W. H. Taytor 
Crystallographic Laboratory, Cavendish Laboratory, Cambridge University 


[Received September 19, 1958] 


In a recent paper, Weiss and De Marco (1958) have used refined x-ray 
methods to determine the absolute scattering factors for crystals of 
transition metals of the First Long Period. The results are then analysed 
on the assumption that the scattering results from the sum of the effects 
of (1) the argon cores of the atoms concerned, (2) the electrons in 3d-like 
states, and (3) the electrons in 4s and 4p-like states. At the Bragg angles 
concerned, the effect of (3) is negligible, so that the number of 3d electrons 
can be obtained by subtracting the known effect of (1) from the experimen- 
tally determined structure factors. In this way, the following results are 
obtained for the numbers of 3d electrons per atom: 


Cu face-centred cubic 9-8 
Ni face-centred cubic 9:7 
Co hexagonal close-packed 8-4 
Fe body -centred cubic 2°3 
Cr body-centred cubic 0-2 


In several recent theoretical papers, these numerical values have been 
taken as well-established experimental facts, and used to justify various 
bonding schemes or hypotheses. It is the object of the present note to 
suggest that great caution is needed before the results of Weiss and 
De Marco are used in this way, and that it is doubtful whether anything 
except a complete 3-dimensional electron—density map is of value for 
comparison with theoretical models. 

The experimental work of Weiss and De Marco was done with great 
care and accuracy. The interpretation involving the subtraction of the 
effects of the argon cores is entirely reasonable, and may be regarded as 
justified by the results obtained for Cu whose Cu+ ion contains 10 (3d) 
electrons and has a nearly spherically symmetrical electron cloud. For 
the remaining elements the interpretation is much less satisfactory owing 
to the fact that only a few reflections of low indices were examined, The 
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critical sentence in Weiss and De Marco’s paper is: “assuming that the 
residual scattering factors (after subtracting the ‘argon core’) for the 
lowest angle Bragg peak are entirely due to electrons with 3d-like radial 
charge densities, then that peak alone is sufficient to establish the number 
of 3d electrons’. Actually, the measurement from a given reflection 
refers only to the corresponding direction in, the crystal, and the method 
is, therefore, assuming that the electron density is spherically symmetrical. 
This is, however, the very thing which is in dispute, and on which the 
different theories differ so profoundly. It will be noted that, in the table 
given above, the apparent discontinuity between Fe and Co is also that 
between body-centred cubic metals (Cr, Fe) and close-packed metals 
(Co, Ni, Cu). The measurements refer to (110) reflections in the former, 
and to (100) (Co) and 111 (Ni, Cu) reflections for the latter, and it may be 
this difference which is revealed in the table given by Weiss and De Marco, 
rather than any discontinuity in the electronic configurations. 

It should be emphasized that the elements Fe, Co, and Ni form a small 
group whose physical properties show such regular sequences that an 
abrupt break in the electronic structure between Fe and Co appears 
improbable. It may also be noted that, according to Deslattes (1958), 
the electronic structures of Weiss and De Marco are not in agreement 
with the details of the K-series x-ray spectra. The treatment of the 
3d electrons as forming an independent group (and not as being in hybrid 
(4s 3d 4p) states) is of course equivalent to assuming what is in dispute 
between the different theories of the transition metals. It seems, therefore, 
that, although the experimental observations of Weiss and De Marco 
are accurate for the directions concerned, a complete three-dimensional 
electron—density map will be required before the electron distribution 
can be analysed with sufficient accuracy to be compared with the 
requirements of the different theories. . 

It should also be pointed out that Pauling (1953) was the first to put 
forward the view, which has recently been revived by later authors, that 
the saturation moment of 2:2 wy for Fe results from atoms with the (3d)° 
atomic configuration, together with a slight contribution from interaction 
between the core and the conduction band. 
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Nucleation of Cracks by the Intersection of Twins in a-Iron 


By D. Huu 
Metallurgy Division, Atomic Energy Research Establishment, Harwell 


[Received August 8, 1958] 


As part of a study of the nucleation of fracture in «-iron, the experiments 
described below have been made on a 3% silicon iron which had been 
produced by vacuum melting. Single crystal specimens 0-8 in. long and 
0:15 in. wide were cut from a large-grained polycrystalline sheet, 0-002 in. 
thick, so that the cube edges of the body-centred cubic structure were 
parallel to the edges of the specimen. The cube face was in the plane of 
the specimen. Before straining, the specimens were electropolished. 

The crystals were strained in tension along the [001] axis and in this 
way the maximum shear stress operated on four {011} slip planes. When 
tested at 78°K the crystals broke without measurable deformation taking 
place. Examination of the specimens after fracture at this temperature 
revealed that {112} twins had formed. No slip was observed. The 
final fracture consisted of {001} cleavage faces and irregular cracks along 
prominent twins. Cracks were not limited to the main fracture and as 
many as ten smaller cracks were observed in other parts of the specimen. 
It was clear that both the main fracture and the subsidiary cracks had 
formed at the intersection of twins. Two methods of nucleation were 
observed and these are illustrated in the photographs, figs. 1 and 2, Pl. 104. 
In the first case, fig. 1, the crack has formed along an {001} cleavage plane 
at the intersection two {112} twins; part of the crack has developed along 
one of the twins in a zig-zag manner, but retaining {001} cleavage facets. 
In the second, fig. 2, the crack has formed at the junction of one of the 
{112} twins with the other. In both cases the {001} cleavage crack has 
grown in both directions from the intersection. 

Initiation of cracks by intersecting twins has previously been observed 
in body-centred cubic molybdenum (Cahn 1955) and in close-packed 
hexagonal zinc (Bell and Cahn 1958), but has not apparently been 
reported for «iron. It has long been known that in tests at low tempera- 
tures on «-iron there is a connection between twinning and fracture, and 
two conflicting views have been taken, firstly that twins were nucleated 
by the high stress concentrations associated with fracture, and secondly 
that the formation of twins actually initiated the fracture as observed in 
the present experiments. It is probable that both effects occur. The 
mechanism by which intersecting twins produce a crack is as yet not 
known, but it may well be connected with Rose’s channels (Rose 1868). 
Miigge (1922) estimated that there was a 50% increase in volume at the 
intersection of two twins on reciprocal planes, e.g. (112) and (112) as a 
result of Rose’s channels in «-iron. He suggested that since this is 
unlikely cracks formed instead, However, Hall (1954) believes that the 
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strain set up is much more likely to be relieved by local slip instead of 
cracking. At very low temperatures and in structures with only a small 
number of slip systems, e.g. close-packed hexagonal, slip is difficult and 
such relief of strain will be less likely. 

Barrett et al. (1937) reported that the stresses to initiate twinning and 
slip in iron and iron-silicon alloys depend on temperature and composition. 
In the material used for the present experiments the stress for twinning is 
less than that for slip at 78°K, and at this temperature the twins interact 
immediately to form cracks. It is hoped to make similar tests at higher 
temperatures, where slip is preferred to twinning, to determine whether 
cracks can be formed at the intersection of slip planes, as postulated 
recently by Cottrell (1958). 
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The Absorption of Sound in Dilute Solutions of Helium-3 in 
Liquid Helium II 


By G. O. Harprne and J. WILKS 
Clarendon Laboratory, Oxford 


[Received July 31, 1958] 


We have measured the coefficient of absorption of sound at a frequency of 
14-7 Me/s in liquid helium IH containing small concentrations of “He. A 
pulse method with a variable path length was used as has been described 
previously (Dransfeld et al. 1958). The results for pure helium a and for 
concentrations of He up to 5-2% are shown in the figure. It is obvious 
that below 1-8°K the presence of the *He considerably reduces the 
ie eee nicnt to discuss these results in the form of reduced es 
equal to the observed absorption divided by the factor 27f?/pc?. ic 
necessary values of the density p and velocity of sound ¢ may be ae 
from the measurements of Dash and Taylor (1957) and Flicker and Atkins 
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(1957) respectively.) The absorption of sound in 3He—*He mixtures has 
been treated theoretically by Khalatnikov (1952a) who showed that for 
not too high frequencies the reduced absorption is given by 

4 De® ($) 


on tnt fat ae) UF 


where 7, and ¢, are the coefficients of first and second viscosity, Disa 
diffusion constant, m, the mass of a ?He atom, and « the mass concentration 
of the 3He. ‘The first two terms are similar in form to those obtained for 
pure helium IT (Khalatnikov 1950, 1952b), while the last term is an 
additional one arising from the presence of the He. Fluctuations in the 
density of the liquid set up by the passage of the sound wave produce 
concentration gradients of *He atoms, and diffusion then results in 
irreversibility and absorption. 


o% 
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The coefficient of absorption of sound of 14:7 Mc/s in liquid helium IT containing 
’He. The figures on the curves indicate the molar fraction of 3He. i 
Abscissae: Temperature °K. 
Ordinate: Absorption (em~). 


The reduced absorption arising from diffusion can be calculated by 
estimating dp/de from the measurements of Dash and Taylor, and using 
the value of D found by Khalatnikov and Zharkov (1957) from the results 
of Beenakker et al. (1952); this contribution to the absorption is always 
small. The coefficient of first viscosity is also known from the experiments 
of Dash and Taylor; this term accounts for only 1 or 2 tenths of the total 
absorption. Thus the observed decrease in absorption can only be 
accounted for if the coefficient of second viscosity is considerably reduced 
by the presence of the ?He atoms. 

Second viscosity arises when the thermodynamic properties of a fluid 
do not immediately take their equilibrium values consequent on an 
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instantaneous change of pressure (see for example Wilks 1958). In liquid 
helium I, the numbers of thermal excitations do not immediately follow 
changes in pressure, and the approach to equilibrium is described by 
relaxation times characteristic of phonon-phonon and phonon-roton 
collisions. It appears that the presence of the 3He atoms acts as a pertur- 
bation which increases the number of collisions which restore equilibrium, 
thus reducing the relaxation time and hence the absorption. 

The variation of the absorption with concentration suggests that the 
number of additional collisions is proportional to the number of ?He atoms 
present. Khalatnikov’s treatment involves two relaxation times and his 
expression for the second viscosity is rather complicated, but for our 
purpose it is possible to introduce a simplification. Thus above 1-2°K the 
viscosity may be described tolerably well in terms of one relaxation time 
7 which varies with temperature in a similar way to the times given by 
Khalatnikov. From the general theory of relaxation processes 


Ce 


wr 
“Tot 
where w is the angular frequency of the sound. In the present experiments, 
at all temperatures above 1-2°, w<1/7 so that ¢, ocr and 1/¢, is proportional 
to the collision rate between the excitations. It follows that if (£,), is the 
second viscosity in pure helium II and (,). is the viscosity for a concen- 
tration e of ?He atoms, then 


(1/2). —(1/La)o] cn. 


where J, is the additional number of collisions due to the presence of the 
3He. An examination of the data shows that to a fair approximation 


[(1/€2).—(1/Ee)o] ace 


and one thus obtains the previously mentioned conclusion that the 
increase in collision rate is roughly proportional to the number of *He 
atoms. We are at present extending this work to lower temperatures. 
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The Strength of Aluminium Silver Alloys 


By AntTHony KELLY 


Department of Metallurgy, Northwestern University, 
Evanston, [linois 


| Received September 26, 1958] 


SINGLE crystals of aluminium—6 atm °% silver quenched from the 6-phase 
field and aged at room temperature show a critical resolved shear stress at 
0°c of 7-8 kg mm". Over the temperature range 77° to 373°K this value 
changes by not more than 10% (Lassila 1958). The value of the critical 
resolved shear stress is two orders of magnitude greater than that of pure 
aluminium single crystals. This increase in strength due to the addition of 
silver is intriguing since lattice strains are likely to be small in this system 
as the Goldschmidt radii of silver and aluminium differ by only 0-7% 
and solution of silver in aluminium produces no change in the X-ray 
lattice parameter (Axon and Hume-Rothery 148). 

Haessner and Schreiber (1956) have shown that amounts of silver less 
than the solubility limit (0-16 atm ° at room temperature) produce no 
change in the plastic properties of aluminium single crystals. Lassila 
in this laboratory has found that aluminium specimens containing 6 
atm °% silver have an initial flow stress approaching that of pure aluminium, 
if tested below room temperature immediately after quenching from the 
solid solution range. ‘Solid solution strengthening’ thus seems to be 
small. 

X-ray studies of these alloys (Walker and Guinier 1953, Gerold 1955) 
show that on ageing at room temperature the alloy contains zones con- 
sisting of spherical clusters of atoms containing a high concentration of 
silver surrounded by a spherical shell rich in aluminium atoms. A large 
volume fraction of the material is in the form of zones which are completely 
coherent with the matrix. Gerold finds that the hardness of the alloys 
depends only on the volume fraction of the zones and not on the zone radius 
or number of zones per unit volume. X-ray examination, of the crystals 
studied here shows the cluster radius is c.17 A and the radius of the zone 
c. SDA, 

Although it cannot be proved conclusively that lattice strains play no 
part in determining the strength of these alloys it seems reasonable to look 
for alternative possibilities. 

If we view the silver rich cluster and aluminium rich shell of each zone 
as forming a region in the alloy containing an equilibrium number of 
Al-Al, Al-Ag and Ag—Ag nearest neighbours, characteristic of a certain 
ageing temperature, the passage of a dislocation through a zone will alter 
the number of nearest neighbours of a given type across the slip plane. 
This will raise the free energy of the alloy and hence a higher stress will be 
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needed to move a dislocation through the alloy containing zones than 
through the random solid solution. This type of strengthening mechanism 
was suggested originally by Griffith (1920) and later by Fisher (1954) 
to account for the strength of alloys containing short range order. 

Since the distance between the centres of zones, c.100 4, is about twice 
their radius a dislocation cannot bow out between them but must pass 
through them approximately as a rigid line. The dependence of the 
strength solely on the volume fraction of the zones can then be explained 
in the following way. The average separation of the zones intersecting a 
given slip plane will vary as (R?/f)"2 where R is the average radius of a 
zone and f the volume fraction of zones. The average length of the 
intersections of zones by a dislocation will vary as R and equating work 
done on the slip plane by the applied stress to the energy of the interface 
formed, a formula for the flow stress (a) of the type 


1/2, 
oa a 


is obtained. Here y is the energy of the interface produced when a zone is 
sheared by one interatomic distance, b is the Burgers vector and a is a 
constant which for spherical particles equals 2/\/7. Thus the flow 
stress is independent of zone size and hardening of this type should depend 
little on temperature at low temperatures, as is observed. 

An estimate of y is difficult. It will depend on the part of the zone 
through which a dislocation passes, on the width of the dislocation, and on 
the angle between a dislocation line and its Burgers vector. However, 
a rough value for a particular model can be obtained from the heat of 
reversion of the alloy (Kelly and Fine 1957). According to the x-ray 
results 10% of the total number of silver atoms in the alloy leave the 
clusters during reversion. Koster and Schell’s (1952) value of the heat 
of reversion (1-04 cal. gm~') can then be expressed as 0-24 ev per atom of 
silver leaving the clusters. To a rough approximation we can write for 
the shearing of the clusters 


_1AE 
1-3 
where AZ is the heat of reversion per atom of silver leaving the clusters. 
Whence y ~ 150 ergsem~*. If we then calculate the flow stress according to 
eqn. (1) considering only the clusters at the centres of the zones and, again 
following the x-ray results, that the number of atoms in the clusters is 
one-third the total number of silver atoms we can take f= 0-3 x 0:06 = 0-018, 
and we find o~ 7:0 kg mm. 

This estimate shows that this type of hardening should be of importance 
in these alloys and perhaps in other age-hardening systems where lattice 
strains are expected to be small, e.g. Fe-Cr and Cu-Co. Further, the 
flow stress should depend solely on the volume fraction of precipitate 
provided the spacing of precipitates is such that a dislocation passes 
through them as a rigid line. 
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REVIEWS OF BOOKS 


Logic Machines and Diagrams. By Martin Garpner. (McGraw-Hill Book 
Od.; Ine.) [Pp.152.]° $4.75; 


THE syllogism dominated logic for two thousand years. The formidable 
weapon which finally brought about its downfall was symbolic logic, developed 
by Boole and others. Logical diagrams and machines also played their part 
by showing, in a non-mathematical and entertaining way, what the new logic 
was allabout. That the subject can still be entertaining is shown in this book. 
Perhaps the machines described now strike us as complicated and creaking 
contrivances for what whey were able to do. However, this part of the subject 
is far from dead; enthusiasts still make logical machines, although they use 
relays of transistors instead of levers and wires. Formerly, problems that 
could be resolved by symbolic logic occurred only rarely in practice; now, such 
problems arise naturally in the design of complex switching circuits and it may 
be that a real demand for logical machines will be created. It is more likely, 
however, that any such demand will be met by the use of general purpose 
digital computers, suitably programmed, than by the construction of special 
machines. 

The author of the book is described on the dust jacket as one of America’s 
best-known writers in the field of science. It is not a popular book in the 
ordinary sense of the term; it would be just the thing, however, for readers of 
The Philosophical Magazine to take away on holiday. M. V. WILKES. 


Quantum Mechanics: non-relativistic theory. By L. D. Lanpau and E. M. 
LirsuiTz. Translated from the Russian by J. B. Sykes and J. S. Bell. 
(London : Pergamon Press.) [Pp.515.] 80s. 


Amongst the numerous treatises on Quantum Theory, this is surely one of the 
most useful. It is not, apparently, meant as a first text. The reader would 
be supposed to know the elements of wave mechanics, and be used to thinking 
in terms of wave-packets and probability. But here he will find a detailed 
account of the subject (up to but not including the Dirac Equation) with elegant 
expositions of all the general principles, valuable theorems, and powerful 
techniques, many of which are only to be found in scattered papers and books. 
It is not clear whether the lucidity of the style stems from the Russian 
original, or from the translators; it could scarcely be bettered. 

The tone is essentially pragmatic and eclectic. Purely mathematical questions 
(What axioms define Hilbert space ?) and philosophical doubts (Is quantum 
theory compatible with determinism ?) are entirely avoided. Von Neumann 
is only mentioned for a theorem on the intersection of molecular terms ! The 
subject is not forced into a restricted formal mould. Every technique is used 
where it is most powerful, so that the amount of mathematical manipulation 
is cut to a minimum. There is a happy balance between wave functions and 
matrix operators—between the Schrédinger and Heisenberg representations. 

What is striking is the completeness with which the whole subject was 
developed between 1926 and, say, 1933. Almost all the references (which 
might perhaps have been given more chapter and verse than an author and 
a date) are from those heroic years. Indeed, this book, basically written 
in 1947, is a little old-fashioned in its choice of applications of the theory. 
We should now expect to include, along with the theory of atoms and molecules, 
the structure of the nucleus, and the electronic properties of solids. The 
many-body problem and the role of correlation and exchange energies, would 
also naturally demand more space. Perhaps these will be dealt with in other 
volumes of this admirable series by these distinguished authors. J. M. Z. 


[The Editors do not hold themselves responsible for the views 
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(b) 


Cc) 
Electron diffraction patterns obtained from a specimen of aluminium—chromium 
by transmission. (a) Aluminium film 200 4 thick; (6) chromium film 
150 A thick; (c) compound film. 


Cb) 


Electron diffraction pattern obtained from a very thin specimen of aluminium— 
chromium. (a) Aluminium film 60 A thick; (6) chromium film 30 4 
thick; (¢) compound film. 
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2 x 108 cycles ( x 354). 
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Strain-time oscillograms. 

(a) Recorded from central gauge. Height of drop: 28in.; Strains corre- 
sponding to calibration lines: 0, 0-92 and 1-29 x 10-3. (b) Recorded from 
series combination of central gauge and gauge mounted in a radial direc- 
tion at 3in. radius. Height of drop: 28 in.; Strains corresponding to 
calibration lines: 0, 0-92 and 1-29x 10-8. (c) Recorded from gauge at 
centre of steel plate (0-330 in.). Impact velocity of steel sphere: 243 
in./sec; Strains corresponding to calibration lines: 0, 0-34, 0-68, 0-97, 
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(0) 


Typical oscilloscope records. 


(a) Spiral oscilloscope record. Positive pulse blanking trace ; negative pulse 
brightening trace. One gyration takes 2 psec. 
(b) Type 58 oscilloscope record. Sweep time 1 psec: 0-1 psec calibration 
. pulses. Eight events, polarity sequence —++——+++. 
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ORIGIN AND DECAY OF ali? HYPERNUCLEUS 


A negative A-meson is captured at rest by an oxygen nucleus, producing a 
fast m--meson, a ,Li7 hypernucleus, a proton and two low-energy 
a-particles. The tracks of the two a-particles are superimposed and 
cannot be separately distinguished in the photograph. 

The hyperfragment decays at rest according to the relation 
ali? >a + Het + He3, 
the He* producing the longer of the two dense tracks. 


D. HULL Phil. Mag. Ser. 8, Vol. 3, Pl. 104. 
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Crack nucleated by two intersecting twins. (x 140.) 


Fig. 2 
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Crack nucleated at the junction of one twin with another. (x 140.) 
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